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ABSTRACT 

The lowering of'dielectric constant of the dispersing medium by adding some non-electrolyte to a 
lyophobic sol shows a decrease in its stability but this being not a general phenomenon, in this paper 
the changes in the dielectric constant of the dispersing medium water has been varied from 78’54 to 45*39 
by the addition of dioxan to the hydrous ferric oxide sols, which were prepared by three different 
methods. The results show that the behaviour of these sols are not identical. While the sols of hydrous 
ferric oxide prepared by peptising the hydrous oxide either by dilute hydrochloric or acetic acid become 
unstable for all concentrations of dioxan, the sol prepared by Kreck’s method is stablised by adding 
smaller concentrations of dioxan, but becomes unstable when the concentration of dioxan increases beyond 
10%. These results indicate that the nature of hydrous ferric oxide considerably modifies the effect oh 
non-electrolytes. It is well known that the coagulum of the hydrous oxide precipitated from Kreck’s sol 
is more compact than those obtained from other two samples of the sols, where the coagulum is highly 
flocky and tends to form gel. 

INTRODUCTION 

The action of non-electrolytes on the stability of lyophobic colloids has 
largely been investigated. 1 Recently Ghosh and coworkers* have studied the effect 
of various non-electrolytes on the positively charged hydrous ferric oxide and nega- 
tively charged manganese dioxide sols. They have shown that besides the various 
factors the free surface energy of the dispersed particles also guides the stability of 
a sol in the presence of some non-electrolytes. While many factors have been 
ascribed to affect the stability of a sol, viz., the preferential adsorptive capacity of 
the peptising or coagulating ions, the amount of added electrolyte and its specific 
character, the chemical interaction between the coagulating and peptising ions, the 
dielectric constant 3 of the medium has also been emphasised to control the stability 
of a sol in the presence of a non-electrolyte. Freundlich 4 opined that the lowering 
of dielectric constant leads to the diminution of the charge on the particles ahd 
hence its stability. Keesom 5 also confirmed this view by noting stabilisation with 
the non-electrolytes which increase the dielectric constant of water. But Mukerji® 
from his cataphoretic measurements could not contribute to this view. 
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In this investigation. dioxan(diethylene dioxide) has been employed to lower 
the dielectric constant of the dispersing medium water and potassium chloride has 
been used as a coagulating electrolyte for the various sols of hydrous ferric oxide 
obtained by the different methods. It has been further observed that the effect 
of dioxan on the stability of sol is also controlled by the time of observation such 
that the effect becomes more pronounced when the time of observation is smaller. 
This had been neglected by previous workers. 

EXPERIMENTAL 

Three different samples of hydrous ferric oxide sol were prepared by peptising 
the hydrous oxide with hydrochloric and with acetic acids and by Kreck’s method. 
These sols were nomenclatured as A, B and G respectively. The sol samples ob- 
tained by one method but having varying amount of peptising electrolyte are indi- 
cated by the suffix 1, 2 and 3. Iron content of the sols samples were kept the 
same (3*6298 gm of Fe +++ /litre) and their pH was adjusted to 2*2*, 3*6 and 4*5 by 
dialysing the sols to different extents. Trie purity of sols were determined by 
estimating the ratio of iron to chloride present in the system, stability of these 
sols towards its complete coagulation by potassium chloride solution was deter- 
mined. As the time oi observation for coagulation greatly influences the precipi- 
tation value of an electrolyte, it was obtained graphically for infinite time. In 
order to avoid the variations of the values with temperature, the experiments were 
carried out in a precision theremostat. 2*0 ml of the sol was taken in each set 
of the test tubes. Different amounts of coagulation electrolyte made to 13*0 ml. 
with distilled water were taken in another set and kept in the thermostat till they 
attained the temperature of the bath. The contents of the two were mixed and 
left immersed in the thermostat undisturbed and the time of complete coagulation 
was noted at the point of visual separation of the clear liquid at the upper surface 
of the test tube. Similar studies were done in the presence of a known amount 
of the dioxan keeping the total volume the same, i.e, 15 ml. The dielectric 
constants of the mixtures of dioxan and distilled water were also measured but that 
of mixtures of electrolyte and sol could not be ascertained as measurement of the 
same in the presence of electrolytes is not ordinarily possible. 

In the following tables we are reproducing some of our observations recorded 
at (30 ± 0Tj°G. In tables 1-9 the values of reciprocals of the time of coagulation 
are given for different amounts of coagulating electrolyte with varying % of dioxan. 

TABLE I 


Sol sample A x 

(Purity of the sol Fe+++/Cl“ •= 3*6296/0*96424 = 3*7642) 


2-ON KOI 
(ml.) 


lit X 10^ in Sec. 



0% 

•Amount of Dioxan added 

5% 10% 20% 

30% 

40% 

6*5 

5*263 

8*621 

12-500 

21-28 

83-33 

100-0 

6*0 

4*484 

6-061 

10-750 

16-13 

68*39 

71-43 

5*5 

4*484 

5-063 

74-630 

15-62 

50' 00 

55-56 

50 

3*086 

4*386 

5-556 

9*524 

40*00 

45-45 

4*5 

2*165 

3*817 

4-065 

75*19 

35*71 

40-0 

4*0 

1*773 

3*460 

3*425 

6-579 

28*57 

34-48 

35 

1'418 

2-762 

2*959 

6 211 

20*83 

27-78 

3*0 

1T00 

1*912 

2.506 

5-618 

16*95 

25-0 

2*5 

0*9156 

T548 

1*869 

4-405 

13*89 

21-28 

2*0 

0*8019 

1*176 

1*393 

3690 

12*05 

18-52 


* pH of »ol sample B 3 was 2*9 as it could not be obtained of lower pH value. 
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TABLE i 
Sol sample A% 
(Purity = 1 1 .6348) 


0-2N KC1 
(ml.) 



1 /t X 10 4 in Sec. 





Amount of Dioxan added 





0% 

5% 

10% 

20% 

30% 

40% 


6‘5 

3-546 

8*772 

14*29 

43*48 

86 96 

200*0 


6'0 

2*801 

7*042 

11*11 

31*25 

71-43 

166*7 


5*5 

2-50 

5*495 

9*709 

27*03 

52*63 

142*9 


5*0 

2*053 

4-566 

75*19 

19*05 

41-67 

111*1 


4*5 

1*862 

3*472 

66*23 

15*15 

35*71 

90-91 


4*0 

1*555 

2*817 

62-50 

12*50 

28*57 

71*43 


3*5 

1-250 

2-506 

46-51 

10*00 

22*94 

55*56 


3*0 

0*9615 

2*141 

41*67 

7-937 

19-92 

44*44 


2*5 

0-8621 

1*949 

28*57 

68*03 

16*95 

37*04 


2'0 

0-6993 

T686 

22*32 

5*556 

14-08 

31*45 





TABLE 3 







Sol sample A z 







(Purity = 19 6897) 




O' IN KC1 
(ml.) 



1 [t X 10 4 in 

Sec. 






Amount of Dioxan added 




0% 

5% 

10% 

20% 

30% 

40% 


6*5 

2-618 

3*086 

3*802 

6*250 

22*22 

62*50 


6*0 

1-972 

2*506 

3*155 

5*236 

20*00 

52*63 


5*5 

1-701 

1-972 

2*667 

4*630 

17*86 

42*55 


5‘0 

1-399 

1*618 

2*123 

3*906 

14*71 

34*48 


4*5 

1-096 

1-404 

1*770 

3*378 

11*79 

2778 


4*0 

0-8475 

1-188 

1*486 

2*657 

10*00 

22*22 


35 

0-7246 

1-000 

1*233 

2-591 

8*197 

18*18 


3*0 

0-6135 

0-813 

1-101 

23*98 

7*042 

14*93 


2*5 

0-5291 

0*6757 

0-9119 

1*9610 

5*848 

13*51 


2 0 

0-4587 

0*5747 

0*7780 

1-623 

5*00 

11 76 
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irAeLfe 4 

Sol sample 


2'0N KC1 
(ml,) 



l/t X 10 4 in 

Sec. 




Amount of Dioxan added 


40% 

0% 

5% 

10% 

20 % 

30% 

6‘5 

2*193 

3 - 906 

8-772 

11-63 

22*73 

62*50 

6*0 

1-972 

3-333 

6*897 

10-0 

19*23 

52-63 

5*5 

1*637 

2*817 

5*291 

8*197 

15*87 

41*67 

5*0 

1*364 

2*392 

4*405 

6*757 

12*50 

33*33 

4*5 

1*100 

1-742 

3*333 

5*618 

10*0 

27*78 

4*0 

0*8929 

1*486 

2793 

4*274 

8*333 

22*22 

3*5 

0*7246 

1*233 

2*222 

3*472 

6*944 

18*52 

3*0 

0*5556 

1*101 

1*901 

3*012 

5*618 

15*62 

2*5 

0*4386 

0-8772 

1*585 

2*500 

4*975 

11*90 

2-0 

0*3690 

0*6944 

1*233 

1*972 

3*968 

10*0 





TABLE 5 






Sol sample B s 




0*5N KC1 
(ml.) 



1 It X 10* in 

Sec. 




Amount of Dioxan added 




0% 

5% 

10% 

20% 

30% 

40% 

6*5 

3*333 

5*587 

7*407 

15*87 

31*25 

58*82 

6*0 

2*793 

4*762 

6*250 

12*90 

25*0 

52*63 

5*5 

2*193 

3*906 

4-902 

10*99 

21*28 

45*45 

50 

1*862 

3*257 

4-237 

8-772 

16*67 

37*04 

4*5 

1*585 

2*660 

3*472 

7*042 

14*08 

31*25 

4*0 

1*321 

2*169 

2*778 

5*882 

11*49 

27*03 

3*5 

rioi 

1*759 

2*294 

5-00 

10*0 

21*28 

3*0 

0-8772 

1*399 

1*957 

3*922 

7*813 

17*54 

2*5 

0*7463 

1*143 

1*661 

3*311 

6*494 

14*08 

20 

0-6061 

1-000 

1*383 

2*801 

5*556 

11*63 


[ 152 ] 



TAfeLfc 6 
Sol sample B 3 


OTN KOI 

(ml.) 



lit X 10 4 in Sec. 



0% 

Amount of Dioxan added 

5 % 10 % 20 % 

30% 

40% 

6*5 

7*143 

WOO 

18*52 

40*0 

83*33 

250*0 

6*0 

5'556 

8*475 

15*62 

33-33 

71*43 

200*0 

5*5 

4*587 

7*042 

12*35 

27*78 

58*82 

153*8 

5*0 

3*906 

6*024 

10-26 

21*74 

50*0 

125-0 

4*5 

3*311 

5*00 

8-547 

17*86 

40*0 

100*0 

4*0 

2*710 

3*922 

6-944 

16-67 

33*33 

83*33 

3*5 

2*190 

3*311 

5*495 

12*35 

26*67 

66*67 

3*0 

1*748 

2*770 

4*587 

10*26 

22*22 

55*56 

2*5 

1-453 

2' 398 

3 966 

8-772 

18*52 

47*62 

2*0 

1*233 

1*972 

3*472 

7*937 

15-62 

40*00 




TABLE 7 






Sol sample C t 






(Purity = 5-6003) 




2‘0N KG1 



1 It X 1G 4 in Sec. 










(ml.) 


Amount of Dioxan added 




0% 

5% 

10% 

20% 

30% 

40% 

6*5 

16*78 

9*709 

12*90 

19-53 

33*33 

69*44 

6*0 

13*33 

7-813 

10*64 

16-23 

26*32 

53*19 

5*5 

10*87 

6-536 

8-621 

12-90 

21*98 

41*67 

5*0 

8*13 

5-319 

6*993 

10*64 

17*70 

34*48 

4*5 

7*246 

4*464 

5-714 

8*850 

14*71 

28*17 

4*0 

5*848 

3*597 

4-630 

6*897 

11*90 

22*42 

3*5 

4*762 

2*924 

3*636 

5-882 

9*524 

20*00 

3*0 

3*802 

2*398 

2*924 

4*762 

7*463 

14*93 

2*5 

3*012 

1*946 

2-398 

3*846 

6*024 

11*90 

2*0 

2*506 

1-562 

1*972 

3*125 

5*000 

10*00 
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TABX J6 & 
Sol sample, C a 


(Purity = 17*0644) 


0'2N KCl 
(ml.) 



l/t X 10 4 in Sec. 



0% 

Amount of Dioxan added 

5 % 10 % 20 % 

30% 

40% 

6*5 

22*22 

1429 

14*93 

28*17 

37*74 

62-50 

6*0 

17*54 

10-99 

11-90 

22'22 

2985 

50*0 

5*5 

13*99 

8*85 

10*00 

17*99 

25*00 

40*82 

5*0 

10-87 

7*246 

8*333 

14*39 

21-28 

33-33 

4*5 

8-772 

5*495 

6-536 

U-76 

16 67 

27*03 

4*0 

6*667 

4*386 

5*495 

9-524 

14*08 

20*0 

3*5 

5-618 

3-401 

4-464 

7*407 

11*17 

18*52 

3*0 

4*651 

2-817 

3*597 

6*024 

9*347 

14-93 

2*5 

3*788 

2*381 

2*941 

4*808 

7*519 

11*76 

2*0 

3*155 

1*976 

2*506 

3*968 

6-250 

1000 




TABLE 9 






Sol sample C 3 





(Purity = 34-1281) 






1/1 X 10* in 

Sec. 


■ 

0*1 KCl 







(ml.) 


Amount of Dioxan added 


. 


0% 

5% 

10% 

20% 

30% 

40% 

6*5 

27-03 

12-90 

15*62 

17-86 

29-41 

71-43 

6*0 

22-22 

10-65 

12*82 

14*08 

23-26 

55-56 

5*5 

21-57 

8*696 

10-87 

11*43 

20-00 

4651 

5*0 

13-89 

7*042 

- 8*772 

9*524 

16-67 

38-46 

4*5 

11-43 

5*714 

6*944 

7-937 

13-33 

29-85 

4*0 

9-259 

4*762 

5*882 

6*667 

10-87 

22-73 

3*5 

7-519 

4-000 

4*762 

5*556 

8-850 

18-52 

3*0 

6-024 

3-077 

3*788 

4-651 

7-463 

14'81 

2*5 

4-762 

1-984 

3*067 

3-788 

6-024 

12 '05 

2*0 

3-968 

1-821 

2*500 

3*125 

5-00 

10-00 
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The coagulation value for infinite time may be taken as a measure of the 
stability of a sol. Same has been obtained by extrapolation and are recorded in 
the following table. 

TABLE 10 


Values of coagulation in ml of KCi for infinite time at 30°G for sol samples. 


%0f 

Dioxan 

used 

Ai . 

A* 


Bx 

B s 

b 3 

G, 

C 2 

C 3 

2*0N 

0-2N 

Strength of potassium chloride required 

0*1N 2*0N 0*5N 0*1N 2*0N 

0’2N 

O' IN 

0 

1-20 

1*35 

1*40 

1*45 

1*40 

1*50 

T10 

T25 

1*35 

5 

1*05 

1*20 

1*30 

1*25 

T30 

J*35 

T25 

T35 

1*50 

10 

0*90 

1*00 

1*10 

TOO 

1 05 

T20 

T15 

T30 

T4Q 

20 

0*80 

0 85 

0*90 

0*85 

0*80 

1*05 

,T05 

1*20 

T20 

30 

0*70 

0*75 

0*80 

0‘65 

0*60 

0*95 

0*80 

1*00 

TOO 

40 

0*50 

0*60 

0'65 

0*50 

0*45 

0*80 

0*60 

0*80 

T85 

TABLE 11 


Values of dielectric constant of water and dioxan mixture at 25°*C. 


Dioxan in 
ml. 

Distilled water 
in ml. 

% of Dioxan in 
Mixture 

Values of dielectric 
constant 

0-0 

15*0 

0 % 

78*5400 

0-75 

14*25 

5 % 

72*8448 

1-50 

13*50 

10% 

69-9200 

3*0 

12*0 

20% 

61*3368 

4*5 

10*5 

30% 

51*3473 

6*0 

9*0 

40% 

45*3900 

15*0 

0*0 

100% 

2*2090 


RESULTS AND DISCUSSIONS 

It is evident from tables 1-9 that with the lowering of the dielectric constant 
of the medium, the sols of the hydrous ferric oxide prepared by peptising the same 
, with hydrochloric and with acetic acid became unstable. In the case of the sol 
obtained by Kreck’s method there is first a stabilisation upto 10% of dioxan used 
and then there is a regular sensitisation, when the dielectric constant is further de- 
creased. The extent of sensitisation or stabilisation increases with the increasing 
concentration of the added nonelectrolyte. The effect of dioxan is rather more 
pronounced, when the time of coagulation is small than that observed when the 
time of coagulation is large as brought about by adding smaller amounts of coa- 
gulating electrolyte. Again the degree of sensitisation or stabilisation is more re- 
markable with the sols having greater amounts of peptising electrolyte i,e, with one 
having lower pH values than one with higher pH. Thus the sols of greater 
purity shows comparatively lesser sensitivity towards the action of non-electrolyte, 
dioxan. 
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The very nature and the behaviour of the various sols prepared by different 
methods differ remarkably in various respects due to their varying internal 
structures as we have reported earlier. 7 It is evident here that the sols obtained 
by Kreck’s method behave differently than those prepared by peptising the hydrous 
oxide by an acid. 

It may be pictured that with the decrease in the dielectric constant of the 
medium, the double layer around the colloidal units becomes mobile and there may 
be self neutralisation of the charge of the particle leading to the unstability. The 
surfaces forces operating in aggregation process is also of importance in defining 
the stability of a sol, which is different for sol particles obtained by different 
methods. The structure of the colloidal units and the surface of the different sol 
samples is obviously not the same as the coagula of hydrous ferric oxide sol preci- 
pitated from Kreck’s sol is more compact than those obtained from other sol 
samples, where the coagula are highly flocky and have a tendency to form a gel 
like structure. 

Further work is in progress in these laboratories in this direction and will be 
reported in subsequent publications. 
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ABSTRACT 

The need for irrigation iri a state like Uttar Pradesh and the importance for making the optimum 
use of irrigation waters therein has been explained* The latter has been advocated by the correct 
determination of water requirements of crops. An account has been given regarding the experimental 
studies conducted for three principal crops of Uttar Pradesh, As water requirements vary from one 
region to another and from crop to crop, the initiation of intensive research by field studies, which can 
be checked up with formulae given by the United States Department of Agriculture, is recommended 
not only for a better utilization of irrigation waters but also for the maximization of food production 
which is the burning problem of the country. 1 

INTRODUCTION 

Irrigation waters are used for fulfilling supplemental needs for crop production 
not met with by natural precipitation. While underirrigation brings down the 
crop yield, it may not be so well known that over-irrigation also affects agricultural 
production adversely. Wherever water is available for irrigation to the 
cultivator, he invariably over-irrigates his fields under the mistaken notion that the 
higher the application of water for irrigation the greater the yield. This practice 
is causing a double loss by reducing, the yield (besides affecting the quality of the 
crop) and also by the waste of valuable water for irrigation which could have been 
utilized for irrigating other fields and thus raising the yield. 

In view of the chronic food deficit and limited water resources for irrigation in 
the country, it is imperative that every drop of water should be put to best use for 
crop production. Prerequisites for the optimum use of irrigation waters will be a 
knowledge of the correct water requirements of crops grown in any area and 
dissemination of the same among the cultivators. With this object in view, a series 
of experiments were laid out for the determination of water requirements of three 
principal crops of Uttar Pradesh, namely, wheat, rice and suaar cane at a number 
of experimental stations. 

It may be mentioned that the; maximization of agricultural production 
depends upon proper land and water management. The present series of experi- 
ments undertaken in Uttar Pradesh were related to proper water management 
alone, while regarding the aspect of land management local cultural practices were 
adhered to as far as possible. For proper management of water which was the real 
factor experimented upon, the object of the experiments was to work out the best 
combination of the level; of irrigation;(or the depth of irrigation treatment) and the 
frequency of irrigation (interval between irrigation treatments) for obtaining the 
maximum yield. : . .. 

The aspect of singling out the factor of proper water management had not 
been touched upon in, previous experiments on the subject in the country. In 
Uttar Pradesh, some experiments were carried out previously at the Sugarcane 
Jte$earch Statlbn, ShahjaH'anpur and Rice Research Station, Nagin^ for the 
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determination of water requirements of sugarcane and rice respectively, but the 
watering needs were studied at these stations along with varietal trials and tnanurial 
responses. Similar experiments had also been undertaken at Padegaoa (Poona) for 
wheat, Gattack and El/ lcrabad for rice. But, as mentioned earlier, no attempts 
had been made therein to single out the aspect of watering alone. 


FUNCTION OF IRRIGATION WATER 

The ideal in irrigation is to keep a uniform soil moisture condition during the 
growing period. However, a water supply that will support this ideal is almost 
never attained because the period of ample supply is seldom coinc dent with that 
of greater demand by plants. 

Where water is applied for irrigation r the soil in the proximity of the water is 
immediately saturated, and the gravity portion* of the water seeks lower levels under 
the influence of both gravity and capillary attraction. The wetted soil reiains a 
sufficient quantity of moisture to satisfy its capillary capacity, and all the surplus 
continues to sink, filling the soil to the capillary capacity as it goes until ground 
water is reached, when the surplus is consumed in augmenting the ground water 
table or drains away in whatever channels are available. 

All water that is in excess of the optimum water supply and that escapes to tfie 
ground water table or passes off as drainage is usually wasted, for it carries mith it 
in .solution plant food. This also results in the rise of the ground water table with 
its accompanying problems of waterlogging. 

At the end of each irrigation, the upper zone of the soil is saturated to a 
cpnsinerable depth, and the gravity Water in this zone slowly descends to lower 
levels. The moisture of the surface soil soon gets reduced and an upward movement 
begins of capillary water to make up for the losses from evaporation from the 
surface. Thus, after each irrigation, a condition is soon reached, in which a narrow 
zone a . short distance below the surface is in astate of capillary saturation, and the 
percentage of moisture diminishes both upwards and downwards, at a nearly 
uniform rate. If at this time the moisture at the lower limit of the root zone is just 
a little above the wilting coefficient, the water has been economically applied, and 
none wasted. It is therefore imperative that irrigation should be applied before the 
soil moisture reaches the wilting point throughout the root zone, otherwise, the 
plants will be severely injured or killed. 

It is thus clear that under-irrigation does not fulfil the water requirements of 
crops and thereby damages the production of the. latter. On the other hand, over- 
irrigation means the application of water in excess of the watering needs. The 
latter is thus merely a waste of water for irrigation, apart from the fact that the 
yield and quality of the crop are damaged. This waste of water has other harmful 
effects as well. The water table keeps on rising steadily owing to the percolation of 
the excessive quantity of water, which in the course of some years may result in 
waterlogging with its accompanying disasters in the form of the creation of swamps, 
damage to public health and sanitation by the spread of malaria, salt encrustation 
with resulting loss in soil fertility. 

It may therefore be repeated that in any. agricultural planning, it is essential 
that both under and over-irrigation should be averted, which is possible only by 
the .optimum use of irrigation waters by the determination of correct Water 
requirements of crops. 
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experimental APPROACH 

With the objectives set forth earlier, experiments on water requirements of 
wheat, rice and sugarcane were laid out since 1942-43. The factors studied in these 
experiments were the level and frequency of irrigation, while other factors such as 
type of seed, type and quantity of manure and other cultural practices were made 
to conform to local practices as far as practicable. The experiments were taken up 
at a number of stations which were typical of varying climatological and 
topographical features of the State of Uttar Pradesh, namely, Bara Banki, Bulandshahr, 
Nagina, Tisuhi (Mirzapur), Atarra (Banda), Bahadrabad, Dhanauri. They were 
replicated over a number of years in order to do away with the effects of climatic 
variations. 

The experiments were conducted with parallel yield trials of different irriga- 
tion treatment combinations of level and frequency of irrigation, based on 
randomized block system. In this method, the experimental area was divided into 
a number of sub-plots on each of which different dosages and frequencies of irriga- 
tion were applied in accordance with a set design. Randomizations in the designed 
experiments. were altered from station to station and from year to year for each 
crop. The size of sub-plots was governed by the area available for experimenta- 
tion at each station. Experimental area ranged between 2 and 4 acres at different 
stations. Subplot sizes were 80' by 16-5', 40' by 22' or 75'by 15', 25' by 11', 25'by 11', 
25' by 11 ',73' by 15',73'by 17.33' for Barabanki, Bulandshahr, Tisuni, Atarra, Nagina 
Bahadrabad and Dhanauri respectively. Type of soil was loam at Barabanki and 
Bulandshahr, black cotton or hard clay soil at Tisuhi, parwa at Atarra, clayey at 
Nagina, loam mixed with sand and clay at Bahadrabad and Dhanauri. 

Particular care was taken to .achieve the maximum accuracy in the measure- 
ment of depths of irrigation. The discharge at the head of the irrigating gul was 
usually measured by a 90° V-notch of standard specifications, as detailed in the 
pamphlet B. S. S. 722-1937. The discharge over a 90° V-notch is calculated by the 
following equation: — 

Q, = 2-49 Ha-*", 

where Q, = discharge in cubic feet per second, 

H observed head in feet. 


EXPERIMENTAL RESULTS 

Yield data for each station was analyzed for each crop every year. A pooled 
analysis was also carried out at the end of each series of experiments, and the best 
comoiration of irrigation treatments for obtaining the maximum yield was thereby 
derived for each crop for the different experimental station: 
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(A) Wheat 

tor one irrigation treatment the best yield was obtained with three inches for 
BaraBanki, Bulandshahr, while this depth was found to v be four inches for Atarra. 



‘Depth dHAijga^bhihilhefei;’* :: - - 

Graph No. 1. Shovving . Yield of Wheat at Atarra fo-r -one Irrigation ~ 

; : ,_\-V ,1 Treatment- for 19.17*48 . • ' 


From Graph No. 1 for the yield of wheat of Atarra for r 1947*48, it will be clear that 
the yield goes on increasing with increasing depths li’jfioHt inches above which the 
yield again goes down, which brings out clearly the harmful effects of both under 
and over-irrigation. ; ...... 

With two irrigation treatments, the optimum depths were , found to be the 
same with an interval of three to four WehkS'i Table -No. l gives an idea of the 
yield data for Dhanauri for the year 1957-58. 

With the limited range of experiments, it has not been possible to study 
the water requirements for various types, of soils prevailing in Uttar Pradesh, 
although it has generally been found that the requirements were higher for sandy 
sous as compared to clayey soils, indicating that the soil-moisture retentive capacity 
does play a role in the watering needS: ; ~This is in conformity with the . findings of 
similar experiments carried out in detail by the United States Department of 
Agriculture. It Was also foundjthat the climate or the humidity did not have any 
pronounced effect on the water requirements of a particular crop, although the first 
two were the governing factors in the selection of the cropping pattern in any area. 
It is, however, the rainfall which has a significant effect on water requirements. 
This has been elucidated below to some extent. 
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' *. it was brought out that in the drier tracts : of Western Uttar fradeshVrfepreSented 
by Bulandshahr with an annual rainfall of 26 inches, it was essential to have two 
irrigation treatments, whereas it was not of such primary importance in a wetter 
place such as Bara Banki with a rainfall of 36 inches* The average percentage 
ration of the maximum yield due to two and one irrigation treatments was 111 for 
Bara Banki, whereas this percentage was 147 for Bulandshahr. 

A higher optimum depth of irrigation requirement for wheat in Southern Uttar 
Pradesh (four inches against three inches for Eastern and Western Uttar Pradesh) 
could only be explained by the important effect of the distribution of rainfall* winter 
rainfall being practically negligible in Southern Uttar Pradesh unlike that in 
Eastern and Western Uttar Pradesh. It was also evident that it was not the total 
quantity of rainfall which mattered most but rather the distribution or pattern of 
rainfall, for the water requirements for Bara Banki with an annual rainfall of 36 
inches were not less for wheat than those for Bulandshahr with an annual rainfall of 
26 inches, This could only be explained by the fact that winter rainfall during 
the growing season of wheat was equally favourable at both the places. It may be 
mentioned here that in some places of Western Australia it has been possible to 
grow wheat with an annual rainfall of 10 inches only -because of a favourable 
distribution of rainfall during the growing season. 

The effect of water table conditions was not detectable on water requirements 
of wheat. This factor has not been mentioned in any of the studies conducted by the 
United States Department of Agriculture, 


(B) RICE 

For this important crop, for which experiments were initiated .only' for Tisuhi 
and Atarra and could not be continued over a number of years; a pooled analysis 
for Atarra for the years 1947 to 1950 indicated the optimum depth of" irrigation to 
be nine inches, no conclusive results being obtained regarding the frequency or 
interval between irrigation treatments. ~ - 

(G) SUGAJk CANE ' " ' " * ~ : 

The experiments were carried out at Bara Banki arid Bulandshahr, ' Four 
inches^ depth at three weeks .gave the best yield for Bulandshahr, while the 
maximum yield was obtained, with., five inches^ depth at intervals of two to three 
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weeks for Bara Bank*. Graph No. 2 for the yield of sugar cane at Bulandshahr for 
1948-49 shows clearly how for the best frequency of 3 weeks, the yield is considerably 
reduced not only by reducing the depth of irrigation but also by increasing the 
same from 4 inches up to 6 inches. Appendix 1 gives an idea of the analysis of 
variance of the pooled data for Bulandshahr for the years 1946-47 to 1947-48. 



Periodicity in Weeks 

Graph No. 2, Showing Yield of Sugar Cane Go 421 for 1948-49 at 
Bulandshahr 

The results were interesting in that the irrigation requirements for Bara Banki 
which is actually wetter with 43 inches of rainfall as compared to 28 inches for 
Bulandshahr were found to be higher, indicating again that the total quantity of 
rainfall was not the most dominant factor, the pattern of rainfall playing a more 
important role. 

empirical approach 

Apart from the experimental approach indicated above, aft empirical 
approach for an assessment of the water requirements from climatological date, as 
advocated by H. F. Blaney of the United States Department of Agriculture, was also 
attempted for wheat for Dhanauri for the years 1955-56 to 1957-58* vide Table No. 2, 
which is self-explanatory. The months of November to February, representing 
the growing season of the (wheat) crop, were taken into consideration* The value 
of K, the constant in the equation U — KF, was assumed from the experimental 
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data of the United States Department of Agriculture as 070 which is recommended 
for grains with a growing season of 4 to 5 months, while field irrigation efficiency, 
E in the formula I = (U-R)/E, was assumed to be 7Q}' # Irrigation requirement at 
the head of the field was thus found to be 1049 inches. From actual experimental 
results, it was found that two irrigation treatments of 3 inches each gave the 
maximum yield. If to this be added a paleo (pre-sowing) irrigation of 3 inches, the 
irrigation requirements work out to 9 inches which is not much different froth that 
obtained by Blaney’s approach, the more so if the seepage losses are also taken into 
consideration. 

CONCLUSIONS 

The need for the correct dosage of irrigation for the maximization of crop 
production has been emphasized. The experimental approach is recommended for 
typical stations representative of various climatological ' and topographical 
classifications. Where this is not possible, the empirical approach of H. F. Blaney 
may be utilized for an assessment of water requirements from climatological data. 
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' \ . t i APPENDIX No, 1 

Pooled Analysis of Yield Data of Sugar Cane Experiments conducted at Agricultural 
School Farm,' Bulandshahr during the years 1946-41 and 1947-48 with 3, 4 # 5, and 6 inches 
of irrigation at intervals of 2, 3, 4, 5, and 6 weeks. 

/ f The. final analysis of variance of Yield of sugar cane of pooled data in (Seers)* 
is given m Table I, m which items significant at 5% and 0*1 % are marked with one 
and three asterisks respectively, 40 seers being equivalent to 82 lb. The significance 
of intervals has been tested by years X intervals; that of depths with years, X depths, 
and the interaction of intervals and depths by years x depths X intervals, 

TABLE I 

Pooled Analysis of Variance of yield of Sugar Cane in (Seers) 11 


Variation due to 

Degrees of 
freedom 

Sum of 
squares 

Mean. 

square. 

Years 

1 

2,979,270 

2,979,270 

Blocks 

4 

25,724,1666 

6,431.0417 

Intervals (Weeks) 

4 

24,994,2167 

6,248-5542 

Depths 

3 

65,109-0000 

21,703*0000# 

; Week X Depths 

12 • 

92,750-9166 

7,729-2430*** 

. Years x Weeks. 

4 

5,447-8833 

1,361-9708 

- ,, Years x Depths. 

3 

6,807-1334 

2,269-0445 

Years x Weeks X Depths 12 

2,603-7834 

216-9820 

Error 

76 

42,281-1677 

556-3311 

Total 

119 

3,244,989. 



The following inferences may be drawn from the above pooled analysis ;~ 
Weeks : Testing this with years X weeks, we find it to be insignificant. ’ 

level 5pmbabliity. tia Now S w forSThe crfe^ff 6 h j S si S nificant at 5% 

sub-plots. “ the cntlcal difference for depths based on 30 


Critical difference = V 2269-0445 x 2 


30 


X *5’/ 


(3) 
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TABLE I. 1 


Summary of Results of Mean Yield in descending order with regard to Depths in the pooled Data. 


4" 

5" 

3" 

6" c. d. at 5% level 
of probability 

Mean yield in seers 
per sub-plot 401 2 

370-9 

359-8 

336*4 39*11 seers 


Comparing this critical difference (39-11) against the mean yields due to effects 
of different depths, we find no signficant differences between the mean yields due to 
four and five inches of irrigation. 

Depths xWeeks : — The variation due to the interaction of depths and weeks, 
when compared with the same due to years X weeks X depths, is highly significant. 
Forming the critical difference at 5% level of probability, we get 

Critical difference = v 2 16-9820 x 2 t r 
g X 5% 

= 18-52 seers per sub-plot 


TABLE I. 2 


Summary of results of Mean yields due to the Effects of Interaction of Depths x Weeks in the 

Pooled Data. 


Week^/ Depths 

3" 

4" 

5" 

6" 

c. d. at 5% 

2. Weeks 

356 

406-7 

428-3 

325-8 


3, Weeks 

356-3 

451-1 

360-5 

307 


4. Weeks 

359 

439-2 

365-5 

374 

18-25 seers 

5. Weeks 

356-2 

377-3 

367-2 

335-2 

per sub-plot 

6. Weeks 

371-8 

331-5 

333*3 

360 



Testing the critical difference (18 52) in the above table against the mean 
vidds due to the effect of this interaction, we see that there are no signficant 
differences between the mean yields due to 4 inches of irrigation at intervals of 
3 and 4 inches. Both are eqaully effective. 


CONCLUSIONS 

Four inches of irrigation, given at an interval of 3 weeks, may be considered 
to be the best combination in producing the maximum yidd of sugar cane followed 
by the same depth at an interval of four weeks for optimum yield. 
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Showing yield of wheat crop (Pb-591) at Government Agriculture Farm, Dhanauri for the year 1957-5$ . ■ ■; 

1 ’ " ^ Actual dates ofirrigation 
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Afotes : — 1. Dates of sowing 6-11-57 to 9-11-57. 

2. Dates of harvesting : — 15-4-58 to 20-4-58. 

3. Rainfall during the crop season: — 5*45 inches* 
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F = sum of monthly consumptive use factors, (/) for growing or irrigation season. 
I = (U — R)/E = irrigation requirement at head of the field. 

Here F = 17.06 inches, so U = KF = 0.70 X 17.06 = 11.94 inches. 

I = (U - R)/E = (11.94 - 4.60)/0.70 = 10.49 inches. 
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ABSTRACT 

The Brinkly and Kirkwood theory of propagation of shoek wave has been extended to the case of an 
inhomogenous medium. The mechanical decay of shock wave of Mach number five has been studied for 
three different ratios (0*1, 1 and 2) of the scale behind the shock front to the scale height of the 
atmosphere. The shock wave is found to decay in the first two cases. In the last case, after decaying 
over a short distance it begins to increase in strength. 


INTRODUCTION 

Brinkley and Kirkwood (1947), gave a theory of propagation of a shock wave. 
The chief feature of this theory is the decaying character of the shock wave as it 
propagates in a homogeneous medium. This theory has now been extended to take 
into account the variations of density and pressure in the undisturbed region. A 
similar extension has also been done by Ono et al (1962) and Kogure and Osaki 
(1962). The decay of a shock wave in an inhomogeneous stellar atmosphere in 
some special cases has been studied here by integrating the resulting differential 
equations. A shock wave of initial Mach number five has been assumed to 
propagate in an inhomogeneous atmosphere. However, the effect of gravitational 
field has not been taken into account. The variation of the scale height in the 
atmosphere has also been neglected. Three different cases specified by the ratios 
0T, 1*0 and 2‘0 between the scale behind the shock front and the scale height of 
the atmosphere have been considered. It is shown that in the first two cases the 
shock wave decays as it proceeds in the atmosphere, but in the last case it starts 
gaining strength after decaying over a short distance. It is observed that for a 
shock to decay, according to this theory, in a stellar atmosphere, it must propagate 
in a region where pressure and density graiients are relatively small, or in other 
words the scale height of the atmosphere is larger or of the same order as the scale 
behind the shock front, 

BASIC EQUATIONS 

Consider a plane, one dimensional shock wave propagating in an inhomo- 
geneous medium. Neglecting the effect of gravitational field, the equations of 
hydrodynamics, at the shock front, can be written in a convenient form as : 


_P 

Po 


8 u ■ 1 dp „ __ 

0* ^ >6* ' dt 


o, 


( 1 ) 


3L + _L. 4 - J— . = £>, 

3 1 t’o 3 * Po dx 


( 2 ) 


where u is the particle velocity, p the pressure excess over the pressure P 0 of the' 
undisturbed medium at rest. P the -density at any instante and p 0 , the density of 
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the uhdisti'ubed medium, x denotes the Lagrangian coordinate of the shock fron 
and C the velocity of sound. The undisturbed pressure and density are functions 
of space coordinate x only. 

The physical quantities on the two sides of the shock front are related 
by the Rankine - Hugoniot relations : 


P = p'o • u.U, 

(3) 

P{U-u) = P 0 U, 

(4) 

AH = l( — L + J_ \ 

2 ' P 0 ’ p / ’ 

(5) 

where U is the velocity of the shock front and AH is 
experienced by the gas in crossing the front. 

the specific enthalpy increment 

Applying the operator, 


d _ 3 , 1 a 

dx dx' U ‘ dt 

(6) 


to eqn. (3) one gets 


| jj on k dp ^ k dp 

& dx P 0 U ‘ dt ~ Q ’ 


where 


= '—u —- . dp 0 u dU dp Q . uU dp 0 

d P o * dx dP Q dx P Q " dx ' 




The fourth equation in the Brinkley-Kirkwood (1947) theory is 


dt + u 


where 


D (*) =Jp’.u'.dt (10) 

*<>(*) ■ 

CO 

= jp 0 .h.d | 0 ( 11 ) 

X 

and A, is the specific enthalpy increment of an element of the fluid on crossing the 
shock front and thereafter returning to pressure P 0 . g is the displacement at time t 
of the gas element whose undisturbed position is x . v, is a very slowly varying 
function. The value of v vanes between 2/3 and J. In the case of a homogeneouf 
medium one can assume that v is independent of*, which is equivalent to imposing 

[ 170 ] 



Sknilatity restraint, This assumption may not be exactly true in tlae present 
However, in view of the limited range of its variation one can simplifity the problem 
and use a constant value. In the present work a value 2/3 is used. It corresponds 
to the asymptotic behaviour of the decay curve. p> and u‘ are the excess pressure 
and the particle velocity respectively behind the shock front. P 


From eqs. (1), (2), (7) and (9) one can eliminate du/dx and 0a/9i and determine 
dpldx and 3 p/fit Substituting these in the equation 


d P = dp , 1 dp 

dx 3 * ‘ U ' dt 

one gets, 

dp v G . P 3 

~Tx ' 15(4" ' 2(4+ 1) -G ' ~pjw~ 


1 

2(ft+l)-G 


[{( 2 - 0 - 



3_£ 1 dJP„ 
3 P a i dx 



The eqn (10) gives 


d D 
dx 



( 12 ) 


d IL+ 2 ± ’> jgnl 

dP 0 dP 0 J dx J 
, ^ . (13) 

(14) 


isothermal jump across the shock front 

Sd far eqs, (13) and (14) are exact. For application to actual cases, 8ne has 
to make some assumption regarding the path of the gas element behind the shock 
front. As the temperature across the shock front is extremely high, the radiation 
flux would be very large and fast. As in earlier works (Odgers and Kushwaha, 
1959 ; Bhatnagar and Kushwaha, 1961 a and b ) it is assumed that due to this 
radiation flux the temperature is equalised instantaneously on the two sides of the 
shock front. After this the gas element returns adiabatically to the undisturbed 
pressure P 0 . Thus the pressure and density across the shock front may be related 
as follows : 


— 1 4“ s » 


*o P 0 

where s , the shock strength is defined by the relation 

P = .... . . : 


With the help of the Rankine - Hugoniot equations one gets , 

' u — j, , Co ‘ A 

A/—*' <i+U* ’ - 

o = -A • (H-U* * 

V Y 


(15) 


(16) 

(1.7) 


( 18 ) 
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jU i _ £ __ 

2(*+I) 

G = 1 L_ 

and 

Ca = C o a = y . - A 

Po 

Introducing these in eqn (13), one gets 


0 $) 

( 20 ) 


( 21 ) 


p„ . * 3 , y* a +y* dlnP 0 

D(x) 2ys+'iy+l *+l ^ 27 ^+ 37+1 hT~ 

27^+37J+y+t dlnP Q 

27M-37+1 ‘ dx ' ( 22 ) 


h, can be written as a sum of the specific enthalpy increments A H, on crossioe 
the shock front and A H that on returning adiabatically to the pressure P a . Then 


h = 






y - 1 V p* 0 


p_ 

p 


)■ 


where P 0 * is the final density, 
tbe shock front, therefore 


As the gas element follows an adiabate after crossing 


P 
P o 



and hence 


y = i+. . 


(23) 


dD 

dx 


= - p 0 k 




(24) 


INITIAL VALUES 


„wv^° Study bow th ? sbock stren g th s > defined by the relation (16), changes as the 
for lnv r n n P r a Hr at i eS ’ ° De haS , t0 5ntr g rate the two equations (22) and (24) numerically 
iv! C ! n aSC UDder con9lderadon ‘ Suppose that a shock wave of Mach 

efluifibrium^stat ^ St3rtS pro P a ° atin " in an inhomogeneous atmosphere. In 
equilibrium state, the pressure at any point x is given by 


P 0 = 


Pi 


s*m 


( 25 ) 
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where ft is the value of at * - 0 and H 18 the scale height of the atmosphere. 
Strictly speaking H should vary in space, but for the sake of simplicity as an 
approximation this variation is neglected. From the definition of the scale height 
(Lamb, 1939) one can show that if the variation of H is not considered then 


dlnP 0 _ dlnpg 
dx dx 


(26) 


Introducing the dimensionless quantities, £ and y defined as follows 

D(x) = P 0 • H . £ , (27) 

and 

x — H .y , (28) 

the eqs. (22) and (24) reduce to 

y y^+v-1 * a j. Yi a +2yi+s+y+l . /nm 

dy £ 2yr+3y+l' f+1 2yj+3r+l K ' 

and 

(30 > 

The second term in eq. (29) and the last term in the eqn. (30) are the 
contributions due to the variations of undisturbed pressure and density. The 
steeper the gradients of pressure and density, larger is £ and hence slower is the 
decay. We can take 7= 5/3 for a monatomic perfect gas and v = 2/3. With the 
help of eqn. (17), initial value of s = 40 is taken for a shock wave of Mach number 
five. The initial value of £ is found below : 


1 oo i co i j 

ppr J p 1 .u' .dt = i rf w . dt = it * 

a*) *<>(*) 


where p 1 has been averaged as l . P 0 ; l being a constant. Following Odgers. and 

r° 

Kushwaha (1959) we arbitrarily take l = 10. The integral j u' . dt = d gives the 

<o(*) 


scale behind the shock front or the distance travelled by the material. It can be 
determined by integrating the observed radial velocity curve for any particular star. 
We have considered three values of the ratio of d (the scale behind the shock front 
to H (the scale height) namely, O’ I, 1 and 2. These give the initial values of £=1*0 
10-0 and 20’0 respectively. 


DISCUSSION OF THE RESULTS 

The eqs. (29) and ( 36 ) have been integrated numerically following Harm and 
Schwarzschild Scheme (1955) for the above mentioned three different ratios of the 
pcale behind the shock front to the scale height of the atmosphere. The values of 
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£5 s f u/C 0 and t)/G 0 obtained in this way as the shock propagates have been * given 
in tables 1, 2 and 3 for the three cases respectively. u/C 0 and U/C 0 have been plotted 

TABLE 1 

Decay of a Shock wave in an Inhomogeneous Atmosphere — (d/H = OT) 


y 


S 

u/Co 

U/Co 

x x 10 xo cm. 

0 

1-000 

40-000 

4-84 

4-96 

0 

, o-oi 

0-835 

35-324 

4-54 

4-67 

0-50 

0-02 

0-691 

30-995 

4-25 

4-38 

1-00 

,.0-03 

0-567 

27-007 

3-96 

4-10 

1-50 

0-04 

0-459 

23-357 

3-67 

3-82 

2-00 

0-05 

0-368 

20-039 

3-39 

3-55 

2-50 

0-06 

0-291 

17-048 

3-11 

3-29 

3-00 

0-07 

0-228 

14-373 

3-84 

3-04 

3-50 

0-08 

0-176 

12-007 

2-58 

2-80 

4-00 

0-09 

0-134 

9-936 

2-33 

2-56 

4-50 

0*10 

0-100 

8-147 

2-09 

2-34 

5-00 

; o-n 

0-174 

6-624 

1-86 

2-14 

5-50 

0-12 

0-054 

5-344 

1-64 

1-95 

6-00 

0*13 

0-039 

4-286 

1-44 

1-78 

6-50 

' 0-14 

0-028 

3-424 

1-26 

1-63 

7-00 

0*15 

0-020 

2-731 

1-10 

1-50 

7-50 

0-16 

0-014 

2-179 

0-95 

1-38 

8-00 
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tABLE 2 

t)ecay of a Shock wave in an Inhomogeneous Atmosphere — (djU = 1*6) 


y 

C 

5 

ujGo 

UjCo 

* x 10 10 cm. 

o-oo 

10*000 

40*000 

4*84 

4*96 

0*00 

0-20 

8*412 

34*360 

4*48 

4*61 

1*00 

0-40 

7*049 

29*440 

4*14 

4*28 

2*00 

0*60 

5*887 

25*165 

3*81 

3*96 

3*00 

0*80 . 

4*900 

21*467 

3*51 

3*67 

4*00 

1*00 

4*067 

18*282 

3*23 

3*40 

5*00 * 

1*20 

3*369 

15*554 

2*96 

3*15 

6*00 

1*40 

: 2*785 

13*225 

2*72 

2*92 

7*00 

1*60 ' 

2*300 

11*250 

2*49 

2*71 

. 8*00 

1*80 

1*900. 

9*580 

2*28 

2*52 

9*00 

2*00 

1*570 

8*175 

2*09 

2*35 

io*oo 

2*20 

1*299 

6*997 

1*92 

2*19 

11*00 

2*40 

1*077 

6*012 

1*76 

2*05 

12*00 

2*80 

0*744 

4*500 

1*49 

1*82 

14*00 

3*20 

0*518 

3*441 

1*27 

1*63 

16*00 

3*60 

0*359 

2*685 

1*08 

1*49 

18*00 

4*00 

0*246 

2*126 

0*93 

1*37 

20*00 
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{’ABLE 3 

Computation, of 1 the Shock Strength of a Shock wave as it frogagates in ah 
Inhomogeneous Atmoshphorcr— (rf/H » 2*0) 


y. ■ 

. C 

s 

m/C 

U/Co 

x X 10 10 cm. 


' 20-000 

40*000 

4*84 

4*96 


0*2 

20-368 

39*139 

4*79 

4*91 

0-50 

0-4 

20*900 

38*495 

4*75 

487 

1*00 

0*6 

21*609 

38*058 

4*72 

4*84 

1*50 

0-8 

22*632 

37*821 

4*70 

4*83 

2*00 

1*0 

23*632 

37*781 

4*70 

4*83 

2*50 

1*2 

24*993 

37*932 

4*71 

4*84 

3*00 

1*4 

26*629 

38*276 

4*73 

4 * 86 . 

3*50 

1*6 

28*579 

38*816 

4*77 

4*89 

4*00 

1*E 

30-892 

39*553 

4-81 

4*93 

4*50 

2*0 . 

33*624 

40*494 

4*87 

4*99 

5*00 

2*4 

40*646 

43*026 

5*03 

5*14 

6*00 

3*8 

50*397 

46*509 

5*23 

5*34 . 

7*00 


against distance travelled by the shock wave in the figures l, 2 and 3. It is seen 
that the shock decays rapidly in the first case, rather slowly in the second, but in 



Fig. 1- The theoretically computed velocities have been plotted against the distance for d/H 0*1. 
Open circles give U/Co and dark circles give u/ Co. 
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Rgi 2. uj Co and U/Co have been plotted against the distance travelled by the shock wave From where 
it first becomes observable. For rf/H =r 1. Open circles represent U/Co and dark circles re- 
present w/Co. 



Figi 3* m/Go and U/Co versus x curve for d/H ss 2. Open circles represent U/Co and dark circles re- 
present u/Co. 

tb« third case it decays only over a small distance and thereafter starts gaining 
strength. From this it is obvious that for higher values of these ratios the growth 
would be still faster. Odgers and Kushwaha (1959) and Bhatnagar and Kushwaha 
(1961 a, and b and 1962) have used the Brinkley-Kirkwood theory to compute theo- 
retically the variations of the observed velocities of a /8 -Cephei star BW Vul treating 
the atmosphere to be homogeneous. This roughly corresponds to the first case 
considered here. Usually the scale heights are determined from the gravitational 
acceleration. Whitaey (1956 a and b) has shown that the effective value of gravity 
is reduced due to propagation of a shock wave. Recently Irosbnikov (1962) has 
shown from analysis of the behavour of , on the basis of a model in which the 

shock wave progagates in a uniform gravitational field that the approximate value 
of the effective gravitational acceleration is of the order of 40 cm/sec a for RR Iyrae. 
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In our present work we have not taken into consideration the gravitational 
acceleration. If we do so perhaps the second case will also give sufficiently rapid 
rate of decay and assuming the scale height to be given by this case for BW Vul, 
i.tf. equal to the distance through which the shock propagates, the effective value of 
gravitational acceleration for this star will be of the order of 50 cm/sec a . However 
there appears to be some uncertainty about the scale behind the shock front. The 
value adopted here was given by Odgers (1955) which appears to be larger by a 
factor of about five, Abt (1959) has given the mean total expension for this star 
to be twice the scale height of the atmosphere (calculating the scale height from 
gravitational accelerations). If Whitney’s reduction factor is applied to Abt’s value, 
then scale height of the atmosphere may be taken to be of the order of the scale 
behind the shock front. If it could be shown that in the relevant region of the at- 
mosphere such a scale height (of the order of the distance travelled by the material) 
exists, that the shock model can be usefully employed to explain the observed radial 
velocity variations. 
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ABSTRACT 

* The rate of production ofNe ao from 0 16 (a,y)Ne s ° reaction through its 5*64 (3-) and 5*80 (i-) 
Mev excited states is investigated and found to be considerably low in comparison with other helium 
burning reaction rates in the temperature rage 1-2X1 0 8 °K and p=10 5 g/ c,c. On the other hand N 1 * 
(a, 7) F 18 (j5+v) O 18 reaction rate through 4*651 Mev excited level of F 18 as calculated and found to be 
faster at these temperatures. This depletion of N 14 will facilitate heavy element synthesis through neutron 
capture because N**is a rapid absorber of neutron through N 14 (n,p) CW reaction. 

At temperatures in the range of 3-6XlO s °K the rates of formation of Ne ao and Mg 2 ^ are 
comparatively faster. At still higher temperatures the photodisintegration of Ne 30 becomes important. 
The photo dissociation and other heavy-ion reactions leading to the destruction of Ne s 0 are considered. 
Rates of production of Ne aa and Mg* 6 through successive helium captures are considered at these 
temperatures. 

INTRODUCTION 

Helium burning reactions start in the core of the red giant stage of stellar 
evolution after the exhaustion of hydrogen as the energy producing material. The 
first reaction that takes place in the helium core at T~10 8 °K is the triple alpha 
collision forming G ia in its 7*66 Mev excited state. Then subsequent captures of 
helium by C 1 * produce O 18 and so on. Hayakawa et al (1956) calculated the 
reaction rates leading to the formation of C 12 , O 16 and Ne aa under the condition at 
temperature ~ 10 9 °K. The rates of forming C 12 and Ne 20 were found to be larger, 
Nakagawa et al (1956) made detailed calculations on the helium capturing reactions. 
Many other authors (Burbidge et al, 1957 ; Salpeter, 1957) studied the fusion of 
helium as a mode of element synthesis in stars. The helium burning process^ 
were found to be responsible for the synthesis of C ia , O 16 , Ne a0 and perhaps Mg* 4 - 

Recent experimental investigations and results thereof warrant a re-investiga- 
tion of the reaction rates during the helium burning phase of stellar evolution. 
The production of Ne^ 5 was expected to be due to the 4*97 Mev excited state 
of Ne*° whose spin and parity where supposed to be 2+. The formation of Ne a ° 
through O 16 Ne ao reaction requires that this level should have either odd^spin, 
odd-patity or even-spin, even-parity. But experiments (Litherland et al, 1961 ; 
Almqvist et al 1961) showed that this level has even-spin and odd-parity and hence 
might not be formed by a — capturing reaction in 0 16 . There are other two leavels 
ofNe a0 inthe relevent energy regions of Ne a ° having odd-spin and odd-parity- 
Thus it is expected that these two levels will contribute to the resonance formation 
of Ne a ° in stellar interiors. 

In the present paper attempt is made to examine the rates of formation of 
helium burning products with the newly acquired experimental informations. 
Helium bruning reactions at different temperature radges are considerd. Tempera- 
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tures In the range 7s* 5x10* °K to 6X 10* °K and higher may be expected to be 
prevalent in the late type giants or pre-supernova stages and it is believed that Ne®* 
will be produced at these stages. 

It was suggested (Cameron, 1957) that hydrogen from the envelope may be 
mixed to the expanding helium core in globular cluster star at the tip of the giant 
branch. Considerable amounts of carbon cycle products would then take part in 
the helium burning reactions. Neutrons, produced through C ia (a, n) O 1 * which 
begins even at slightly lower temperature (<^8 x 10 7 °K) will soon be depleted by 
N 14 (n, p) C 14 reaction and enough neutrons will not be available for heavy element 
synthesis. To facilitate heavy element synthesis it is required that N 14 should be 
-used up by N 14 (a, 7) F 18 (jBf+v) O 18 reaction (Burdidge et al, 1957), so that neutron 
consumption may not take place on a lavish scale. N 14 («,7) F 18 Q3 + v) O 18 resonant 
reaction, therefore, is considered in the present note and similarly O 18 (a, 7) and 
Ne®* (a ,7) non-resonant rates are also studied. 

The non-resonant rate of formation of Ne aa is quite low and so is the Ne** 
(<x,7) Mg 28 reaction. Considering that at higher temperatures (r g >5) several 
resonance levels are present within the ‘Gamow peak' in Mg a4 , the reaction cross- 
section for Ne®* (a s 7) Mg a4 is obtained as an average over many resonances and 
the reaction rates are found to be quite faster* 

The photo-dissociation rates at higher temperatures for Ne a ° are calculated. 
At higher temperatures when heavy-ion collisions can take place the destruction 
of Ne a ° through C la +Ne® > and O^+Ne 2 * 0 are also considered. 

Protons and alpha-particles that are produced in the heavy ion reactions are 
being utilized for still farther heavier element synthesis. 

RESONANCE AND NON-RESONANCE REACTIONS 

The rate of formation of C^, Ne ao , and O 18 (through F 18 ) are considered to be 
through the resonant reactions because at relevent energy regions, excited states 
of these compound nuclei are found to exist. 

The resonant rate of formation of the compound nucleus is given by (Reeves 
and Salpeter, 1959) 

p% t* r 

Log,/ =« 12*69 4- logjo • — —jt - - --- 1 log M (AT 8 ) - Er'Scc., -1 (1) 

where, the widths Ty and V and resonance energy E& arc all measured in 
Mev, and tempetature is in 10* units. A is the reduced mass in AMU of the two 
interacting particles and A a that of helium nucleus. T a and Ty are respectively 
the alpha-particle and radiation widths. The total width + Ty g)=(2J+1)/ 

(2S+1) (21+1), where J is the spin (in the unit of A) of the compound state formed 
.and S and I are the spins of the two interacting nuclei. 

Wherever the experimental values of <aT a ty/r are known their use has 
been made in this paper. Otherwise, considering that T a << Ty ^T, it istakcu 
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that the second term on the right hand side of-(l) is reduced to «r c The alpha- 
particle width r« is obtained by using Wigner and Isenbud’s (1947) disperson 
relation of nuclear reactions which gives 

r«=2K a P,y* a , . ( 2 ) 

where K a is the wave number of the incident particle, Pf is the barrier penetra- 
tion probability for particle of angular momentum l and y ^ a is the reduced width 
of the particle in the level under consideration. The particle width is represented 
with a subcript « throughout the following because alpha-particle reactions are 
considered here 

P/ ^ G , (for Low energies) (3) 

where G / is the irregular solution of the Schrodinger wave equation for a particle 
in a Coulomb field. 

At very low bombarding energies the S - wave penetrability is (Cameron, 
1959*) 

P Q s G 0 ~ a ^ 27T7?0 o a exp .( - 2iry), (4) 

where y =* 0*1574 Z Q Z 1 A* / E*, E is the bombarding energy in C — M system in 
Mev. 

==* x% Kj 2 Qc) and x — (fipy)* 

where K x (#) is the modified Bessel function (Burbidge et al, 1957) given by 

K 1 (*) 21 ex P ( - *+ l A*) . • — 1 (5) 

and p *»KR. R, the interaction radius is taken to be 

R se 1*45 x IQ" ^ A 0 » + A x ^ ^ cm, 

y ^ is the reduced width given by (Blatt and Weisskopft, 1952), . 

y* / D = 2 X 10“ 14 cm. (6) 

a 

By knowing T5 from the energy level scheme one can determine y^ and thus T a in 

a a 

(2) Will be completely known, y in (2) is always smaller than y which is 

called the single particle limit of the reduced width A/m R, m being the reduced 
mass of tb c interacting nuclei). 
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For the non-resonant thermonuclear reaction between particle of type 0 and 
1, the number of reactions per nucleus of type 1 is (Burbridge et al, 1957), 

P..X = 4-34 X105 — -2- S (AZ 0 Z,) e _C sec- 1 (7) 

Aq 

In the present case particles of type 0 are helium nuclei, p is? the density in g/cc ; 
x 0 is the concentration by weight of alpha-particles ; A is the reduced mass in 
AMU of the two interacting nuclei and S, the cross-section factor in Kev-barns 
is given by 

S = <r E exp ( 3 1-28 Z 0 Z X A* E~* ) Kev, barns, (8) 


where <r is the reaction cross-section in barns and E is the particle energy in the 
CM system in Kev. E. can be determined from (Burbidge et al, 1957) 

E = 1-220 ( Z 0 * Z x a AT 6 « )V« Kev. (9) 

Consequently for different kinetic temperatures v can be determined by the 
usual methods. ^ in (7) is given by 


r 42*48 ( Z * Z x * A 


( 10 ) 


where T Q is the temperature in the unit of lO^K. 

For M* 44 , the cross-section is determined as an average over many resonances 
with the help of expression (Cameron 1959b) 

< > * ~j£i~ < * - — > X 10» barns E 11 ) 

<e 

where D is the average level distance and the other terms have their usual meanings 
With certain assumptions this expression can be written as 1 


< a > 




P 0 <«>- < ?a/D > X 10W barns. 


( 12 ) 


Ey knowing < a > from (12), one can calculate with the help of (7) and (8) the rate 
of formation of Mg s *. <<o> is taken to be equal to 3. 

For the non-resonant O 19 («, y) Ne 2a and Ne* 4 (a, y) Mg 29 first the cross- 
sections for (a, n) reactions are obtained (Cameron 1959a). Since T« < < r„ £5 Tj 
< <r > a >(a, n) TylT a . We have taken this ratio to be of~0’01. 

Thus < <r > (a, y) i are obtained and with the help of (8) and (7), the reaction 
rates are calculated, s 
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PHOTO- disintegration reaction rates. 


The rate of photo-dissociation involving any one nuclear resonance is given by 
(Cameron, 1959a) 


\ t - 1-50X1 o 15 -22^ y exp 11*61 $Q c .fe> sec., - ! 
®oi T 0 


03) 


where * 0 , and o> 01 are the statistical weights of the two interacting nuclc 0 and 

1 and their product respectively and 


V - - r « r y / r 

Er is measured in Mev and widths are in ev. Temperature is in 10 9 °K unit. A" 1 
gives the mean life-time of photo-disintegration rate for a given energy level. " 


CALCULATIONS 

The helium bruning reaction rates are calculated at different temperatures. 
Ne* production rate is calculated from the two excited levels separately. The" 
quantity «> r ffl T y / T for the 5-64 Mev and 5-80 Mev levels is calculated by Gove et 
al (1961) and Kuehner et at (1961) respectively. These values have been used in the 
present calculations. The quantity for 5-80 Mev level is found to be little less than 
0*15 ev and in the present case a value of 0-10 ev is taken and for 5*64 mev level 
«* F« r y / T » 0-003 ev. For the N 14 (a, y) F t8 resonant reaction rate, the 
different parameters for the 4*651 Mev excited state are taken from the work of 
Cameron (1959c). (y a =0-01, Y amax = 2-3X10“° ev. cm.). It is, however, assumed 
that this level will be formed due to the reaction under consideration. Er => 0*23 
Mev. The particle width is calculated by the method outlined above. T is 
found to be 4*0 X lO' 11 ev. Taking P = 10 5 g/c.c. (pure helium) all the reaction 
rates at these temperatures are calculated. 

The 3a-»C 14 reaction rates are calculated with a slightly modified value of 
Fy due to Alburger, (1961) for the second 7'66 Mev excited state of C 13 . It is 
assumed that the second excited state decays to the ground state through cascade 
y— ray transition, through the emission of 3- 23 Mev gamma-radiation. Fy is 
taken to be 2’5 X 10 ~ 3 ev with an accuracy of ~ 50%. The equation that is used 
for this set of calculations is taken from Burbidge et al (1957). 

The other source of neon in stars may be due to 0 18 which is produced by a — 
capturing reaction in N 14 . The non-resonant reaction rates are calculated by the 
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i \ , , , , , j * ; . ? j < ^ ! 

methods outlined above. The 0 16 formation rate is taken from the work of Burbidge 
et al (1957) and shown in Fig. 1. 



Fig. 1. Helium burning mean life-times vrs temperature. ( t is in years and temperature is in 
10 8o K-unit. px a = I0 5 g/c.c.) ( 1) and (2) against 0 I 6~>Ne 2 ° give the mean life, times 
of this reaction for the 5*64 and 5*80 Mev levels respectively. 


The reaction fates are shown in Fig* 1 for the temperature range 1 — 2 
X 10 8 °K* The formation rates of Ne 22 a Mg a6 , Ne a ° (from the two levels under 
consideration, separately) and Mg a4 in the temperature range 3-6 X 1 0 a °K are 
shown in fig* 2. 
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Fig. 2. Helium burning reaction rates as a function of temperature. (1) and (2) against 0**-* 
show the rates of 5*64 and 5 ‘80 Mev levels separately. 

Table I gives some of the mean life-times of reactions at a few different 
temperatures. 


TABLE I 

Logarithm of the mean life-times (£) in helium burning (yrs). 


r 9 

3He 4 -»'C ia C' a -»0 16 * 

0 16 -*N e 20 
(1) 

0 16 ~»Ne a ° 

(.2) 

(/3» +v )0 19 

0 ia ->Ne M Ne^Mg* 6 

1*0 

7*49 

7*57 

29-05 

36-09 

3-92 

19-06 

26-20 

1*2 

4-59 

5-94 

21-78 

27-39 

2*11 

- 


1*4 

2*56 

- 

16-60 

21-19 

. 1-09 

13-86 

20-10 

1-5 

1-76 

4-09 

14-53 

18-71 

0-32 

- 

- 

1-8 

-0-08 

- 

9-72 

12-95 

-0-85 

10-34 

15-99 

2-0 

-0-98 

1-91 

7-33 

10-08 

-1-75 

8-95- 

13-38- 


*This result is taken from Burbidge et al (1957), 
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At higher temperatures when V-rays acquire sufficient energy to disrupt 
the ground state of nuclei, photo-dissociation becomes a very important pheno- 
menon, Photo-disintegration rate depends on the disintegration energy which is 
comparatively smaller (4*75 Mev) in case of Ne a °. Destruction rate ofNe a °is 
quite considerable at temperature ^10 9O K which is shown in Fig, 3. 


' Fig. 3. 



The photo-disintegration reaction rates vrs temperature. The heavy-ion collision times 
are also shown against temperature. 


j°Ai l i 5 f e Tv?«rt t f UCU0n t rates tl}rou &h heavy-ion collisions such as 

G^+Ne^ and 0 16 +Ne ao have also been considered, Ne 5 * 9 collision with itself takes 
place at still higher temperatures and is therefore not considered here. 

RESULTS AND DISCUSSION 

The helium burging time-scale is supposed to be ^lO 7 years (Burbridge et al, 
1957). From the present calculations it seems that during this period negligible 
amount of O 16 is converted into Ne a ° due to helium burning. After a time ^ 10 7 
years has elapsed the reaction 0 16 («,7) Ne ao begirs at T X 10 9 °K. Therefore, an 
equilibrium situation with O 10 («,7) Ne aa and Ne a ° (<*,?) Mg a4 as constituents cannot 
be envisaged. Since nothing can be said about the mode of evolution of stars 
after this stage it is hardly possible to estimate anything about the helium burning 
reactions leading to the production of Ne^ 3 , Mg 24 , etc* which occur at temperatures 
higher than 2 X 10 8 °K. 
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At higher temperatures, supposed to be prevalent in the core of the late 
giants or pre-supernova stage of a star, the rate of production of Ne* 3 becomes 
quite faster. But, since life history of the star at this stage is not known the ratio 
of abundances is difficult to estimate with any certainity. However, the reaction 
rate of Ne*° (a,y) Mg* 4 is definitely an over-estimate at lower temperatures. But 
as helium is exhausted at lower temperatures, further helium burning reactions 
do not proceed unless a fresh supply of helium is available. At this stage, heavy 
ion reactions such as C 1 * (C ia ,a) Ne*° can provide a helium source ; but it takes 
place at temperature of ^ 6 X 10 8 °K (Cameron 1959b). Otherwise, the helium 
burning reactions (fig, 2) may be envisaged to be taking place at the supernova 
helium burning shell sources where the temperature is considerably higher. 

From table I, one can see that the reaction rates at T 6 =120, for all the 
reactions 3 a C 1 *, C 1 * (a, 7), are roughly the same, Ne* a production rate would 
however, be too slow by a factor 10 16 . s 

Ne**is formed through Q 18 (a, y) and it depends on N u (a,y) F 18 (j8 + v) 0 18 
reaction. The binding energy of an alpha-particle in F 18 is 4.42 Mev and there 
are excited states of this nuclei at 4M0, 4*651 and 4*74 Mev- If the 4*651 Mev 
level has the right spin and parity (only the level should not have 0+) then a 
resonance rate of formation of F 18 is faster. This F ia emits neutrino (half-life 
211*87 hr) with the release of energy ~0*3 Mev and goes over to Q ia (Gamei^n 
1959c). This rate is found to be faster than 3a C 1 * reaction at the helium 
burning temperatures (Table I). Therefore it can be expected that this reaction 
would compete with 3 a reaction as a source of energy generation at these 
temperatures and density as was also suggested by Cameron (1959c). It may be 
hoped that Ne** is made also, as was pointed by Cameron (1957), through Ne 19 
(a,y) Mg** (/3 + v) Na** (j3 + v) Ne.** Ne**, thus formed, may be destroyed by (j> 3 y) % 
(a } n), (a,y) and (n,y) reactions, Table II shows a cross-section taken from the work 
of Cameron (1959b). The last two values are calculated according to the usual 
method. It seems that Ne** synthesis will be considerably low in these steller 
conditions. 


TABLE II* 


Reactions with Ne** and their reaction cross- 3 ections calculated at T a = 6, 


Reactions 

< cr > 

N Na a3 

0*22 micro- barns*}* 

Ne aa (a 3 n) Mg as 

0*2 Micro-micro barns 

Ne aa (a, y) Mg a « 

0*002 Micro-micro barns 


# since Ty < < F Q ^ T, <1 a > (<*,7)/ < cr > («,«) ^ Ty / r n and this ratio 
is taken to be equal to 0*01. 

^Cameron (1959b), 
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DECAY OF INTENSE SHOCK WAVES IN STELLAR ATMOSPHERE 


By 

M. 3, BHATNAGAR and R. S. KUSHWAHA* 
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[ Received on 14th May, 1962 ] 


ABSTRACT 


The mechanical decay of a shock wave taking into consideration the interaction of material gas 
with radiation and including the loss of energy by radiation, has been studied with a view to interpret the 
radial velocity variations in jd -cephei stars. The temperature jump is found to be much less than that 
for a perfect gas. 


INTRODUCTION 

Bhatnagar and Kushwaha, in an earlier communication (1961 a), have 
studied the Brink-ley-Kirkwood theory (1947) of propagation of shock waves for 
„ n n ii ca tioTi to interpret the discontinuous radial velocity 'curve of 0-cephei stars, 
taking into account the interaction of material gas with radiation. Radiation 
nressure and energy density become important greatly before a shock may be con- 

Fidered relativistic For shocks propagating into normal air the shock velocity is 

u 6 tenth thl speed of light at a Mach number of about 10 5 , while radiation 
° ecmal to material gas pressure behind the shock at a Match number of 

Et 2 X \ll Also, radiative transport of energy may be important even though 
T.tinn nressure and energy density are negligible. For equilibrium between 
radiation and air at atmospheric density, radiation pressure is equal to material gas 
rl\nre at a temperature of about three million degrees absolute whereas radiative 
f Et ofenerv P becomes important at much lower temperatures. In order of 
SSnFe (with increasing Mach number the effects. are : radiative transfer, 
radiation energy density and pressure and finally relativistic nature of the shock,; 
if another wofk, we (1962) studied the decay of shock waves -including the loss of 
1 hEdiation for a perfect gas only. The ratio between the gas pressure and 
FbF ^(Ration pressure is of the order of 4 : 1 in the relevant parts of the atmosphere 
S a B 2 star (Bhatnagar and Kushwaha, 1961 a). Therefore in the present note, 
° f a ^i ® of the shock wave is being studied by taking into account all these 
SdiaS eff^ U radiation Soxf radia.ioa eacrgy dcasi.y aad rad.at.oa 

pressure. 

BASIC EQUATIONS 


We have at the shock front : 
(a) Equation of notion 


du i ^ 

W +X 


= o, 

3* 


(0 


* Present address ; Department of Mathematics, University ofjodh pur, Jodhpur. 
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(b) Equation of continuity 


du , __1 

9* ^ Pc * 


W + 
8 1 


au 

x 


* 0, 


( 2 ) 


where « «= ^ -M ^ is the particle velocity, p(— p t + />>•) is the total pressure 

excees over the pressure P 0 (=* P« 0 + Pr 0 ) of the undisturbed fluid. The subscript 
g is used for the material gas and r for the radiation, p is the density at an instant 
and P 0 that of the undisturbed fluid. t is the time, and i- the Eulerian coordinate 
at time t of an element of the fluid with Lagrangian coordinate x. c is the velocity 
of sound. The variation of P D with respect to x has been nenglected in comparison 
with that of p. The coefficient a takes the values 0, 1 and 2 for a plane, cylindrical 
and spherical shock wave respectively. 


The Rankine— Hugoniot relations can be written as 


P ■■ 

P{V - u ) 

AH 


P o . «/ u, 
P a U, 

_ 

2 \ t 


+ 


) 


( 3 ) 

(4) 

( 5 ) 


Here U is the velocity of the shock front and AH is the gain in specific enthalpy in 
crossing the shock front. Quantities ahead of the shock front are given by the 
subscript 0 and no subscript is used for the quantities across the front. On applying 
the hydrodynamical operator. 


dx 


_ _9. 

dx 


+ 


l 

u 


JL 

dt 


to eqn. (3) one gets, 


m 

dt 


+ U 


du 

dx 


K 

Pn 


1L 

dx 


K. 

Po U" 


JA = o, 

dt 


where 


k = P a U JPU = 1 - JL 


dp 


U 


d\J 

dp 


( 6 ) 


( 7 ) 


( 8 ) 


A fourth equation is necessary in addition to the three equations (1), (2) and 
(7) given above, to solve for the four partial deriatives of u and p. For this, purpose 
consider the work done per unit area of the generating surface (including the loss 
of energy by radiation), which can be written as 

x x 

a * = f p ° f (e + w ) - f p ° F + 


if «'(p. +p')dt, 

*9 (*> 

t 1?9 | ; 


( 9 ) 



where E and W are the specific internal energy increment of the fluid and the 
radiation energy density increment at P 0 . u' and p' denote the particle valocity and 
excess pressure behind the shock front. f 0 («) is the time of arrival of the shock at *, 
a a is the Lagrangian coordinate of the generating surface of the shockwave. F 
is the net energy loss by radiation per unit mass of the fluid. Introducing, 

h “ E + -r- + P o A (VP), (10) 

* O 

the specific enthalpy increment of an element of fluid on crossing the shock front 
and returning to pressure ? 0 and then comparing it with the equation (9) for the 
limit #-»co , we get 

CO 

D (x) = f fu'p’dt (10) 

l o («) 

=J Po lo "hdt 0 -j r 0 z 0 a ¥d £ 0 (12) 

X X 

In the Brinkley-Kirkwood theory of propagation of shockwaves, increase in en- 
thalpy of the gas element is equal to the mechanical work done by the shock wave. 
However, in astrophysical cases where high temperatures are involved some of the 
energy is radiated away (which escape into intersteller space). 'Consequently the 
difference of enthalpy increment and this energy loss by radiation should be equated 
to the mechanical work done by the shock wave. 

As in our previous work, from (11) (Bhatnagar and Kushwaha 1961 a) one 
can get 


1 x 1 dt) . _ X a puv (x) 

u IT p 8* X D (*) 1 W 

where v(x) is a function depending on the form of the wave under consi- 
deration. 


By eliminating JUL and from eqs. (1), (2), (7) and (13) and solving for 
0* 0^ 


.0 ft 


and we get the equation 
0 * 0 / 


<7 )+ 2 (‘-£) g 


dp x a p* l G otfi 

v * P 0 U* "2l(i+"l) - G 2 x 2(i+l) -G ,(14) 

where, 




( 15 ) 
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From (12) one gcti 


dD 

dx 


x* P 0 h + x a P 0 F. 


( 16 ) 


In the present work, only the case of a plane shock wave, i.e, a =- 0 is consi- 
dered. It has been shown earlier (Bhatnagar and Kushwaha, 1961 b) that the 
introduction of spherical curvature does not introduce any material difference in 
the computed velocities. If one makes a = 0, then eqs. (14) and (1 6) reduce to 

d P __ _ y G (17) 

dx D(*) * P 0 XJ*~ ‘ 2 (ft + 1) - U ’ 

and 

^ = _ PoA+ P 0 F. ( 18 ) 

So far these equations are exact. For application to astrophysical^ problems 
these equations have to be numerically integrated by chosing appropriate initial 
values for the problem in hand. 


APPROXIMATIONS 

The Brinkley and Kirkwood (1947) theory of decay of a shock wave makes use 
of the fact that the irreversible processes all occur at the discontinuous shock front, 
so that these effects can be estimated purely from the strength of the shock, together 
with a simple assumption about the return of the material behind the front to 
the original equilibrium state. On the other hand, the energy carried by the 
disturbance can also be expressed in terms of the shock strength, together with an 
ass umption about the form of the wave, so in this way the dissipation can be related 
to the energy transporated. Thus a combination of thermodynamic reasoning 
with fairly plausible assumptions circumvents the necessity of explicitly working out 
the details of the dissipative processes in the shock transition zone. This method 
was applied to gas dynamic waves in the chromosphere by Schatman (1949) and 
Wevmann (i960) and to the /? -cephei radial velocity variations by Odgers and 
Kushwaha, (1959) and later by Bhatnagar and Kushwaha (1961 a and b). The 
Rankine -Hugoniot relations are applied to connect the physical quantities across 
the shock front with those ahead of it. To obtain the physical quantities behind 
the shock from at an instant when the decay starts one has to make some assumption 
about the path of the gas element in the P — P plane. In our earlier works (1961 
a and b) we assumed that the mechanical decay started after the instaneous equali- 
sation of temperature on the two sides of the shock front by large and fast radiation 
flux, but did not account for this flux in the relevant enery equation. Later (1962) 
an attempt was made by the same authors to include the loss of energy by radiation 
in our equation by neglecting the radiation pressure and energy density. It was 
done for the following two cases on the assumption that after crossing the shock 
front the radiative loss of energy takes place at constant volume till, (i) the gas 
element returns to the original adiabate, which consequently means a return to the 
: original entropy (as the changes in ionisation behind the front were not considered) 
(ii) the temperature across the shock front becomes equal to that ahead of it. 
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The mechanical decay occurring after this was assumed to be due to adiabatic 
expansion to reach the original pressure. We have used the same assumptions in 
the present note also to estimate the loss of energy by radiation flux. We have taken 
the stellar matetial to be a mixture of a perfect gas and radiation in equilibrium. 
We can express the various physical quantities in terms, of the shock strength, 
defined by the relation: 

P- P o (19) 

The relation between the densities on the two sides of the shock front can be 
found as follows. If the expression for enthalpy increment for a mixture of perfect 
gas and radiation, 

A H = — ( Z# list ) + 4 ( — - — 0 ) (20) 

y - i v. p p 0 / v p c o' 

is equated with the corresponding expression in eqn. (5) then one gets 


7 = fo _ M + s 

^ ~ P (N 0 + 8) U + *) - t ’ 

"(2i) 

where, 


M-2[( ^-4)^ + 4], 

(22) 

N = 2 ( r ~-i - 4 )’ 

, (23) 

f> - W . Po - P8o/Po. 

'• (24) 


,(25) 

From eqs (24) and (25) one can easily get 


1 «“ P _ 1 — i®o Cl _L_ ^3 74 

P /3 0 4 ( + } ' ■ 

(26) 


which gives the dependence of /? on s and Z as shock decays. 

Using the Rankine-Hugoniot relations (3) and (4) one can write U , G, k t t 
and c Q as follows : 




r - £ 


d A. 

ds 


]• 


(30) 



(31) 


‘ a = ■ • . eo a .20 +,), 

1 O 


(32) 


r -■ P +- 


(4-3 j8)» (Y - 1) 

/?+ 12(7 - 1) (1 - /?)’ 


(33) 


where, 


+ 


(4 - 3 / 3 n ) 2 (7 - 1 ) 

0O + 12 (7 - 1) (1 - J3 0 ) ’ 


^ _ 'N /? + 7) 1 - M) 4- 1 
ds [(N 0 + 8) (1 + s) -<]* 


(34) 


( 35 ) 


_ N /3 (1 - /5) (1 + f) [(N jff + 7) (3 i + 4 - M) + 4] + 3 (M + s) 

[N (ft + 8) (1 + j) - ip ‘ (1 + ,) (24 p - N j8> - 32) + (4 - 3 j8) t 

We can determine the energy loss by radiation for following two cases. 

{ Case A : 

On crossing the shock front the steller material radiates at constant volume 
till the entropy on the two sides has the same value Le. it returns to the original 
adiabate specified by T 0 and thereafter it expands adiabatically to the original 
pressure. 

Case B : 

The radiative loss of energy behind the shock front takes place at constant 
volume till tlie temperature on the two sides of the front is equalised and thereafter 
p ex P£ nc > 3 adiabatically to the original pressure. We first evaluate in detail h and 
F for the case A and then give the corresponding expressions for the case B. 


CASE A 

crossing the shock front the internal energy of the gas and the radiation 
and g fiti Slt rl lncrea j?\ Then the material radiates some energy F per unit mass 
Therefore, W ° rk ^ e *P andin S back t0 the ^iginal pressure! 


p = A' - A w 


(36) 


masJoncrossin^ft? 3 ^ 11 ^ 11 * nterna l energy and the radiation energy per unit 

S e shock front. A to is the work done per unit mass of the stellar 
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matter in adiabatic expansion to the original equilibrium pressure. The terms on 
the right side of eq. (36) can be expressed. 


A V - (l + ») . 1 ■ [ f-fr + 3 <' -« ] - -%-[ Tfr , +3 c - w ] • < 37 > 

and 

a»--% ■ + - «]— $-[ A + 3(1 - 

Therefore, 

F=i- (!+<)«.[ 5^! +3 0 -«] - r " t ^,+ 3 0-W] • (*>> 

The enthalpy increment can be written as 

A = A H + A H* — 

^ ' r ° [ ^ + 4 (1 - ^ : 

+ -^[^---1 &>+4(l - &)], (40) 

where, A H* is the change in enthalpy during the adiabatic expansion to the final 
pressure P Q . - — 

Using eqs. (27) to (40) in eqs. (17) and (18) we get 

-$-*= - - - - P 0 . . (1 - Z) • LiL + fl n .. _ j SL" X *_ , (41) 

dx X) (^) >■ ^ i 

r ( i + s)[ 2(i - Z) - * •-=% ] + Z* 

and 

^ = Po •[ - ~0 + 2) - ~--K ^ + 3}+<e 1 “ r °] .(42) 

It is convenient to introduce in eqs. (41) and (42) two dimensionless quantities 
£ and y defined by the relations 

* = R._>> . (43) 

and “ 

DW=P„.R.{, (44) 
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where Ris. the stellars radius. In terms of these, equations (41) and (42) re- 
duce to 


di 

dy 


and 




r (i +»> ( i - - £jL 

dZ 


(45) 


r (1 + ,) [ 2 (1 - Z) - < • -f ] + Zs 


dy 


+^) ~ + Z (1 + P + 3 ] + ^ " l ° • (46) 


CASE B 

Following the similar procedure we find that in this case, A® remains the 
same, but &w becomes 


A to 



a,, r ^ 

ft L V - 


f3 0~ft) ] . [ 1 - ( -|^X Z ) l ' 1/Tl ]• (47) 


where 


and 


ft = 



ft 

ft d — Z) + z 

B . (4 - 3 ftv> (y - \) 

^ /?, + 12 (r - 1) (1 - j8 a ) ' 


(46) 


(49) 


Therefore, 


^* F= _Po_ 


d + O.'^C^y+Ml - |51 ]--A.[ ^r+3(l - ftV] 


- A. 


+ 3(1 -^ ] • [ 1 - ( -f: - ^ ] • 


A 


vi-1 /ft 


(50) 


and 


T-i.V|,i+o 


So 


s 0 ' ft L r- 1 


[ ft + 4(l-ft)1. 


1 '*• 


-t)‘ l,r *]. (51) 
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The eq*. (41) and consequently (45) remain the same but eqn. (42) be- 


comes 

dD 


dx 


+ «( 1 + .){^ I +3(l_ W }_{& 1+3( ,_ w }] i 


which can be written in terms of £ and y as 


d£ m _ £ 




(52) 


4 * ( ! + s) [ y~Y + 3 (1 - p) ] - [ £a- +3 (i-/? 0 ) ] (53) 


INITIAL VALUES AND NUMERICAL INTEGRATIONS 

For application to actual problems one has to integrate numeracally the 
equations (45), (46) for the case A or (45) and (53) for the case B. In the present note 
we compute theoretically the radial velocity variations of (3 -cephei star BW Vul. 
For this star we take y = 5/3. The initial values /3 0 = 0-8, T 0 = 1-5 and £ = 1 
are as before (Bhatnagar and Kushwaha, 1961 a), v is a function depending on the 
form of the wave under consideration. Its initial value varies between 2/3 and 1. 
For asymptotic behaviour its value was given by Brinkley and Kirkwood to be 2/3 
which has been adapted here. Using the initial value of particle velocity - 130 Km/ 
sec. (Odgers and Kushwaha, 1959), we solve equations (21), (26) and (27) simultan- 
eously to determine the initial values of s, Z and /?. These give s ~ 48-3; p/p o - 
5 - 5340 and /3 =* 0*3228. The equations (45), (46) for case A or (45) and (53) for t>he 
case B, have been numerically integrated following Harm and Schwarzschild 
(1955) scheme At every step and Z are determined from eqs. (21) and (26) for the 
value of s. 

DISCUSSION AND CONCLUSION 

The value of £, s } uJG 0 and UJC 0 obtained fiom the integrations as the 
shock propagates are given in tables 1 and 2 for the two cases under considera- 
tion. 
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TABLE 1 


Decay of a shock wave for case A 


y 

C 

$ 

u/C, 

u / c 0 

0*00 

1*0000 

48*3000 

5 * 1 1 

6*24 

0-01 

0*9384 

43*7086 

4*85 

5*94 

0*02 

0*8830 

39*7272 

4*62 

5*68 

0-03 

0*8331 

36*2546 

4*40 

5*44 

0-04 

0*7880 

33*2128 

4*21 

5*21 

0-05 

0*7471 

30*5338 

4*02 

5*01 

0 - 06 . 

0*7097 

28*1654 

3*86 

4*82 

0*07 

0*6756 

26*0613 

3*70 

4*65 

0*08 

0*6444 

24*1867 

3*56 

. .. 4*49 

0*09 

0*6158 

22*5096 

3*42 

4*34 - 

0*10 

0*5895 

21*0036 

3*30 

4*20 

0*11 

0*5653 

19*6472 

3*18 

4*07 ' 

0*12 

0*5429 

18*4222 

3*08 

' ' 3 "95 

0*13 - 

0*5222 

17*3121 

2*97 

r 3*84 

0*14 

0*5030 

16*3037 

2*88 

3*74 

0*15 

0*4852 

15*3850 

2*79 

3*64 , 

0*16 

0*4686 

14*5460 

2*71 

3*55 

0*18 

0*4388 

13*0744 

2*55 

3*38 

0 ; 20 

0*4127 

1 1*8306 

2*41 

3*23 
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fAfeLfe 1 

Decay of a shock wave for Gas# B 


y 

; 

s 

»/c 0 

u/c 0 

o-oo 

1-0000 

48-3000 

5-11 - 

6-24 

o-oi 

0-9282 

43-6858 

4-85 

5-94 

0-02 

0-8629 

39-6461 

4-61 

5-67 

0-03 

0-8033 

36-0941 

4-39 

5-43 

0-04 

0-7488 

32-9575 

4-19 

5-19 

0-05 

0-6988 

30-1732 

4-00 

4*98 

0-06 

0*6528 

27-6977 

3-82 

4-78 

0-07 

0*6104 

25-4851 

3-66 

4-60 

0-08 

0-5712 

23-5018 

3-50 

4-43 

0-09 

0-5350 

21-7188 

3-35 

4*27 

0*10 

0-5013 

20-1107 

3*22 

4-12 

(HI 

0-4701 

18-6567 

3-10 

3-98 

0*12 

0-4410 

17-3384 

2-98 

3-84 

0*13 

0-4139 

16-1402 

2-86 

3-72 

0-14 

0-3887 

15-0486 

2-76 

3-60 

0-16 

0-3429 

13-1396 

2*56 

3-39 

0*18 

0-3027 

11*5340 

2-38 

3-20 

0-20 

0-2673 

10-1738 

2-22 

3-02 


If one plots U/C 0 against y y the relation is found to be almost linear for 
C 0*15 and wc can relate UJC 9 and y by a linear equation, which on integration 
gives time, Th« values of computed velocities as a /unction of time are given in 
Table 3* 
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Gase A 


TAbLe 3 


Case 


t (days) 


u (km/sec) 


t (days) 


u v km/sec) 


0 


111-52 


0 

111-82 

0-004270 

106-81 

0-008654 

102-20 

0-013159 

97-94 

0-017795 

94-01 

0-0225.67 

90-36 

0-027483 

86-97 

0-032553 

83-82 

0-037787 

80-88 

0-043196 

78-13 

0-048791 

75-55 

0-054587 

73-14 

0-060597 

70-88 


0-005313 

106-38 

0-010863 

101-56 

0-016671 

97-08 

0-022762 

92-90 

0-029165 

88-99 

0-035913 

85-19 

0-043047 

81-89 

0-050613 

78-66 

0-058666 

75-62 


In fig. 1 we have plotted the observed radial velocity from Ojgors (1955) ob- 
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stations. It is seen that for this particular cycle the initial plotted value of 
velocity after applying corrections for fore-shortening and limb-darkening is of the 
order of - 112 km/sec. Hence we have shifted the time scale of our computed 
velocity curve in such a manner that time scale starts at a velocity ^ — 112 km/sec. 
We 6nd that the computed radial velocity curves are some what lower than the 
observed one and time of decay is a little larger than the observed value. For this 
a qualitative explanation can be given as follows. We have supposed the shock 
front to be transparent to the radiation below so that the observations refer to some 
sort of integrated motion of the matter behind the front. After an element of gas 
has received its initial outward velocity it then moves under elective gravity (not 
negligible in these stars but not included in our considerations). In the lower layers 
of the atmosphere where the shock wave is strong the instantaneous velocity im- 
parted by the shock dominates the observed motions, since the even lower layers 
are largely obscurred by the much increased opacity. However, at higher levels 
where the energy of the wave has largely been dissipated the observed motion con- 
sists mostly of material under the effective gravity forces only. On this view as 
the shock progresses, it should rapidly decay due to gravitational forces. The 
computed velocities which are having large negative value will be reduced and 
may give the required agreement. We have neglected the variations of pressure 
and density with space in the undisturbed atmosphere. As the total distance tra- 
versed by the shock wave is large, this inhomogenity may also have some effect on 
the computed rate of decay, The study of the effect of inhomogenity is in hand. 
Finally we have estimated the loss of energy by radiation just intuitively without 
giving any specific mechanism for it. In the previous work (-1962) the initial jumps 
in temperature and density across the front were T/T 0 1*6 and P/P 0 £z3*8. In the 
present case they are changed to T/T 0 ^ 3*6 and pjp Q zz 5*5. This change is due 
to the effect of radiation pressure and density being taken into account in the 
present work. These changes in jumps are in closer agreement with the ideas of 
Abt (1955) and Save doff (1954). In figs. (2), (3) and (4) we have plotted P/P 



Fig. 2. Pressure versus distance curve across the shoek front. X give the case A and 0 give Case B. 
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pjp a and T/T s against y (which have bean tabulated in Tables 4) to study roughly 
the variations of pressure, temperature and density across the shock front with space. 


TABLE 4 


Case B 

Case A 

y 

PlP 0 

H 

o 

P/P 0 

p\p B 

T/T 0 

P/Po 

0*00 

5-5340 

3-5946 

49-30 

5-5340 

3*5946 

49-30 

0*01 

5-4434 

3-4846 

44-69 

5-4437 

3-4854 

44-71 

0*02 

5-3505 

3-3824 

40-65 

5-3533 

3-3836 

40-73 

0-03 

5-2576 

3-2842 

37-09 " 

5-2632 

3-2888 

37-25 

0-04 

5-1645 

3-1918 

33-96 

5-1733 

3-2000 

34-21 

0-05 

5-0713 

3-1052 

31-17 

5-0839 

3-1161 

31-53 

0*06 

4-9788 

3-0225 

28-70 

4-9975 

3-0369 

29-17 

0*07 

4-8862 

2-9432 

26-49 

4-9116 

2-9629 

27-06 

0-08 

4-7943 

2-8673 

24-50 

4-8263 

2-8938 

25-19 

0-10 

4-7030 

2-7948 

22*72 

4-7438 

2-8273 

23-51 

0*10 

4-6126 

2-7266 

21-11 

4-6642 

2-7645 

22-00 

o-n 

4-5230 

2-6609 

19-66 

4-5851 

2-7059 

20-65 

0-12 

4-4320 

2-5985 

18-34 

4-5086 

2-6499 

19-42 

0-13 

4-3465 

2-5394 

17-14 

4-4326 

2-5975 

18-31 

0-14 

4-2600 

2-4818 

16-05 

4-3592 

2-5474 

17-30 

0-16 

4-0910 

2-3745 

14-14 

4-2176 

2-4549 

15-55 

0-18 

3-9263 

2-2778 

12-53 

4-0850 

2-3713 

14-07 

0-20 

3-7680 

2-1875 

11-17 

3-9588 

2-2963 

12-83 


It may be noted that there is not much difference in the computed values for cases 
A and B. Density and temperature gradients are almost linear. 
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ABSTRACT 

The problems on General Uniqueness and successive approximations was considered by F. Brauer 
(1959) in a sufficiently general way so as to include the results of F. Brauer and S. Sternberg (1958) 
which are generalisations of Coddington and Levinson and E. A. Kamke. In this paper we shall prove 
a pair of theorems under conditions which are generalisations again. 

§1. We know that the uniqueness of solutions and the convergence of succes- 
sive approximations are logically independent in the theory of Differential equations 
cf. [2] and [6]. In this paper, we shall prove a pair of theorems generalising the 
results of [2] and [5] which are generalisation of Coddington and Levinson [3] and 
Kamke [4]. The generalisations here are in a different approach as pointed out 
in [1]. 

§2, Consider the differential system 

(1) x'(t) =* /(* (<), t) *(<>)-• ( ' ” ) 

Where x and / are rt-dimensional vectors. 

Let $ ( x(t ), Q be a function, defined for vectors x and real t with non-negative 
real values, which is continuous in (#, t) has one sided partial derivatives with 
respect to t and the components of a? such that $ (x, t) = o implies x — o> Let 
<ji X be the partial derivatives of $ with respect to t and x respectively. 

Let u>| (r, t) (z = 1,2) be continuous, non-negative functions defined for 

o < t < a, r ^ o which are monotone non-decreasing in r for each fixed t. We 
will always assume that / satisfies the pair of conditions : 

4>t ^ (*i — j 1 1* 4- (*i - ^ (#i *— *a) > f I" 

(2) [ /(*!> 0 - /(**, t) ] < «■ [ <t> { (*i - *a)“.<}.']i = l,2 

where o > 1 for (< x x x ) and (<, x 2 ) in the region o<t< a, \ x\ <b 

TheoremA : — Let J(x,t) be "continuous and satisfy (2) in the region o < t < a, 

| * | ^ b. Suppose A(f) and B(/) are functions on o ^ t < a with A(o) = B(o) = O, 
such that 

lim A(t) _ q 
t-*o B(0 “ 


( 3 ) 
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Suppose also that all solutions u(t) of 

(4) m'CO = 0 

with «(0)--0 obey u(t) < A(<) ono</<4 and that the only solution 
v(t) of 

(5) v'(t) = o, 2 (v(t), 0 

on « ^ t < a such that 

lim p(Q q 

B(0 

is the trivial solution, Then there is at most one solution of (1) on o ^ . t < a 

Proof Suppose ^(f) and # a (0 be two solutions of (I) on o ^ t a, if 
possible, let 

/*(*) = #i(0 “ *a(0> so that 

(6) ‘ - ; /*'(*) = V(0 - *s'(0 

Let m(t) --■= <p ( p* (0, i ) , m*(t) = lim Sup [ m(t) - m(< - A)_ ] „ 
x ' / z — ^0 

then 

^-J). * [ <f> ( > * (0. 0 - A),*.) 

-) h $ ( (< - A), < ) -$(/*“(*- A),* - A ) ] /A 

. - Since $ has one sided partial derivatives, the left hand side pf the above 
equation is bounded by a sum of the form 

{« “ 1 (0 . <p x (V ( t ), t ) + * x )[ M(0 - K‘ - *) ]/ h 

+ ^ ^ (t), / ^ 

where s 9 s -> q as h v 

Now allowing A —> O, we obtain in view of (2) and (6) that 

(7) m*0> < «i (m(0. 

Let us assume that there exists a £, o < g < a such that 

(8) K£) > A (|) 
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. . ^® n t * iere 7S a solution u(i) of (4) passing through the point (£, m(£) ) and 
existing on some interval to the left off As far to the left of £ as Vt} exists it 
satisfies w 

(9) «(<) < m(0 

To prove (9), consider the solution u(t, e) of the equation 

(19) u' = Oij («,-<) + s «(£) = m(|) 

which we know exists for sufficiently small e > o, as far to the left off as alt) exists. 
Moreover 

lim u(t, e ) = u(t) [4] 

£ ~* 0 

Thug to prove (9), it suffies to establish that 

(11) «(/, e) < m(t) 

for all e > o and all solutions of (10). If this were false, there must be a least upper 
bound V of numbers t < $ for which (1 1) is not true. Since 'm(f) = «(£) = e ) 
and the functions m(c), u(t, e) are continuous, it follows that 

(12) m(V) = u(v, e) - m*(y) > u'(v, s) 

Then 

“lOb m (V) ) + e — “iC 7 ?. «(7, «) ) + £ 

— U'(V, £ ) 

< m*(y) 

< m(y) ) 

using (7), (10), (12), This contradiction proves (11), which implies(9). 

The solution a (t) can be continued to t = O. If u(c) = o for some e, o<c<£, 
we can effect the continuation by defining a(<) = o for o < / < c ; otherwise (9) 
ensures the possibility of continuation. Since m{ O) = O, lim u(t) = c and we 

t — * o 

define u(o ) - O. Now we have a solution u(t) of (4) with u(o) = O and by hy- 
pothesis u(t) < A(t). However, in view of (8) and the definition of u(A, W e have 
> A(|), a contradiction which proves 

(13) "»(0 < A (0 o < t < a 

To prove the uniqueness of solutions, we must show that m(t) vanishes identi- 
cally on o < t < a. To complete the proof, we must show that this implied by (13). 
Proceeding as before. but using <a 2 in.the place of «„ we obtain . . 

m*(t) < (m(0, t) 

The assumption that m{ f) > a for some F, o < T < a yields by the same argument 
as before, a solution v(t) of (5) on o <1 t ^ r such that r(T) = m(r), o ^ v(t) ^ m(t), 

■ - ' . . v(o) ■=■ 0 
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fhcfr 

O <; hm v(i) . lim w (f) iim A (i) _ ^ 
t -* o B(<) ^ t o B(ty % t -> o B(0 ~ 

using (3) and (13;. But by hypothesis, this implies that v(t) is identically zero, which 
contradiction v(T) = w(T) > o and therefore m(t) vanishes identically o no (I < a 
which completes the proof of the theorem since $(*, t) = o implies x »= o 
Suppose, if wc take 

(r, t) = a» (r,t) 

(r, 0 = 2 M 

where M is a bound for | /(#,<) | in o < * < a, | a; | < b 

we gef the^^sults'off^ = * ^ let be 3 solution of ( 4 ) with A(o)=A'(o)=o, 

of,JiirT^,?,Uqui«on i LfS', “_ d I.' S, “’ , “ 8 ***** “ “ rli " " f ° r 8 SyS " m 

thcor.ifofKamke 8 4 |. ° “ 1 * 1 ' m obtai ” ,he r “ ulK of g“« a > uniqueness 
§3. The successive approximations to the solution of (1) are defined by 


(14) 


£ 

* j+ 1 (0 =/ /(*■ (S), S) d s 


*o (0 = 0 


(.j o> 1>*»«) 

guarantee^th^™ Section ’ f we shal ! im P 0Se two additional conditions on $ to 
5S approrin.auons. Wc a 5!u „c ,ha, A, 

# 1 ( s > [ / K-+1 (S)> S) -y (*.: (S), S) ] d s-t ) 

t-h J 7 

... * t 

(15a) </$, C^(S),S)dS+ / (S)$,(/*J(S) > S)[y(* J+1 (s )>S ) 

■ • “/(*,■ (S),(S) ] dS 

* («« (h* '»), Xj+ l + X. < $ (*. + 1 , f,) + $ («. , h ) " 

Where a ^ .1 for any continuous functions of f and m . — x- .. 

>■ j+i *y 

Cl 5 b' 2 f\i, L f et f- be as *? §, 2 ’ w . ith the additional properties (15a) and 

(15b;. Let f(x, t ) be continuous m the region o < t < a , | a | < t and bounded in 
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and 

(15b) 



n birth by M in this region. Suppose the hypothesis of theorem 1 are satisfied, then 
the successive approximations (14) converage uniformly on the interval O ^ t < ft. 
where (3 — min (a, bf M) to the unique solution of(l). 


Proof : — It is easy to verify that since / is continuous and hence bounded, the 
sequence^. (/). of successive approximations is uniformly bounded and equi-contin- 

uous in the Euclidean norm on some interval. It follows that there is a subsequence 
x j k ^ w hich converges uniformly on this interval to the solution x (0* Since 

t 


7*+l 


(0 = / f(* jk (S, S)dS 


Xj (0 converges uniformly to a solution #*(*). We shall prove 

(16) <K* j+1 (0 - (0.0-»o 


pn this interval which imply #(£). = #*(*) so that x(t) is a solution of (1). Since the 
solution is unique, every convergent subsequence converges to x(t) and it follows that 
the original sequence converges to x(t). Since this sequence is uniformly bounded 
and equi-continuous, the convergence is uniform in the Euclidean norm which 
implies by the continuity of that $ (#(0 — (0> 0 converges uniformly to zero 

on some interval. 

To prove (16} s since 

(i t ) = xj j ^ l (t) ~ xj (/)> therefore 

and let m(t) = lim Sup <£ as j co 

J 

then m(o) = O and m(t) is continuous, since it is the upper limit of a uniformly 
bounded cqui-continuous sequence of functions. We have ' : ’ **" v 

t 

t — h 
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~ ' ' i . 

= f « **“+ 1 (Sj [ /(* i+1 (S), S) - /(*. (S), S) ] S 

/-A 

' " ; 

••• - * (*J +l (0. - /*; + L (i - A), 0 - $ ( f « <*“.“} (S) [ /* >+ J (S), S) 

£ — h 

-f(*j (S),(S) JdS.t ) 

£ £ 

< J $, (S), S) rf s + J « 1 (S) ^ (/*“ (S), S) [ f(x. +l (S), s, 

f t - h 

-f(x.(S),S)]dS 

t 

(17) < ( m l [ $ (S), S), S ] d s fr«m (15a) and (2) 

* - h 

given any S > O, there exists an integer N($) independent of s and j such that 

(18) $ ■'/»“ (S), S) < m(S) + 8 > N(S) 

This follows from the fact that m is uniformly continuous and'/ A J is equi-con- 

tinuous. Since Ml ( r , t ) is assumed monotone non-decreasing in r, it follows from 
(18) that 

( t 

f », [ $<m“ (S), S), S ] d S < f ai [»(S) + 8, S ] d S 
i - h t ~ h 

using (17), we obtain 

£ 

(191 $ (M* + j (/) - 1^ + !«-*),*) < J" [ W (S) + 8, S ] d S | > N (s) 

i — A 
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(20) 


From the definition of m(t ) and (19) and (15b), it is easy to verify that.. 

t 

m(t) - m(t -A) s? f w , [ w(S) + 5 , S } d S 
t-h 

Since ^ (r, t) is continuous in r 


“1 £ ffj(S) + 8, s ] ->#,[ m( S), S J as 8 = O 

and this together with ( 20 ) yields 

t 

(21) m{t ) - m(t - h) sg; J'a, (m(S), S) dS 

t-h 


This follows that 


m*{t) = lim Sup f m(t) - m{t - h) 1 ,, 
h ->o L J 


"1 [ m(0. ' ] 


The argument used in the proof of theorem l starting with ( 7 ) proves 
( 22 ) m(t) sg A (t) 


To complete the proof, we must show that (22) implies m(t) == O. The argu- 
ment is much the same as the last stage of the proof of theorem 1 , Suppose r»(r)>o 
for some I\ O < T < /3. A repetition of the first part of the proof using in th e 
place of <* v gives a solution v(t) of (5) on O < t ^ T such that r(T) = m(r), O ^ 
•(0 < m(t), v(o) = O 


Then 


o< 


lim 

t 


v{t) 

*B(0 


lim 
t 0 


m(0_ 

B(0 


^ lim 
t -> 0 


A(0 

B(r) 



using (3) and (22). By hypothesis this implies that z;(f) vanishes identically, con- 
tradicting o(T) = m(T ) > o and it implies m(t) = O 9 o ^ t < /? which proves the 
theorem. 


Suppose if we take 

®i 0 = « (r, 0 
(r, 0 = 2 M 
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where M is a bound for | /( x, t) | in O < t < c, | x | ^ b 

Further if we put B(i) = t and let A {t) be a solution of (4) with A(o) =» A'(o) 
and a = 1 we get the results of [2]. 

Again putting t) = | x | , we obtain the successive approximations theo- 
rem of Coddington and Levinson [3]. 

JVote . — For « => 1, the above theorems were pointed by [1] but our results are 
more general in the sense that a takes any positive real value a ^ 1. 

REFERENCES 

1. F. Brauer, Can. J. Math., 11, 527 (1959). 

2. F. Brauer and S. Sternberg, Amir. J. Math., 80 , 421 (1958) 

3. E. A. Coddington and N. Levinson, J. Ind. Math. Soc. , 16, 75 (1952). 

4. E. A. Kamke, Leipzig, 83 (1930). 

5. V. Lakshmikanth and M. M. Subramanyam, Proc. Nat. Acad. Sci. India, 29A 

(1960). 

6. M. Muiler, Math. Zeit., 26, 619 (1927). 


f 3612 3 



LATITUDE DEPENDENCE OF AMPLITUDE AND PHASE OF 
DIURNAL VARIATION 

By 

R. S. YADAV 

Gulmarg Research Observatory, Gulmarg, Kashmir 
[Received on 23rd October, 1962] 

ABSTRACT 

Tiie latitude dependence of amplitude and phase of the diurnal variation, of meson intensity has 
been studied by analysing the data of I. G. Y. stations, Gulmarg (\ =24’7°N), Ottawa (\rz 57ojsj) t 
ChruchiU (\ :~70°N) and Resolute [\ = 83°N) for the year 1958. 

INTRODUCTION 

For understanding the energy spectrum of the primary anisotropy responsible 
for the daily variation one of the very important tool, is the study of the latitude 
dependence of the amplitude and phase of daily variation of cosmic-ray intensity, 
Dorman has calculated the energy spectrum for the period 1937-51 using the data of 
meson intensity recorded by ionization chambers. During the recent I. G. Y. f 
Cosmic-ray meson telescopes have been operated at several stations in the world by 
various research groups. This has provided a unique opportunity for studying the 
problem of latitude dependence of the diurnal variation of meson intensity. 

Thus with, this objective in view, we analysed the pressure corrected meson 
intensity data for I. G. Y. stations Ottawa (57°N), Churchill (70°N) and Resolute 
(83° N) along with that of Gulmarg (24*7°N) for the period January to August, 1958, 
to examine the latitude dependence of the diurnal variation. 

EXPERIMENTAL RESULTS 

The values of the amplitude and phase of the diurnal component of meson 
intensity, recorded with meson cubical telescope at Gulmarg, . Ottawa, Churchill 
and Resolute, were determined with the help of harmonic analysis using 12 bihourly 
values. The results obtained are shown in the tables (!,■ 2) and the figure L 

r 2.13 i 



table i 


Geomag. Lat. 

X 

Station 

a d 

Year 

57°N 

Ottawa 

0-39 ± 0*008 

1958 

70°N 

Churchill 

0*25 ± 0*008 


83°N 

Resolute 

0*06 ± 0*008 

it 


TABLE 2 


Geomag. Lat. 

A 

Station 


Year 

24*7°N 

Gulmarg 

212° 

1958 

57°N 

Ottawa 

194° 

>» 

70°N 

Churchill 

O 

CO 

•'T 

l-H 

»} 

83°N 

Resolute 

124° 

)» 
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blSCbSSION AND CONCLUSION 

Thompson (1938), and Thambyahpillai and Elliot (1953) studied the latitude * 
dependence of the diurnal variation of cosmic- ray intensity, but they found no 
significant dependence on latitude. Firor et al (1954) have also examined the daily 
variation of the nucleonic component at two different latitudes viz., Huancayo 
(A =0 3 ), Climax (\ = 48°N) arid found that the ratio of the peak to peak variation 
at these two stations is 1*4 ± 0*02. The analysis of the data of Huancayo (A =0°), 
Christchurch (a = 48°S) and Godhavan (a ~80°N) where the hard component data 
have been averaged over eight years, brings out the following facts (Dorman 1957) : 

(1) That the variations at Cheltenham and Christchurch which are sym- 
metrically situated with respect to equator are similar and the hour of maximum is 
also approximately the same. 

(2) That the amplitude of variation at Godhavan (A » 80°N) which is 
situated near the pole is smaller than at the other three stations and 

(3) that the diurnal variations are apparently of irregular character. 

The study of the amplitude of the diurnal variation of the meson intensity 
recorded with meson cubical telescope at Ottawa (A ^ 50°N), Churchill (A = 70°N) 
and Resolute (A =*83°N) reveals that it does not remain constant over the various 
latitudes considered. It is apparent from the table (1) that the amplitude does not 
show any systematic variation with latitude. However, we see for the station 
situated near the pole (a = 83°N), the amplitude is appreciably smaller than others. 

Leaving Gulmarg (9,000 ft.), the other stations of the table (2) are situated at 
sea level. Sarabhai et al. (1955) have shown that the phase of the diurnal 
component is earlier at high altitude for the same latitude. It means that the phase 
for a station which would be situated at the same latitude as Gulmarg but sea-level 
will be somewhat greater than what has been shown in the table (2). So our 
conclusions regarding the latitude dependence of phase, will not be effected because 
of the difference of altitude of Gulmarg from other stations. 

Form table (2) and the figure (1) it appears that the phase of the diurnal 
component of meson intensity is earlier at higher geomagnetic latitude or we may 
say that the phase of the diurnal component decreases as latitude in creases. 
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SUPERPOSA3ILITY AND STREAMLINE PATTERN 
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§ 0. ABSTRACT 

Kapur [2] has deduced the relations existing between two (a) Plane (6) Axisym metrical (Poloidal) 
Self-superposable and mutually superposable flows. This paper examines the stream patterns for such 
flows. They are either circles or parallel straight lines, when (i) the fluid is non-viscous and iao-vels 
coincide with stream lines* («) the flow is viscous and self-superposable. (iii) Self -superposable flows 
which are mutually superposable. 

The steady self-superposable poloidal flows, which are mutually svperposable are axial. 

§ 1. (A). We shall obtain possible stream patterns for plane flows whose iso- 

vels and iso-curls coincide with streamlines, i. e. 

(1) q=q(P) and (2) £=£(*) 

The condition (i) is necessary and sufficient for stream lines to be parallel curves 
[£]. Also condition (2) is the same when the flow is self-superposable [2] or motion to 
be steady when the fluid is non-viscous [6]. The curvature of the stream line is 


( 3 ) 



grad ft 
| grad $ | 


r 1 0/0 Z 0/0 


dq (jp 


This shows that the streamlines are concentric circles. Also proceeding as 
Ghildyal [1] from (3) we get 

either 4'n = o or £(f) = q $) - ~ 

either way it may be observed that the streamlines are parallel straight lines. Hence 
'Vtfc have theorem 
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1. Theorem : The stream patterns for a *iass of plane rotational flows whose 
iso-curls and iso-vels coincide with streamlines arc a system of concentric circles or 
a set of parallel straight lines 

(B). We shall examine the nature of the flow whose iso-curls coincide with 
streamlines and velocity c q’ is constant throughout the region, 
i.e. (4) £ = £ 0/0 and (5) q = constant. 

In this case (3) becomes 

(6) P-i = C (*) g- 1 

Again proceeding as Ghildyal [l] we have either = o or £=0 
In either case motion is irrotational. Hence we have 

// Theorem : There can be no plane rotational flow, whose velocity ‘q’ is 
constant in magnitude throughout the region of flow. 

(G). We shall examine the stream patterns for a self-superposable rota- 
tional, viscous steady flows. It has to satisfy (2) in addition to integrability 
condition. 

(7) V 2 c m = V^-=o or (8) £" (f) P f £' (W C « = o 

(Therefore by theorm II { q’ is constant on w . All the conditions of theorem I are 
satisfied. Hence we have 

“All rotational self-superposable flows when the fluid is viscous and motion 

steady have the stream patterns as stated in Theorem I.” 

(D). Consider the stream patterns for all rotational self-superposable flows, 
which are mutually superposable. 

Kapur [2] has obtained the conditions for two suGh flows (if, £,) 
and (if 2 , (£ 2 ) as : — 

(9) £ a =£(£i) and (10) = 

(ID and (12) £,=/,(*,) 

(«) Stream patterns will be the same as stated in theorem I, when the fluid is 
viscous and motion steady. 

(j3) When the fluid is non-viscous. 
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(13) v 2 { ft (W } — H" (^) ?1 * +H' (W/i (W 

From (13) and theorem II it follows that ‘q J must be constant on^j. Again the 
stream patterns for (jp lt £ J} ) are the same as in theorem I. Also from (10) it follows 
that (^ 2 , g 2 ) will have the same patterns as ($ v g x ) 

From (a) and (/?) we can conclude : — 

“The only possible stream patterns for all rotational self-su per posable flows 
which are mutually superposable (whether the fluid is viscous or not) are the same 
as stated in theorem I.” 


§ 2. AXISYMMETRICAL (POLOIDAL) FLOWS 

(A) We shall consider self-superposable rotational steady flows. 

The self-superposability condition [2] is 

(14) 

If iso-vels coincide with streamlines we have 

(15) co q WO = + ** J 4 - | grad ^ | 

The expression for curvature becomes 

(16) - <0 { m- 1 ! m - —P } = “ m 

say which shows that the curvature varies inversely with the distance from the exis of 
rotation along the stream line. Again we can prove as in § l that no rotational 
poloidal flow exist if ‘q’ is constant throughout the region of flow. We have also 
proved [5] that if the fluid is viscous motion steady and self-superposable, the flows 
are axial. 

(B) We shall determine the stream patterns for all rotational self-superposable 
axisymmetrical (Polodial) flows which are mutually superposable. 

Kapur 1 2] has obtained the following conditions for two such poloidal flows 
Ci) and (^ a > £a) t0 be mutually superposable as 

(17) = G (^) and (18) = H (or 1 Ci) 

(19) (ft) and (20) »-i g a -/, (ft) 

(a) If the fluid is viscous, they are axial by the result stated above. 

(/?) If the fluid is non-viscous, operating E a on (17) and using (18), (19) and (20) we 
have 
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(21) H {A (A)} = G" (fc) rt + G 1 (« 

Hence ^ is constant only on « constant. Also 

(22) /> A p^ 1 = — ^ + / w" 1 ^ + constant, which is a function of ^ ne. 

velocity, pressure and or 1 £, are constant along individual streamlines, for the 
first flow which according to Prim [4] is axial. 

From (17) it follows that the flow (ft ±> £ a ) is also axial. 

Hence we have 

111 Theorem : — All rotational self-superposable steady poloidal flows which are 
mutually superposable (whether the fluid is viscous or non-viscous) are axial, 
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kinetics of oxidation of mercurous nitrate by 

POTASSIUM PERSULPHATE IN AQjUEOUS SOLUTION. 
PART I — UNCATALYSED REACTION 
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[ Received oo 19th September, 1962 ] 

ABSTRACT 

The kinetics of oxidation of mercurous nitrate in dilute nitric acid solution by potassium 
persulphate has been investigated in aqueous medium. The results are fairly reproducible when the 
necessary precautions about the purity of the materials and the medium are taken. The over- ail order 
of the reaction is one, being first order with respect to persulphate and zero order with respect to mer- 
curous ions. The reaction is attended with an induction period, the extent of which increases with 
increasing concentration of mercurous nitrate and decreasing with increasing concentration of persulphate. 
The temperature coefficient for the reaction comes out to be 3T685 for 0 01 M concentration of the 
reactants. From this the energy of activation, frequency factor and the entropy of activation have been 
calculated. There is no appreciable effect of HNO 3 on the rate of this reaction. 

INTRODUCTION 

During the systematic study of the oxidation of different reducing cations by 
persulphate ions, it was observed by us that mercurous ion is slowly oxidised by 
persulphate ion. Much work has been carried out on the kinetics of silver catalysed 
redox reactions involving persulphate ion. However the uncatalysed reactions 
involving persulphate ion have not received much attention so far. 

The present paper deals with the kinetic study of the hitherto unstudied oxi- 
dation of mercurous nitrate by potassium persulphate in the absence of the catalyst. 
This study has been undertaken to find out if this reaction follows a similar beha- 
viour as other uncatalysed reactions of persulphate ion or not. 

EXPERIMENTAL 

Potassium persulphate G. R., E. Merck, mercurous nitrate E. P., E. Merck, 
nitric acid A. R., Basynth, have been used. All other chemicals' used were of A. R., 
B. D. H. quality. Potassium iodide was of E. P., E. Merck quality. 

The standard solution of potassium persulphate was prepared fresh daily by 
dissolving the required amount of the salt in redistilled water and its strength was 
checked by the iodometric method of Szabo, Csyani and Galiba 1 , as modified by 
Khulbe and Srivastava^. The same method was used for studying the progress of 
of the reaction (as described below). 

The solution of mercurous nitrate was prepared in the minimum, concentra- 
tion of dilute nitric acid to prevent hydrolysis. The strength of mercurous nitrate 
solution was determined by titrating with standard KBr solution using bromophenol 
blue as adsorption indicator, according to the method proposed by Zombory and 
modified by KolthofF and Larson 3 . A glass stoppered Pyrex conical flask blackened 
outside with black japan and wrapped in a black cloth was used as the reaction 
vessel. The progress of the reaction was studied as follows : 5 ml. of the reaction 
mixture was pipetted out at different intervals of time and added to 5 ml. of a 10% 
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solution of sodium chloride, t’he mixture was centrifuged for about 10-15 minutes 
to remove the precipitated mercurous chloride and then 5 ml of the clear liquid 
was pipetted out from the centrifuge tube and added to 1 ml. of a mixture of 1%, 
CuS0 4 and )% FeS0 4 solution in an iodine .flask, and then a 15% solution of Kt 
was added The flask was kept for about 5-6 minutes for the complete liberation 
of iodine and titrated against a standardised sodium thiosulphate solution. The 
blank reading with l ml. CuS0 4 and FeS0 4 mixture was taken which was subtrac- 
ted from the above titre value to obtain the amount of potassium persulphate 
remaining in terms of titre value of sodium thiosulphate. 

Preliminary experiments showed that unlike Sn^-SaOg" reaction studied by 
Gupta and Srivastava 4 earlier the uncatalysed reaction is a slow reaction and has a 
measurable velocity at 40°G. and above and at OOl M or higher concentration of 
potassium persulphate. Hence the measurements were carried out at 50°G., unhss 
otherwise stated. 

RESULTS OF THE MEASUREMENTS 

First of all the reaction was repeated several times at different equimolecular 
concentration of the reactants to see if the results are reproducible or not. It was 
found that the results are fairly reproducible when necessary precautions about the 
extreme purity of the materials and the medium are taken. The following table 
gives the results in two typical runs each at equimolecular concentration of the 
reactants. 

TABLE 1 


K a S a O g = Hg 2 (NO s ) a — 0*01 M ; HNO s = 0T54 N ; Temperature = 50°C. 

Initial reading = 5‘02 ml. 


Times in minutes " 

Vol. of N/100 hypo used 

Expt. I Expt. II 

k unimolecular x 10 3 

10 

4-86 

4-91 


20 

4-72 

4-78 

2-695 

30 

4*44 

4-49 

4-479 

45 

3*95 

3-98 

6-001 

60 

3-52 

3*51 

6-716 

90 

2-76 

2-77 

7-158 

120 

2-11 

2-15 

7-511 

150 

1-63 

1-63 

7-877 

180 

1-36 

1-37 

7-506 

210 

1-07 

1-07 

7-621 


it U fn Th ^ !w v dita shows , the reproducible nature of the results obained. Further 
™ Und that the reaction has a fairly long induction period. After the induction 
period is oyer, there is a fair constancy in the unimolecular rate constant values, 
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where log c (c in terms of titre values in ml.) is plotted against time. The curves 
are nearly linear showing the unimolecularity of the reaction. Further it is 
seen that the induction period almost disappears at higher concentration of the 
reactants. 

EFFECT OF POTASSIUM PERSULPHATE CONCENTRATION 

In order to determine the effect of K 2 S a O g concentration on rate and thereby 
to determine the order of the reaction with respect to S fi O e r/ , the reaction was carried 
out at three different initial concentrations of K a S a O s and at constant concentra- 
tion of H gl (NO,) ai the results of which are represented in fig. 2. The curve of 



10 ao 40 sa 6a 70 30 30 

— time. (min) 

Fig. 2 : Variation of log c with time at 50 P C 


Gone. ofHg 2 (N0 3 )., : 0 01 M 

Cone, of K 3 S 3 O g : A 0.‘0 1 M ; B 0'015 M j C 0 025 M 
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log c vs. t in each case is linear showing the unimolecularity of the reaction. 
Further it is seen that the extent of the induction period decreases on increasing the 
concentration of potassium persulphate. The value of the specific rate constant as 
extrapolated from the curves comes out to be 8.093 x 10~ 3 min" 1 , 7T08 X 10~ 3 
min" 1 and 10*02 X 10' 3 min -1 for 0*01 M, 0*015 M and 0*025 M potassium persul- 
phate respectively. 

When the time for the same amount of potassium persulphate to decompose 
(corresponding to T5 ml. of hypo solution) at the three concentrations of potassium 
persulphate 0*01 M, 0*015 M and 0*025 M is extrapolated, it comes out to be 59, 41 
and 23*5 minutes respectively. It shows that the rate depends on the potassium 
persulphate concentration. 

The time taken for the same fraction of K a S 2 0 8 (l/4th.) to decompose as read 
out from the curves, is 53, 49*5 and 39*25 minutes for 0 01 M, 0*015 M and 0 025 M 
K.,S s 0 8 respectively which suggests that the order of the reaction is higher than 
one. However, this may be due to the fact that the extent of induction period 
depends on K a S a O a concentration. 

1 o confirm the order of the reaction with respect to S 2 O a " ion the differential 
method has been applied and for this the values of c (in terms of titre values) have 
been plotted against t (time in minutes) in fig. 3 and from the curves thus obtained 
the values o i-dcjdt for 0*01 M, 0*015 M and 0*025 M initial concentration of K a S a O 0 
have been extrapolated at 45 minutes (after the end of the induction period). On 
substitution of these values of - do jdt in the expression 

n * lo S (rdcjdt)- 1 og(dc 2 /d t) 
log c x - log c a 

the value 'of n comes out to be 0*97 and'T05, confirming the unimolecularity of the 
reaction with respect to S a O a ". 

EFFECT OF MERCUROUS NITRATE CONCENTRATION 

The reaction was carried out at four different initial concentrations of mer- 
curous nitrate, keeping the concentration of potassium persulphate consta 
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Fig. 3 : Variation of c with time at 50°C 
Cone, of Hg,(NO :i ) 3 : O'Ol M 

Cone, of K 2 S 3 Oa : A 0‘01 Vt ; B 0'015 M ; C 0‘025 M 
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results of these experiments have been plQtted in fig. 4, where the titre values alter 
conversion into mercurous nitrate concentration have been plotted against itime. 
From these curves the time taken for the same fraction (1 /4th) of Hg a (N0 3 ), to 



Fig. 4 : Variation of c with time at 50°G 
Cone, of KoS„ 08 :0*Ql M 

Cone, of Hg 2 '(N0 3 )a : A 0*01 M : B 0*015 M C 0*02 M ; D 0*025 M 
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oxidise comes out to be 53, 88*5, 98 and 130 minutes for 0*01 M, 0*015 M, 0*02 M 
and 0025 M Hg 2 (NO a ) a concentration respectively. On substitution of these values 
in the expression 


« - i + 

log tJd 

the value of n obtained is -0*26, 0*11 and 002, showing thereby that the reaction 
is zero order with respect to Hg.^ + . The value of n by the differential method 
from these curves also comes out to be -0*002, 0T0 and -0*02 confirming the zero 
molecularity with respect to Hg a ++ . Further it is seen that the extent of induction 
period increases with an increase in Hg 2 (N0 3 ) 2 concentration. This suggests that 
probably some thing present in mercurous nitrate solution is responsible for the 
induction period. 

EFFECT OF TEMPERATURE 


In order to determine the temperature coefficient and thereby the energy of 
activation, the reaction was carried out at four different temperature from 40°C. 
to 55°C., the results for which are represented in fig. 5 where log c (in terms of the 



Fig. 5 : Variation of log c with time 

Cone. :,K a S 2 0 8 = Hg 2 (N03) 2 ~ 0*0 i M 
Temperature : A 40°C ; B 45°C ; C 50 °C ; D 55fQ 
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titre values) is plotted against time. From the curves at these four temperatures the 
the value of specific rate has been evaluated and thereby the value of temperature 
coefficient, energy of activation, frequency factor and entropy of activation have 
been calculated. These values have been tabulated in the following table. 

TABLE 2 


K 2 SjO a = Hgj(NO s ) a = 0-01 M ; HNO a = 0-1540 N 


Tempera- 
ture in 
°K 

k unimole- 
cular 
x 10* 

Tempera- 

ture 

coefficient 

1 Energy 

I of 

activation 

Frequency 
factor in 
seer 1 

Entropy of 
activation 
in E. U. 

313 

2-413 

... 

... 

1 X 10 12 

- 3-727 



3*353 

24-32 Cals. 



323 

8-091 

... 

• • • 

1-047 X 10 ia 

-3-698 


318 4-826 

2-984 22-68 Gals. 1-122 x 10 12 -3-530 

328 14-400 ... ... 1-072 X 10 ia -3-683 


Mean values 


3-1685 


23-5 Gals. 


1-06 x 10 12 


- 3-659 


The reaction has a very high temperature coefficient and the curves further 
show that the extent of induction period at higher temperatures is less than that at 
lower temperatures. 

EFFECT OF NITRTC ACID 

In order to decide whether the concentration of HNO s affects the rate, the 
reaction was carried out at three different initial concentrations of HN0 3 , keeping 
the concentration of the reactants the same in each case. It was seen that there is 
no appreciable effect of nitric acid concentration on the rate as is evident from the 
following data. 


TABLE 3 

K. 2 S 2 0 8 = Hgj(N0 3 ) a = 0-01 M ; Temperature = 59°G. 


Concentration of HNO a 

k unimolecular x 10 3 from the curves 


0-1540 N 

7*598 min** 1 


0-4780 N 

7*278 min- 1 


0-8020 N 

7*015 min" 1 
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The specific rate is found to decrease by hardiy 1% for a more than five fold 
increase in the concentration of HNO a . 

CONCLUSION 

Thus it is seen that this uncatalysed oxidation of Hg a ++ with S^Og" is similar 
in behaviour to other uncatalysed reactions of potassium persulphate where the 
reaction is found to be unimolecular with respect to S^Og^ and zero molecular with 
respect to the reductant. 

Further work to elucidate the exact mechanism of the reaction, to determine 
the cause of induction period and the effect of catalytic ions is in progress. 
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ABSTRACT 

Meijer gave two generalisations of the classical Laplace transform in the form 
co 

O) = J-| / JP* Krn(pt)f«)dt (1.1) 

O 

03 

and 4>(p) =p e~^ W k+lm {pt) f (t) dl. (1.2) 

O 

In the present paper some of the properties of the transform (l.l) have been given, which are 
useful in evaluating infinite integrals. One property of the transform represented by (1.2) has also been 
established at the end. 

1. Meijer (1940) gave a generalisation of the well-known Laplace transform 
viz; 9 

03 

£ ip) = P J •’*/(*) dt > R ( P ) > °> (l.l) 

o 

by means of the relation given by 

03 

4>(P) f W* k m Wf(.t)dt, (1.2; 

O 

He again generalised the Laplace transform in (1941) in <he form of another 
transform now called the Meijer transform given by 

co 

4>(p)=P f iP‘)~ k ' 1 • " & W k+i,m W dt > (1.3) 

o 

For k = ± m, fl.l) reduces to the Laplace transform, 

to 

4> (p) —P j' e ~^ l ‘ f (0 dt. 
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Varma has generalised the Laplace transform, in the form 
HP) $ e ~ lpt W Km ( pt)f{t)dt , 


(1.4) 


Since W (*)==x“"'” + ^ e (L4) reduces to (1.1) when k = - m + 

- 

The object of this note is to investigate some of the properties of the 
generalised Laplace transform given in the form (1.2). At the end I have given a 
property of the transform given by (1.4). 

2. Theorem : — If <f> {$) is the generalised Laplace transform ofg(tf) and g(x) is 
the cosine transform of /(#), then 
co 

(*) — f/U) fU) d y> 

o 

where ^ (j>>) = - % 2 2r ~^ T x ( Y + fi” 1 )y 2r < 

^ r= o (2r) ! 

provided that the integral exists. 

Proof : — Since (s) is the generalised Laplace transform of #(*)> it will be 


given by 

CO 

4> (*) = A 2 

f («)■ 

V TT 

J 

o 

Also, g (x) = 

CO 

d~ f 

''TT J 


CD CO 

Therefore, (s) J ( sx) ^ k m (sx) dx J cos (xy)f(y) dy . 


Changing the order of integration, which is permissible, we get 

CD CD 

0) = — vT j f(y)dy jW cos (*y) (**) <fx, 

o b 

- | vT f / 0) * 1 ,»J , 2 '+» M * 


A’oie : f* (a±6) = f ( a + £) r ( a - b) 
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j t • f y * * * , ? „ 

Evaluating the integral, with the help of the known integral, vik. 

J / 1 k m (x)dx = 2^ r # ( Y ± 1>) . ± m)>o. 

O 

2 « (-l) f o-(2f+l) 2 2r-^ r ( r+ l ± Z L^x 

we get <f>W = - ^ Z *\ , + 4* 2 / 

CO w 

x J y* f ( y ) = J $ (y) f U) dy> 

O o 

-«*W-| r l o ^- (2f+5 > 2 2r '* r. ( ,+i ±| )/’. 

Thus the theorem is established. 

The change in the order of integration will be justified if 
oo 

(*) J7(*) *m (*) dx 

o 

Co 

and (ii) j V~ k m (*) cos (*J*) 

o 

are both absolutely convergent. Assuming (1) to be absolutely convergent, (it) will 
converge absolutely if ^ r+^| ) - >0 > 

Example 1. Consider the relation 
J's - * 2 cos (xy) dx = V| e - ^ 4 

O 

CO 

We have, Jn» x (y) fy 

o 

(-1/ (2r+ 3/2) 2 2r-i r / r+ 3 , m \ ^ 2r 

w ., = o ( 2r ) 1 \ r+ 4 £ 2' J 

Evaluating this, we have 

J / W x W * - I l *-^ + : 3/21 r» (r + l ±?> 2 ' - 3 ' 2 *+B 

o 

Example 2. We know, (Erdelyi 307 — 4), that 
os 

J cos (st) t” ~ Erfc (at) dt 
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r c -L+ 

2 2 


v7r a 


rp / V V + 1 . 1 V . — \ 

a MT ’ 'T~’ T ’ 2 + 1; 1^: 


Now 


CO 

J/(0'x(0 


[ £ (»)>o, R (a) >6.) 


dt 


= f 1 + 3 ' 2 > 2 2 '-* r, ( r + ± ±m) f^-‘ Ei!c . M <, 

r ° o 

On using the result (Erdelyi 306-1), 


r (•£+« 


uj 2 

f/- 1 Erfc (ax) dx = — , ( R (P) >o, R (a) >o. ] 

J 7 rp a r 


o 

we get 

CO 


Jf (t) x (i) dt— ^ X s C 2 '' + 3 / 2 ) 2 2r “2 r 0+ J-M) r x (r+ |±^) 

7 T (2r+z>) a 2r-J-& 

3. Theorem 2. If <#> ( 5 ) is the generalised Laplace Transform of g (a;) and g (a:) 
is the sine transfora of /(*), then 

co 

(*) = J/OO X (y) dy, 


9 CO 

where x O) = ~ % 


( _ l)r > y 2r+ 1 


r,2 r + j 


r, :fr+ ■» ± ” ■' 


7 T f=0 ( 2 r+l) l 2r+3/2 

S’ . _ i ■ . v ■ 

The proof is after the manner of the above theorem. 

4. Theorem 3* If <j> ( s ) be the generalised Laplace transform of 2 (x) and 2 (aA 
is the Humbert transform of /(a?), then * } 

co 


> (*) = jf x(y)f(y)dy. 


where 

. (,) = /r 1 ( - , i rr »<T +,+ l ± "> 


^ ,=< 


p+2r+8/2 5 -(o+2r+5/2) 


~o r (r+3/2) r (r-f-v+ 3/2) 

Proof : Since <p (s) is the generalised Laplace transform it will be given by 
.(*)= J- f k m ( sx > S (*) dx, 

V 7 X J 

o 
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00 

3. g (*) = J l xy H v (xy)f(y) dy. 


Therefore 


00 CO 

, <f> (i) = JJL J" (sx) 1 k m (sx) dx j" J xy H 0 [xy)f(y) dy. 


Changing the order of integration, which is permissible, we get 

CO CO 

<*> (?) = v / I s y f(y)dy j\j * h v <*y) k m (sx) dx, 


f 00 (“I)*- (|) y+2r+l 00 

fJy f (y) dy s ^ — 

/ * s' V r = 0 r(r+^) r(r + c+^) 


o+2rH-3/2 


( J *) 


"/£ j .. V ' T *= V . ( y V+2T + t f(y)dy, 

V,r f ° T (r + 3/2) T (r + 0 - 1 - 3/2) s v+2r + 5l<L / 

CO 

= J x (y)f(>) d y> 

o • 

where x (,) - /T ! C , < [mill / + 2f + 3 ' 2 

1 r (r+3/2) r (r+»+3/2) S ! + 2r + 5 ^ ' 

Thus the theorem is established. 

The change* in the order of integration will be justifiable, if the integrals 
co 

(0 / / 

o'* /(*) k m ( sx ) dx 


■W T x < 2 + ' + ?±*) 


f/ + 2r + 4 /(j>) )^ 


00 

<”> fj~. 


H v (*X k m (**) dx 


s/ * 


are absolutely convergent. Assuming (i) to be absolutely convergent, the second 
converges absolutely as term by term integration is permissible. 

Example ). Consider the relation (Erdelyi 16T19) 


x m ~ 3 r aK H v (xy) V *J dx 


v + 3/2 f (i> + m+2) 


2" ,"+'+ 2 y- K.+|) 


» / i m+z> | i m+fl + 3 . 3 ,3 y % ) . 

a^a^l, — 2 — ^ ^ ’ 2’ 
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where k (a) >o, J?(m + o)> -2. 


00 

Now, (fU) XU) dy. 


1° % 


C-I) r r * ( 2 } 


Vw r = o r (r + 3/2) r(r+o+3/2J 


- (o 4- 2f -f- 5/2^ J' 


-aj> m + o+l+2r 


. nr r j. r j. 5 j. e~(»+2r+5/2) r.m+«+2+2r) 

fS-co (~ l ) r r x (2 + r + 5 ± 2) 6 - 

V'Vrio r (r+3/2) r (r+a + 3/2) a »»+»+ 3 /2 + 2r 


5. Theorem 4. If </> (0 be the generalised Laplace transform of / (*), then 


P jkS£L-. ds = 2~P r *(t " | ± 2 } 

J sp TT 


CD 

o Jf 1 ' 


provided that the integrals involved are convergent. 

CO 

Proof We know that </> (*) = J~ f v / i y^ £ 

v v J i 


m (sx) f (x) dx. 


Therefore, 


J ~HT‘ Jv J {sx)k kfn (l * )/( * } 


Changing the order of integration, which is permissible, we get 
J dx — Jj* J x^ f ( x ) dx J k m (sx) dx. 


On evaluating the integral by a well-known formula, we get 


J-±$- Jx 




P f ( ^ P , »2 \ 

J * -4 2 ± 2 ) 


C /(*)<** 

J x'~P 


'fhus the result is established 
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1 he charge in the order of integration will be justifiable if both the integrals 


co 

(0 j /(*) (sxj dx 


O 

CO 


and 


(»') / 


. w 


k m ( <x ) d*, 


converge absolutely. Assuming (i) to be absolutely convergent, (ti) converges 
absolutely if _/>±^~^>o. 


Particular Case : (1) When /> = i , we have 


os 


±(s} 


/(*) 


f^A = ilL / 1± » ) f — & 

JO v/2ir *\2 2/Jv'x 


6. Theorem 5 : — If (s) be the generalised Laplace transform (given by Varma) 

of / (#) then 

provided that the integrals involved exist. 

Proof fallows from theorem 4 given above 
Particular Cases : — 

(j) When ^ = 1 , we get 


? Hil dx = 

T x (m 4- 1 ±m) f 

J S 

0 

r (m — A+ 3/ a ) J 

(ii) When k — 

- m + i > then 

T *w Js- 

J s p 

0 5 

r < 2 -» / f f-p 

0 * 


yrhere <f> ( s ) and / (#) are the Laplace transforms of each other, 
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(in) When k = - rn + J and p~ t , we get 

f T ./M. dx 

J s J X 

o Q 

Example — Let <p (s) = 2 a k a s m "' ^ A:2 m ) 

. r , . — alx k — m - % , 

and /(#)=* ' * 3 


Then, by the above theorem, we have 



© 


iWl+l fe 2m C 2 ^ 5 ) 




r(2~p) r ( 2-^+2»i) 

r' v m -/> - &+ 5 ^a) 


/ 

o 


-ajx m+ k — 3/2 

e» <v 

7^ 


dx 


On evaluating the integral on the left hand side, we get 

< 2 a~% j s m -P +l k 2m ( 2 Vasj~ 'ds = a k +P ~ m - 2 ^ ( m _^ + 2 ±m). 

O 

The right hand side integral 

J ,-«/* x k +P- m ~ V l dx = a k ~ p ~ m ~ 2 1 •*(«+■§-*-/>)• 

O 

Thus the left hand side = right hand side, and the theorem is justified. 

The author wishes to express his respectful tbanks to Dr. Brij Mohan for his 
generous help, suggestion and guidance in the preparation of this paper. 
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NOTE ON THE BOUNDS ON THE NORM OF A SOLUTION OF A 
SYSTEM OF NON-LINEAR DIFFERENTIAL EQUATIONS 


By 

JAGDISH CHANDRA 

Department of Mathematic t, S. V. University, Tirupati 
[Received on 14th May, 1962] 

ABSTRACT 

The author has discussed the asymptotic behaviour and obtained criteria for the uniqueness and 
stability of solutions of a system of differential equations using the notion of maximal and minimal 
solutions of a scalar equation. Continuing this approach, in this note, the author obtains the bounds 
on the norm of a solution of such a system. 

In this note we obtain upper and lower bounds on the norm of a solution 
of a system of non-linear differential equations. The results are essentially 
analogous to [2], however, the method of approach is entirely different and depends 
on a lemma established by the author in [1], 

1. We consider the system 

4 (u) 

where % denotes a ««dimensi mai vector and / (*. x) is a given vector field which 
is defined in the product space A = I X R tt and is continuous there, where I is 
the interval o ^ t ^ co and R n is Euclidean space of rc-dimensions. We now state 
the lemma proved in [1] in a suitably modified form so as to suit our requirements 
in § 2. 

Lemma I : 

Let V (f, x) > o be a scalar function also defined on A* We assume that 
V satisfies a Lipschitz condition in x locally. We define 

v* *) = lim SU P -L-["v «+ Ji r,ei) - V (/, *) 1 

h-> o + A L J 

Let w (f, r) 'be measurable in i for fixed r and continuous in r for fixed t. 
Suppose 

V* (*, x) « (f, V (* f *) ) - « V (i, *) 
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where a > o is a suitably chosen number. 

Further let 
| f{t,x) [ < V(i,*) 

where | x | indicates the Euclidean norm of x, 

Let r (/) be the maximal solution of 

r‘ — m (l, r ) 

with r (t 0 ) = r oJ r 0 > o. Then if x (t) is any solution of (1.1) such that 

V (t, x (t 0 ) r 0 > we have 

V (f, * (t).) sS r (t) for t > t 0 (1,2) 

We shall also require the following lemma, which is essentially an extension 
of iemma I, to prove the Tneorem 2 in § 2. 


Lemma II : 

Let V (/, ^ o be a scalar function defined on . 

we assume that V satisfies the inequality 

V (t, Xj) — V (/, x%) ^ — L | — x 3 j 

where L > o is a constant, Fu iher let 

l /(*,*) 1 < V(t. x) 

As in lemma I, we define V* (t, x) and 0 (/, p ). 

Suppose 

V* (t, *)> - “ [t, V (t, *) ) + 0 V (,<, *) 
where fi > o is a suitably ehosen constant. Then 
whenever o < p 0 < V U 0 x ( t 0 ) ), we have 
P 0) < V (/, x it) ) for O t 0 
where P (/) is the minimal solution of 

^ " “ (<» p ) wit h P (t B ) = P 0 , f'o > o 


(1.3) 

(1.4) 


0 . 5 ) 


( 1 . 6 ) 
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The proof is quite similar to lemtfta t and we shall merely indicate* ]?roif* 
(1.3), (1.4) and (1.5) we observe that 

lim sup 1 r v (t+h, x + h f) - V (<,*. 1 > - - (t, *> 
h — > o + « J 

If* (t) is any solution of (1.1), defining m (0 = V (f, * (€)), we get 

lim. sup. J r w (i + ft) - w (/) 1 > - o) ( /, wi (0 ) 
h o + h L J 


Now suppose that b (t t s) is a solution of 

p' « - « (t, p) 4- s, e > l\ with P (/ Q ) = P 0 > 0 

Proceeding on the similar lines as in lemma I, we obtain the result of the 

lemma. 


2 . 


Theofem : 1 

Let the assumptions of lemma I be satisfied and let x ( t ) be continuous, satisfy 
| x (t) | > o and be a solution of (1,1) for a t ^ b. Then for a ^ t b, we 
have 1 *(<) | *(*) !)+(<-«)] * (2.1) 

u 

where P (ti) * r {s)]~ l ds, u Q > o (2,2) 

«0 


ftw/ : We have = / ( t, x ), so that 
at 


t 

*(<) = * (* e ) + J/ ( i, * (*) ) * 

/o 


so that 


I # (<) | < I *(*o) I + / I /(*»#(*)) l«fr 

*0 
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p GKO) - p (*(*,))< <-*• 

which gives 

^ (<) < P- 1 [ P (/„))+(<- f 0 ) ] where 
P _1 (u) is the inverse of P (u). Since 
t (*<>} = y (to)* we have 

*(0< P'MP (r«e)J + (/ - O J 

and trom (2.3) and (2.4), we get 

,y(<)< p ’ l [P(.y(fo)) + (<- * 0 ) ]. 

Now changing ( 0 to a we get (2.1 ). 


(2.3) 


(2.4) 
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Theorem 2 : 


Let the assumptions of lemma II be satisfied and let x (<) be continuous, 
satisfy | x (t) | > o and be a solution of (1.1) for a b. Then for a < t 

j££ b , we have 

| X(t) I Q; 1 [ I * ( a ) I ) + c 1 - *o ) 1 ( 2 * 5 ) 

where 


u 

& (a) = - f [ P (*) ] ' x ds 

«o 


( 2 . 6 ) 


Proof : 

We have 


I *(0 I > I * (<o) I -/ v (»,*(»))* 

*0 

Setting j> (i) = I #(0 I .we get 
(0 > y (<o) ~ / v ( * w ) 


Let 


(2.7) 


t 

4 >{ i)=y (< 0 ) - J v (*. * (0 ) 

then (0 — “V (*i * (0 ) 

From (1.6) and (2.6), we have 

j Q(<f> (t)) > | > integration of which between ( < 0 * O 8' vcs 

dt 

& ( 4 (0) “ & ( 4 (*o) ) > ' " *o> which we 


( 2 . 8 ) 
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* (o > a - 1 [ a ( * v 0 > ) + (t ~ 1 0 ) ] 

Since <f> ( t 0 ) = y (f 0 \ we have 

> qtm aoM) + o - *« > i 

From (2.7) and (2.8) we get 

j>( 0 > a-Mao-co) + (/-<„)]■ 

Again changing / 0 to a we get (2.5). 
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UNSTEADY FLOW OF A VISCOUS FLUID THROUGH A 
RECTANGULAR PIPE 


By 

G. TEEKARAO 

Department of Chemical Technology , 0>mania University , Hyderabad ( A . P.) 

[ Received on 25th January, 1962 ] 

ABSTRACT 

In the present note the technique of super posablity is adopted to solve the Navier- Stokes equation. 
Two solutions of the flow of a viscous liquid through a rectangular pipe are jobatined. under different 
pressure gradients. The superposed solution is so adjusted that the velocity is zero for t~o. The 
particular case of the flow between two parallel planes is also discussed. 

Recently Ram Ballabh has studied the unsteady flow of a viscous fluid 
through a circular pipe (l). He has used the method of superposability to obtain 
the solution which vanishes initially and also satifies the viscosity condition on 
the walls of the tube. In the present note we extend his method for a pipe of 
rectangular cross-section. 

The axis of the cylinder is taken along the Z-axis and the velocity in that 
direction is denoted by U * From the equation of continuity it follows that u 
depends only on the x,y coordinates and the time t . We first obtain two particular 
solutions. These two solutions are superposed and the arbitrary constants are so 
adjusted that the velocity vanishes for t=o. 

The boundary of any cross-section is taken as 

-**)(/ - *•) = o (oa) 

The velocity vanishes for * = + a,y = =F b. 

1. The equations of motion are 


vV** - 


du_ 

dt 


l di 

P dz 


(i.i) 


dp _ dp = Q 

a* dy 


( 1 . 2 ) 


From 1.2 we see that the fluid pressure is constant on any cross-section of the 
cylinder perpendicular to the axis. We consider two cases. 


(1) - 

dp 

dz 

- = 0 , (II) — — — = a function of / 

dz 

Case (I) - 

dp . 
02 

= 0 
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The equation 1.1 is solved assuming that the velocity contains an exponential 
decreasing time factor. We therefore take 

U = e~ vAH f(xj) (1.3) 

The function f (x,y) vanishes for x^^fay^^b. 

Substituting from 1.3 in 1.1 we have 

(V a + A*)f{ XJ ) = 0 (1.4) 

with 

/(¥ a,y) = o, f (x, qp b) = o (1.5) 


The regular solution of 1.4 satisfying 1.5 is 

/(*.j) = C»(« cos X cos ^ 1 ) Tty (1.6) 

and 


,44 u _ (2m+l)<M (2 n+\)W 

xx m,n — r~5 T - 


4a a 




The velocity is given by 


u = e cos 


( 


5 -5L+l )„ x co,( 


2a 


The equation 1.1 being linear the solutions are super posable therefore u can 
be taken as 

a = S S e ~ vA * m ’ n 1 cos ( 1 ) * x cos ( (1.7) 

where the constants will have to be so chosen that the double series is uniformly 
convergent and can be differentiated term by term twice. The summation* in 1.7 
ranges over all positive integral values of m and n. 

Case (II). 

We take 

_!!_ = BPe~ at (1.8) 

u now satisfies” the equation 

v-V 2 « - -1“ = Be~ at (1.9) 

d£ 

Let 


u 


e at F (x,y) 


( m 1 



( 1 . 10 ) 


SubsiitutingTn 1.9 we get 


( 





O 


The 1 unction F(x,y) satisfies the same boundary conditions as the function 
/(#*?)» given by the equation 1.5. 

Expressed as a Fourier series valid over the range | y j <6 the constant 
in 1.10 is 


B- *L I <-■>* 

^ 77 = 1 2/1+ A 




( 


( 1 . 11 ) 


We shall assume that the function F (x t y) is of the form 




■ 1656 2 

TT 3 V 


(-1)” 


COS 




(2n+l)3 \ 2b 

1 1 + m ] 

where S(x) satisfies the differential equation 




(1.12) 


S"(x) - ( JV 1 - JL. ) S(x) + + = o 


(1.13) 


with 


4 6 a ft* = (2n+l)M 


S(x) has to satisfy the boundary cinditions (in view of the boundary conditions 
satisfied by F (x,y)]. 

1 + S(a) =5 0 , 1 + S(-fl) = o (1.14) 

the solution of 1.1 3 -is 


S(x) =' A x e^ 1 W - «M * + B l e ~ “ a l v ) x 


+ 


a.- 

v[JSff — a/v) 


(US) 


Determining the arbitrary constants A x and B ± from LI 4 it is easy to obtain 


S x *)~ - 


N l 


cosh \'{JV%-alv;x 


+ - 


- a 

~ CL 


JN* - ajv 


cos h ^(N‘ i -~a/v)a 

[ -2*7 , ] 


(1.16) 



Thus the function P{x,y) is given by 


F(x,y) = 


-4B (-1)" 


w 2(n+i) 


cos 


( l _ COS ft p X \ 

^ c s h P a ) 


c 


2n+l 

2b 


7ty 


(1.17) 


where 


p* = N* 2L 

V 


The equation of motion being linear the two solutions are additive. The 
superposed solution is 


“ = ssc„r v ^‘ 2m+1 


cos 


( 




2a 






wv (2n+l) i 8* 

2. Evaluation of the Constants. 


COS 


( ^ )* (I- S4j*> (■•*« 


We shall evoluate the coefficients C ra , n by assuming the velocity to vanish 
at 2=0. For 2 = 0 we get from 1.18 


f o, j •“+'? r s 


,= 1 


2# 


[ S C„,„ co, C ^5+L ) 

L m—\ \ 2 a ' 


7 XX 


- jjLLz 1)” / ^ cos h fi x \ i 
7rv(2n+])p\ ~ cos h pa~~ /J 


(2.1) 


If this is to vanish we must have 


CO 

2 Cnun cos 
m=l 


( - 2 ™+ l W = jj i ( ~ ? )W ( 1- cos h P* \ (2.2) 

2^ / 7t v (2n-\-\)(3 \ cos h pa ) ^ 

The right-hand side vanishes for x = =F a. The constants C m>n are given by 


m r a 


„r _ «5C-n . , , 

'"""SWP) 1 < 1 


/< 


Cf)S * \ cos J2m ' 1) irxdx 
3a / 


cos A pc 


2a 
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M 


Evaluating the integral we have 

8/?(- i) m +” / 1 

aGm,n = ttv (2n+iyW > ^ 

where 


) 


2 . 3 ) 


2a M = (2m+l) * 
The velocity is given by 


« =— ss 

TTV 


COS 


7TV 


(-1)'”+” r v ( Mt+ W‘ ( \ m \ 

(2n+l)flj8» V7T’ Af a +/? a / 


4B«~“* (-1)” 

2 -k , cos 


(2/1+1)/?* 


If a =0, we have 
B 


u = 


irv 


■ S2 


(■ 


■ D OT+ ” ,-^+jya) / ^ ^ 


COS 


( 


(2«+l)jV s a \M M a +j8 2 

2 ?7) + 1 \ ✓ 2n+l \ 

)^cos( — £-)v ~ 


) 


2a 


1606* ^ (-1) 
7r»v S (2n4-l)» 


coj i 


( ^±j~) 2M) ( l_ COS A ^ 

\ 26 ^ V cos A Ba f 


26 / V cos h fia 

As f co this reduces to the well-known classical solution (2) 

V » /4 [ A* - ^ (— l) n "^* jV~ 3 Sech JVa cos A JV# cos Ny J 

where A is a constant. 

3. Flow Between Parallel Planes. 

As a limiting case we consider the flow between two fixed parallel planes 
y=zb 3 y~ — b> Case (I) gives 


9 *u 

dy* 


dti 

IT 


= 0 


(3.1> 
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U vanishes for y=*b, j>= - b. It is easy to verify that u is given by 


« - S C, cos ( ?!±i ) ,y (3.2) 

Case (II), taking 


u = 


— at 


F(y), 


dp _ 

dt 


PB 


- -at 


the equation of motion gives 


dy 

which is satisfied if 


W ^ — {aF — B) = o 


K B , Bi j a 

t = -h — cos ; ^ 

a « V v 7 


The boundary conditions give 


O = 5 cos ~ a b 
V ■> 


so that 


B 


Bi 


cos 


/V 


V' 


Thus the velocity is given by 


/“ 
COS w 

5 -at t , V v 

U as £ 1 1 - ~ 


a 


( 


cos 


/ or 

V“ 


) 


The superposed solution is 
u — % C n e 


— vN%t / 2n-M 


COS I 


26 


77^ 


(3.3) 


(3.4) 
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, B -at 

+ e 


(- 


cos , / a 

v" v y 

cos V“r * 


) 


To evaluate the constants we impose the condition that w is 
For this we take 


cos 


1- 


l~^r 


/— T- 

cos /_? L t 


A n cos ( ) wji 




where 


= 2(-l) n ( * a ), 4 v 7 -?- 4 


^2 


At *=0 3*5 gives 

u = § £ C n 4* 


2B(-l)» / 1 __ JV 

^ JV ^ 


6a 


jV*- 


>] 




If this is to vanish we choose 


G n = 28 (- 1 ^( -JV'-V) 

6a ' xra & N * 


JV*- 


As a -> 0 it is easy to show that 

r - 2g(-i) n 
« 6v TV 3 

For any other value of f the velocity is given by 


“x <-*> 

ba 


N 


1 \ - vJV 1 / 


JV 2 - 


** JV 


) 


cos 


6* 


+ ^ 1 - cos 

a \ 


/-r^ &c v-T-* ) 


(3.5) 

zero for <=0. 


(3.6) 
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At a = 0 this gives, from 3.6 and 3*3 


" - -2- * » ( ^-) * - fr <’*- **) 

As t t~> co this gives the classical solution ( 3 ), 
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A NOTE ON FUNCTIONALLY BOUNDEDNESS OF SOLUTIONS 
OF NON-LINEAR DIFFERENTIAL EQUATIONS 


By 

M, RAMA MOHANA RAO 

Department of Mathematics , Osmania University , Hyderabad-7 , A, P . 

[Received on 24th July, 1962] 

ABSTRACT 

In this note we shall derive some criteria, for the boundedness of solutions of# f(t t x) as co i n 
a sufficiently general way so a9 to include the results of V. Lakshmikantham, B. Viswanath am and 
A. Wintner as special cases. 

1. Let x and/ (x, t) be vectors with real components (x x it). x 2 (t) ,.x a (t)) 
and — /a) respectively. 

Consider the differential system 

(1.1) x 1 = / (#, t) x(t Q ) = 0 (' = djdt) where J (x, t) be continuous in the 
region t 0 < t < oo > j x ] < co 

In this note, we shall derive some criteria for the boundedness of the solutions 
of (1.1) as t ->co , in a . sufficiently general way so as to include the results of 
[ 3 ], [ 4 ], [ 5 ], as special cases. 

We shall suppose that the function <£ (x it) ,t) and another function o» (r (*)> 0 
respectively possess the following properties. 

(A) (; x , t) is a function defined for vector x and real t with nonnegative 
real values, which is continuous in (x,t) has one sided partial derivatives with respect 
to t and the components of * and <j> (*, 0 = 0 implies x = 0. Let <£t , </>* > 0 
be partial derivatives of <j> with respect to t and x respectively. 

(B) (o (r, t) is continuous, non-negative function defined for r ^ 0, t 0 ^ /<co 
and x (0 the maximal solution of the differetial epuation r' = at (r, t) 
passing through the point (/ o ,0). 

The functional boundedness of solutions of differential equation (1.1) is 
defined as follows: 

A solutions # (t) of the differential equation (1.1) is said to be functionally 
bounded with respect to <f> and L , if it satisfies < 03 whcre the fac- 

tion L (i) > O be continous in the region t 0 sg t < os . 



We shall first* prove the following: 


Lemma — Assume that (A) and (B) hold and let the function / (#> 0 of the 
differential system (1. 1 ) satisfy the condition (1.2) </>t (#, 0 + <£x (#>*)• 1 / (*> 0 I 

^ co [ £ (#,*)>< J tben ^ (*> 0 < X (0 
Proof: — Let m (t) * </>. (*,/) 

Taking right hand derivatives with respect to t , we get 

| m' (t) | 4 > t (*, /) + <j> x (*, t). | | 

di 

I »»' (0 I < 4>i (*. 0 + (*> 0 I / (*» 0 I 

By condition (1.2) | m‘ (t) | < a £ 4> [x, t), t J 

^ O ^ ffl (i) , ( J 


Suppose x (t> E ) is the solution of r' — <o (r, t) + s r(t 0 ) = 0 where s is 
a positive quantity. The maximal solution of r f « w (r, /) + e is given by 


lim 


x(^«)-x(0 [2] 


We shall first show that 

m (0 < X (*> «) < 1 

Suppose at a point t > 0, the inequality is not satisfied. Let [ t 09 / a ] be the 
interval where 

m (f) > x (b *) 

At we have m (f 0 ) = x (*o> e ) 

Hence by taking right hand derivatives at t ot we obtain 
m OoO > X' ( f o» £ ) 


From this obtain 

Cl) (^o) ] ^ w [ ^05 X (/o> e )] e 

This obviously a contradiction. Hence m (/) ^ x (*> e ) 



Making s 0, we get 

m OX X (0 

Note: - If we put $ (x, t) = | x j , we obtain the results of [ 3 J and [ 4 ] 
which are the generalisations of Bellman’s lemma c. f. [ 1 ], 

Theorem 1: — Assume that (A) and (B) hold and let the function /(#,<) of the 
differential system (1.1) satisfy t he condition 

(1.3) *< ( | * ! , o +**(1*1,0- I / (*, 0 I <•[*(*. 0. 0 and if x (0 - 

0(L(£)) as t -» co , then every solution of (1.1) is functionally bounded as co. In 
particular if x (0 = 0 (L {t) ) as co t then each component of every solution of 
(1. 1) tend to zero as co . 

Proof i — Let m (£) = * ( | x | , t) 

Taking right hand derivatives with respect to t , we get 

I m* (0 1 ^ *t ( I x | , t) + <j> x ( | x | , 0- I /(* 9 0 I 

(i.e.) | m' (f) | <o[*(|*|»0»O 

co [ m (f )» 1 j 

By applying the lemma, we get 

(0 

(i«.) *(|*|,0<X(0 

As t & t and with the assumptions of the theorem follows the desired 
results. 

Motci — The results of [ 3 ] are obtained by putting *(*,*) = | x | . 

Theorem 2: — Assume that (A) and (B) hold and let of the differential 

system (1.1) satisfy the condition 

(1.4) fa ( | Xi ~ | , 0 + 1 > 0 * I / 0 * f( x v 0 I 

< « L * ( I *i - #i I , 0. *] 

where x x and x 2 are two solution of (1.1) suppose x (0 8=8 0 (L (0 ) as f~»co then, 
if one solution of (1.1) is functionally bounded, then every solution vector is 
functionally bounded as £-=► . In particular if x — 0 (L (0 ) as co, then 

every solution tends to the same functional limit as t -> co . 
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Proof * — Since x L and x, A be the solutions of (1.1) 

Writing | v (t) | = | x x (t -* a (t) | 

I *' (0 I = I *1 (0 - V (01 =' I /(*! (0. 0 - /(*a (0. 0 I 
Let »t (0 — ( I 0 i > 0 

Taking right hand derivatives with respect to t 

I «' (0 1 < ( I 0 I . 0 + 4>x (1 V | , t). | / (x v t) - / (* 8 , f) | 

< “ [ * ( M . 0, 0 

o> [ m (0, t] 

By applying lemm#, we get 

» (0 < x (0 

i.e. <j> ( | » | , t) < x (0 

As £ — >co , and with the assumptions of the theorem yields the desired results. 

If we put <j> (*, 0 *= | x | , we obtain the results of [ 3 ], 

In a similar way, we can derive the theorems 3, 4, and 5 of [ 3 ] in sufficiently 
general way. 

In the end, I wish to thank Dr. B. Viswanatham for his criticism during 
the preperation of this note. 
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IfttERFERENCE OF KETONES IN THE ESTIMATION OF PERFUMERY 
ALCOHOLS IN SYNTHETIC MIXTURES BY ACETYLATION. PART III. 


By 


JAGRAJ BEHARI LAL 

Chemical Engineering and Tech . Seciion } H \ B. Technological Institute 5 Kanpur 
[Received on 25th September, 1962] 


ABSTRACT 

In the classical method of estimation of primary and secondary alcohols by acetylation with 
acetic anhydride and fused sodium acetate, the presence of acetophenone causes interference and gives 
higher results for esters of alcohol and lower results for acetophenone in the acetylated product. Pre- 
viously this interference was reported to be due to the ketonic group of acetophenone with a methylene 
group of the perfumery alcohol. In continuation of previous work considering that perfumery alcohol 
and the enolic form of ketone undergo acetylation, formulae have been derived for the calculation of 
percentage of individual alcohols persent in a pentenary mixture consisting of two primary or 
secondary alcohols, an acid, one ester and one ketone. 

Consider a mixture containing a , b , c t d and k percent of alcohols A, B, ester 
C, acid D and interfering ketone K respectively . It is assumed that their 
respective molecular weights are hVMi»,M c ,Md, and M k , the ester value of 100 
per cent pure acetate of alcohol A and B and ester C being V a ,V b and V c respec- 
tively. Let Vd be the acid value of 100 percent acid D and v x and o 2 respectively 
represent the ester value of the mixture before and after acetylation and v z and z/ 4 
be the acid value of the mixture before and after acetylation. 


The mixture after acetylation with hot acetic anhydride and fused sodium 
acetate contains k l per cent of interfering ketone* On the basis of original mixture 
let k ± and (fc — k x ) percent be unchanged ketone (ketonic form) and enolic form of 
the interfering ketone which gives ester having ester value Vic. Then, before acetyla- 
tion a -f* b d" 4 c -f- d k — 100 ... ... (1) 


IOOzt , v\ Me 
V c + 561-04 


( 2 ) 


100y 3 v z M& 
561*04 




k 


100 xt 

a* 


(4) 


Where x x and x^ are the respective oxime number of the mixture and pure ketono 
K. Ther, after acetylation with hot acetic anhydride and s.dium acetate, the 
proportioi of acetates of A and B, ester G, anc acid D and aretate of enolic form 
of interfering ketone and ketone K in the ketone form would be as follows : 


A ( 1 + 


42-01 

M* 


>:* ( 1 + 


42*01 

Mb 
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'then, after hot acetylation 


V 2 = 


“ (1 + nar* v * + b (' + ) f> + cn +'(*-*,> <i + Jjj**.) n 


‘ (1 + + S(1 + THC-> + « + ^ + ( t- 1.) (l + i^l) + *, 


(5) 


4201 
fb 


M k 


As (Ma + 42*01) Fa = (M b + 42*01) F b = (M k + 42*01) Fk — 56104 
and 6 = 100 -a — c-d-k 


F„ = 


56104 a 56104 (100 -a-c-d-k) . 56104 (k-k,) . 1ftn 

M. + mT + + 10 ° Bl 


100 + 42-01 a + 42*01 (100 - g - c - d - fc) .. . . 42*01 

Ma Mb + (& “ Al) ~m7 


(5a) 


100t> a 


42*01«y a , 42*01 ( 100 - a ~ c - d - k) v, 

I -i A • -a r ~ 


Ma 


56104a . 56104 


Mb 


Wt .+ 2 
M k 


Ma 


+ ^gL(100 - * - c - rf -j) + 561 04^) + 100|>j 


Mk 


561.04a _ 56104 a _ 42'01ao 2 42*0la» 

° r Ma "Mb Ml + Mb ' 

- 51!^ + ioo,., + «JI 


(100 - C - d - A) r a - 100oi + (A - A,) — 42 ' 0 1 Pg 


M k 


or, 56104a [ 


1 1 


M. Af b ] - 42 ' 0I "'[ ‘Hr ” THT ] “ ,00,J ' ‘“"a - 

W* (56104 - 42 ' 01 '=>- ( 1 °° - ‘u, - ~ *> -(56104 - 42 01,,) 

“• “ (lb ~ A?,] < 56104 " 42 ' 01 '>> “ l00 ’i * X 50 '. ~ (56104 - 42-01,,) 

r 10 0 - < - a - * i (* - *i)i 

*■ M b M k J (*) 
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( 7 ) 


or a = 


M a 


M a - Mb 


£ (100 -- c - d ~ k) 4 


(/c - fci) Mb 


Mk 


(e a - gj) Mb 1 

561-04“ 0-420 1» 2 J 

Next, substituting a = (100 - b - c - d - k)va equation (5) 


56104 


+ ----- + CF 0 + (k - k t ) (1 + 4 -H£L) Fk 

Ma Mb Mu 


, 42-0, (100 - b - c - d -it) , 42-015 ^ 42-01 ( k - fc x ) 

ioo + — — + -ir + Mk 

, 42-01 (100 - b - c - d -ki v % 42-01 bv., , 42-01 (5 - &,)b» 

or ioo„ + i 3JT MT" + M," - 


56104 (100 - b- c- d-k) 561045 , n . n , 56J.04 (5 - fcj) 

= M. iW b 1 Mk 


or 


^2 - *l) 


561-04 - 0-42018, 


(100 — b ~ c - d 

M. 


M, 5 (A ; - *, ) 

Mb M ft 


«““• » “ irr^t < 100 - ‘ 


d - k) + 


■k - *b)M, 

Ms 


- - gi ) M a 1 

561 04 - 0-4208 2 J 


( 8 ) 


With the help of the equation (7) and (8) it is possible to determine the 
percentage of individual perfumery alcohol in the mixture provided the value of 
k x and ( k — fc x ) which are respectively the percentage of total ketone and the 
enolic form of ketone in the original mixture. The percent ot ketone, in the 
acetylated product 


kf = 


100ij 


( 9 ) 


a (1 + 


42-01 

M* 


) + b (1 + 


42-01 

Ms 


■) + £ + d 4- {k — k^ (1 + 


42-01 

M k 


■) + 


Since a J fb+c + d+ k — 100 

lOWq 

k' = 

, nn , 42-0 la , 42-015 , 42'01 (fc - k x ) 

iOVJ -j- , — — i- — Z — T ' 

M a Mb Mk 


( 10 ) 
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Hence, 


]0Qk , + 42-OlaA' + 42'OIM/ + 42-OUifc' 


M a 


M fc 


Ml 


42-OlAiA' 

Mk 


= 100*! 


or 100*, +_ ffiS W j = |00jt' + 42 ' 0M ' + 42 *0*^ + *2-01W 

M k — y> M b M lc 


Ma 


or 100*, [ 1 + = 42-OH' [ ~+ -A- + JL ] + 

L Mk J L M a Mb Mk J 


100*' 


M k J L Ma ‘ M b 1 Mk 

-‘'[ 42 - 01 ( M T + i + A) +H, °] 


Hence A* 


** r 1 + 0-4201 I JL. + 

*• l M. 


(100 - a-c-d-k) , * 


M b 


M k 


}] 

-(11) 


r 1 + 0-4201*' -j 

L Mk J 


The^ equation () 1) enables one to calculate k t the percentage of ketone in 
the ketonic form in the original mixture provided a is eliminated from it. 

The acid value of the acetylated mixture v & 


dv A 


a ( 1 + igi) + 4 ( , + 4 -g‘) + , + „ + (t - U + igj +il 


dv a 


a + + 6 + 42 ' 01 * 


M« 


, , j . , , 42-01A: , 42*01*! , t 

Mb M k 1 Mk 


Since a + * + c + rf+*=100 



dv d 


100 + 42-01 I JL. + _A + (* - h ) 

1 M a Mb Mk 



100 + 42-01 


a 100 — a — c - d —k . (k — k ± )] 

MT + M7 + ~ Mk t 
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42-0lat> 4 

Mb 


4 i*01 ( 100 - c - d - k > 4 
+ Mb 


or dvi 


I00a 4 + 


42-0 laa 4 , 


42-01 ( * - Aj)o 4 
M y 


1 1 1 _ dvi _ (100 - c — d - fc) _ (* - * 4 ) _ 100 

AIT ” ) = 42*0 1 v i -Mb My 42’01 


Now equation (10) can be rearranged as follows : 


*1 



1 + 04201 




a \ (100 — e — d — k) * | "I 

/ Mb Af k j J 


0-4201*' 

My 


Substituting by the value of ( ~ M ) * n a ^ ove 


(100 - c - d - k) 

Mb 


1 + 0 - 420 . { 1 ^- - 

, (100 - c - d - k) ^ k \l 

+ Ml + My j J 


(* - * 4 ) - 100 


My 4201 


1 + 


0-421*' 

My 


*' [ 1 + 0-4201 


dv a 


42-0 lr, 


100» 4 ^ *i 


My 


] 


1 + 


JO-4201*' 

My 


v [ 1 + °' 4M1 { + -Jfc-1 ] 


1 + 


0-201*^ 

My 


( 12 ) 


(13) 


! 201 ) 



The equation (13) would enable one to calculate k x the percentage of ketone 
in the ketonic form in the original mixture. 


l, T i _l °’420 \dv & 

k L 1 + 42 - 0I v i 

k x = 

, , 0'4201/t' 

1 + n 


1 + 


0-420 U, 
M k 


] 


P k' dv a 

, 0-420 \k t V -j 

L 100i> 4 

‘ M k J 


or *i “ - ... ... (14) 

, , 0-4201 k’ 

1 + Mk 


or k x 


k* dvj 
100 ^ 


(15) 


The equation (15) enables one to calculate k x from d the percentage of acid 
D having acid value V& and k f percentage of ketone in the acetylated mixture 
having acid value v ±9 Thus, knowing the percentage of enolic form of ketone 
in the mixture it is possible to calculate the percentage of the two individual 
alcohols A and £ present on a penteznary system consisting of two alcohols, one 
acid, one ester and one interfering type of ketone. 

Experimental work on the verification of these formulas is in progress. 
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up£er and lower bounds of the norm of solutions of 

NON LINEAR VOLTERRA INTEGRAL EQUATIONS 

By 
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ABSTRACT 

Consider the system, of iaiegial equations of the form 

u (t) = f 0) + J l A 0 — s) / (r, U (s)) ds 
o 

The existence theorems of such equations hive been considered by J. A. Nohel [ 1962]. Assume 
that the solution of the above equatioa exists, we consider in this note the bounds of the solution in a 
sufficiently general way. 

1. Consider the system of integral equations of the form 

(1.1) a (0 = 9(0+ f A(t -s)f (s,u(s)-) ds 

O 

under the assumptions 

(i) (p and / are given vectors with n components each; A is a given nXn 

matrix defined on o co , |«| < co . 

(ii) r 0) , A ( t ) are continuous functions in o < * < co and f (t , u ) 
is continuous in ( t } u) for all u and for all t. Define II » II — S 1*1 

The existence theorms of such equations have been considered in [ 5 ]• 
Assume that the solution of (1.1) exist, we consider in this note the bounds of the 
solution in a sufficiently general way. 

We require the following 1 lemma before we proceed further. 

Lemma\-~ Tf 

r l 

(1.2) z (0 < k (0 + * j © 0, z (s)) ds 

b 

where a > o , some constant ; « (£, t) > o is continuous in both variables in a 
certain region R: o i < oa , £ > o ; and monotonic increasing in l foi e\tery 
fixed value of V ; z ( t ) and k (i) are continuous functions on the same -range, then 

« (0 < * (0 + M (0 
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where M (t) is the maximal solution of 

(1.3) i' (0 - «<■> 0, 2(0 + *(0) 

through the point (o. o) 

Pros/; — Put r (0 = 2 (0 “ * (0 

then the ineqality (1.2) becomes 

r (t) < a j to (i, r (0 + k (0 ) * 

O 

Take r (f) as the zero approximation to the solution of the differential equation 
(1.3) through (o, o) and set up the successive approximations recursively by 

J 

fk+i (0 = « j «> (s, r k (s) + k (s) ) ds 

6 

we ean show by induction that these successive approximations form a 
monotonic increasing sequence, for suppose that r * > fk-j. 

Then 


f k+l (0 - (0 = * / { <*> fa, ric (^) + k fa) ) 

© 

“ <0 r*.j (0 + & (*)) X ^ 

> o Since « (t, z) is monotonic increasing in Z 

Therefore fu +1 (0 > tk (0 

But the zero approximations ^ first approximation. So, the successive 
approximations form a monotonic increasing function sequence, and since they 
are uniformly bounded and equi-continuous must converge uniformly to a solution 
1 // (/) . Further it is clear that d (t) is a solution of the differential equation 
(1.3) through the point (o, o) and r (t) < <p (<)• 

Hence r (t) < M ( t ) where M ( t ) is the maximal solution of (1,3). This proves 
the lemma. 
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tt. Now we are in a position to prove the following result: 

Theorem i: — Let 

(i) the vector function f of the system (i.l) satisfy 

11/0, U (t) ) (1 < (i) 0, (| u ( t ) (I ) 

where e> satisfies the conditions of the lemma 

(ii) || A (/) 1) < ct for all values of t and a is the same as before and let u (£) 
satisfy || u H > o and be a solution of (1,1) in the interval o ^ t < co. Then for all 
t in o ^ t < co , we have 

(2.1) H a (0 U ^ 11 9 (t) 11 + M (0 

and 

(2.2) || u (0 II > l! 9 ( t ) || + m (t; 

where M (t) and m (t) are the maximal and minimal solutions of 
2'(l) = ±«« (t, Z ( t ) + 9 (0 ) 
through the point (o, o), respectively. 

Proof. — Let u (1) be a solution of (1 .1) 

J 

u (0 = 9 (0 + J A (i- s) f(s , u (t) ) ds 

t 

« “(0 II ^ II 9 in || + f II A(t-s) II II /(*,«(*)) II ds 

© 

Thus is view of conditions (i) and (ii) of the theorem, we have 

J 

II « (0 I < I f « B + « / « (*, II a (i) II ) * 

© 

By applying lemma, we obtain 

II « (0 II « II 9 (0 || + Af (j) 

To prove (2 .2), we have to use essentially the same argument as in lemma, but now 
we have to consider the minimal solution of z' (t) ss — « o (£ z (t) -f~ p (() ) 
through (o, o) instead of maximal solution of z! (/) = a os ^ (/) _j_ ^ ^) ) 
through (o, o). This completes the proof. 
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kemark The above theorem includes some of the results of [ 5 ]. Taking 
9 (t) ~ k x , A (t) = k 2 for some constants k v k z > o we obtain the results of 
V. Lakshmikanth [ 4 ] which includes the results of Bellman [ 1 ], Behari [ 2 ] and 
Langenhop [ 3 ] as special cases. 

It should be noted that the results in this note require that the function <o (t, z) 
satisfy a monotone condition with respect to z. In the coresponding results in 
Ordinary Differential Equations as I have considered in [ 6 ], thus rcquirmcnt is 
not necessary. It appears unlikely that this condition can be dropped in the case 
ol integral equations. Except for this restriction our results contain those mentioned 
above lor Ordinary Differential Equations as very special cases. 
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SOME THEOREMS OF INTEGRAL TRANSFORMS 


C. B. L. VERMA 

Maharaja College, Chhatarpur, (M. P.) 
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ABSTRACT 

In some of his previous papers the author has studied inter-relationships between various integral 
transforms. Continuing the study in this paper also, he has established the relationships of the Laplace 
Transform and its generalizations as given by Meijer and Varma , with the H- and the Y- Transforms. 
Four theorems in this connection have been obtained and they have been illustrated by suitable 
examples. 

1. Introduction : The present paper is in continuation of some previous 
ones ( (8) & (9) ) by the author in which inter-relationships between various integral 
transforms have been studied. 

We call (6) 

g, ( p ) = j CP*) 1 (px)f(x) dx (1‘1) 

O 

the T— Transform of order v of/(*) and regard p as a positive real variable. 

Its reciprocal is called the H — Transform (6) defined by 


■ i 

0) = j (pxr H v ( px)Ax ) 


We shall denote (Tl) and (1*2) symbolically as 


Slip) JLf(x) 
v 

respectively. 

The classical Laplace Transform viz*, 


<t> 0 ) = P j • 


dx, R(p) > O 
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has been generalized by Varma (7) in the form 
00 

*1 iP) = P JV ~ 4 pX (f) M " 1 W k m (px)f(x) dx R(p)> o (1-4) 

O 

while its generalizations by Maijer ( (5) & (4) ) were given as 

oo 

(i) </> e (p) = . p J (*#)* K y (px)f(x) dx (1-5) 

O 

CO 

& (») 4> 3 (P) = P fe~ h pX (px) " M W k+hm (px)f(x) dx, R(p)> o (1-6) 
© 

We shall represent (1*3) to (1 ’6) as 

<P ip) fix) 

4>i(p) ==z fix) 

k,m 

<p* ip) JL fix) 

V 

<Pa ip) ^ fix) 

k t m 

respectively. 

It may be noted that (1‘4), (1‘5) and (1'6) reduce to (1*1) when (i) k — m-\-i 

( ti ) v = ± J and (Hi) k = ± rn respectively. 

We now establish the following theorems : 

2. Theorem 1. If / “ 3 / 2 fix) and 4> ixj belong to L (o, co) 
and if <p (p) ===== x 1 * * 3 ^ 2 /(*) 

*(/>)=£; /(*) 

V 

then 

*(#) = - (^ 2 + /> a r 4 ** a;:, [ ] «#> o» ^ 


where 


R(p)>o, R (/*) > I ft (v) I , -K«(v) 
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(2-2) 


Proof : Titchmarsh (1937) gave the inversion formula for (l*2j in the form 
oo 

a /(*) - j wfc ? v i x y) <t> (>) dy, R(y) > 0 , - $ < R{v) < $ 

O 

Substituting this value of j(x ) in 

P iP) = P j e ~ P** 1 ~ 3 ^ 2 fix) d X 

O 

and interchanging the order of integration which is permissible by Fubini’s theo- 
rem under conditions stated, we obtain 

CD CO 

Pip) = P f <f>(y)[ fe~f X x /J ‘~ 3/2 (*)0 3 r v (xy) dx ] dy (2'3) 

O O 

Evaluating the inner integral by (2), p. 105, we immediately get the theorem. 
3. Theorem 2 \ If /( V"3T ) an d </> (a) belong to L (o, oo) 
and if 


tip) 


4> iP ) 


M 

k, m 

H 


/(V * ) 


fix) 


then 

^ ip) = (- o 




\k - m, k + m. 


4 p 


1 1 


1,1 1 1 0 , 1 1 

A — H -o v » ~A ' V, 2k, - r - — 


) 


R (5/2 — 2& ± 2/72 ± v) > o, - 1/2 < (v) < 1/2, « (£) > o, R (y) > 

Proof : Substituting the value of/ (vGT) from (2.2) in 
co 

Pip)=P f {px)- k -i W k m (px)f(\J~)dx 

o 2> 

and interchanging the order of integration and using [ (2) p. 117] 

00 

jx 2X (xyy- e~ X ^ X% W k fi (| * a ) 7 (xy) dx 

O 

- z 1 - 


P (y)dy 
(3.1) 

O 


= (-l) m 2 


/t, k 1 3 k — \ i 
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h - — L 4 \ v, 2T' = -I - 4 V, l = - 1-. - I V 
4 4 ‘ 4 

« (2 X ± 2 m ± v) > - -|- , R Q) > o 

R (y) > o, - | < # (v) < i (3*2) 

we obtain (3*1) 

4. Theorem 3 : If ~ * /(x) and <f> (x) belong to L (o, cc) 


and if ^ (j?) cJLs x^ ^ /(x) 
M 


then 


^ (p) = 


<#> (P) 




=• /(*) 

r( i + L+L^ii 1 ) 


r(- + 7* ) p 


X+v + Va 


X J / +3/a s F a £ 


, 1 i v + \ i f 

'> 1 + § . 


] 


4> (y) dy 


3/2, v + 3/2 P 

R (X + v) > | R (*\ | - 2, R(p) > o, R (y) > o, - % < R [v) <.* (4'1) 
Proa/ : Inversion formula ( (6) ] for (l’l) gives 

' ■ ■' (4 -2) 


/(*) = J (*J>0~ # v (*)>) <#> (j) dy 


Substituting this in 

oo 

* (?) = y—- P j K lx (P*) x X ~ 1 /(#) <** 

O 

and interchanging the order of integration as before, we obtain the result (4‘1) on 
using [ (2) p. 165) J 


co 

jV-i (*,)* ^ (<w) H (#,)</* 
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x , s/o r ( 1 + V _±^±J ) 
_ 2 CT w y + 3 / 2 v 2 / 


+ w + 2 


r (v + 3/2) 


1,1 + 


V + or ± f 1 


L 3/2, v + 3/2 J 

ft (tr + M) > | ft (/*) | - 2. (4-3) 

Corollary : Putting /'■ = ± 4 , the theorem reduces to : 

If ' f(x) and <j> (*) belong to L (o, os) 
and if $ (p) f(x) 

i> (P) —M 


t(x) = P X v " '- . v( 1 4 - v -+ * ± i 

v W w r (v + 3/2) \ 2 


x i y + 3/2 3 p a 


i, i + 


v + x ± 4 


<t> U) dy 


3/2, v + 3/2 


ft (X + v) > - 3/2, ft CP) > o, ft (j) > o, - J < ft ( v ) < | 

5. Theorem 4 : If /(V ^ ) an d 4 1 ( x ) belong to L (o, os) 

• - - and if f (/>) X /(^ ~) 

ft, m 


^ (/>) /(*) 


T 1/2 p~ v / 2 - 3 / 4 r(7/4+'v/a) r(7/4 + v/2+2m) 

2 2m + v+l • f fv + 3/2) f(9/4+« - ft- r/2) 


05 r 

jV^L 


1 , -J + ~ , X + y + 2m 

4 4 4 2 


3/2, v + 3/2, 4 + m - k + ~ 
4 Z 


*(v)4r (5*i) 

v' 
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R ( T + ~2 + 2m ) > °’ > °> R (>) > °> - i < ^ (v) < 

Proo/ : Substituting for J(,\t~x) from (4’2) in 

CO 

t (t>) = P J « ’ 1 ^ (/>*)"' _Tr w k> m (px)j(ST) dx 

O 

and inverting the order of integration, we get on using [ (2), p. 171 ] 

03 

j x 2x (x } ) 1 / 2 e~ ] l^ x ' Z W k M (1/2**) 

o 

_ 2 1/4-a-v/ 2 r(7 /4+v/2 + x ± /*) * + 3/2 

r (V + 3/2) r (9/4 + X - fe^Tv/2) J 

p, 7/4 + v/2+ X ± m. „ a ' 

x p* • , 

3 f 3/2, V + 3/2, 9/4 + X - A + v/2 2 

P (2 X + v) > 2 J P (m) | - 
the Theorem (5T) 

Corollary : Putting# = -mf 1/2, the theorem reduces to : 

If /(*/ T) and </>(*) belong to L (o, co) 
and if $(j>) ^ /( a/ - * - ) 

*00 — /(*) 


* (?) = 


1 / 2 J &~ v / 2 ~ 3 / 4 . r (7/4 + v/2) r (7/4 + v/2 + 2m) 


?2m + v + i 


f (v + 3/2) f (7/4 - v/2 + 2m) 


f f 1, 7/4 + v/2 a 1 

X) v + 3/2 8 F a - U(jO<0 (5.3) 

S' J* L 3/2, v + 3/2 J 

— 1/2 < ze (v) < 1/2, P (7/4 + v/2 + 2 m) > o, P (p) > o, P (>) > o 

6* The above theorems may be utilised to evaluate certain infinite integral!. 

By way of illustration we take a few examples. 

Exampte 1 : In theorem 1, taking 

f(X) = X~ 1 / 2 key ( 2 ^«) 
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/(*) = # v+ 3 / 2 



</> (p) is then given by (2), p. 11 & and to find ip (p) we use Cathie’s result 
(5), vis.; 


r (8 )P 


l ~y 


r 08) r (y + m - k + 1 / 2 ) 


+ 2 » !# ] 


= * > ~ 1 ,fP' 

m L 


y 4- m — A: -f- 1/2 

■ : -# 

§ 


£ (r) > o, i? (y + 2 m) > o 
so that we get on a slight simplification after putting 


7 . v 7 , v , ~ n 7 

4 + 2 = a ’4 + ~ + 2m = P; m - k - - + ~ = V 


CO 

/ 




y + 4 a — 5 


1, a, /3 


a" 3 


, . 
^ 2« - 


_ w 1 / 2 r (2a - 2) r (y 4- 3/2) 2 |S + 3a +y - 6 
r («) r (/3) - 


2, y .+ 2a - 2, 3/2 

2* 


_i| 

4/.J 


* v O’) d > 



R (a) > o, R (j8) >o, jR (y 4 «) > 1- (6-3) 

The author is thankful to Dr. P. L. Srivastava for his valuable guidance. 
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ON another representation of an h- function 


By 

V. K. VARMA 
Fergusson College , Poona-4 
[ Received on 29th November, 1962 ] 


ABSTRACT 

Recently Fox has defined an H-function by a Mellin Fame’s type integral, involving product 
of Gamma-functions ; and has proved it to be a symmetrical Fourier Kernel. He further points 
out that it is a most general function and needs further study. In 1958 the author gave the Mellin 
Barne’s type integral representation of this H-function, and established its kernel property basing 
his proof on the L 2 -theory of general transforms. Moreover the author gave for it another interest- 
ing representation in the form of a multiple integral involving product of Bass el’s functions ; thus 
connecting H-function to a well known function of mathematical physics. In a paper he denotes it by 




m 


,X 


(x) where m + n > p. 


A comparision^ of the forms of H and p functions at once leads to the fact of their being identical. 


The $ *» function, mentioned above, is a subsequent developement of a kernel 


v i>' 


\ (x) 


defined earlier. A number of interesting applications, whi eh establish self-reciprocity of well known 
functions have been made by the author. 


1. In a recent paper [ 1 ] Fox has defined a function 


(i-i) 


H 


(#)= 1/2**. J 


<r + tos 


cr - » co 


? P 

it r (bi + cj s) tt r (ij — «ji) **' 

*— i . ; = i 

: is 

q P 

7T r (bi 4-<i — Cis) r r (/3j - «j + «jj) 

i— l j = l 


where a > 0, i - 1, ... q ; <j > 0, j = 1, ...,p ;q > p + 1, and all the poles of 
the integrand in (1*1) are simple. He has proved in details that the function in 
(1 - 1) is a symmetrical Fourier Kernel. He further points out that it is a most 
general function which needs further study. 


The object of this paper is to give for an if- function another representation 
in the form of a multiple-integral involving product of Bessel functions. 


2. In a recent paper, [2 J, I have defined a function 
(2‘1) <p Vl ,.. , v m , \{x) 


- i^rtrr ) <•’' 


1/M 


dy i 

*1— i-i 
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where 

fO o ; Vfr 4* fc > o 5 k « 1 2, m ; A + £ > o. 

I have established that the function defined in (2*1) is a symmetrical Fourier 
Kernel. 

Later on [2], the function in (2*1) has been generalised and the generalised 
function has been defined in the form 


( 2 - 2 ) 


^ v i 


An 


An 


(*) 


V* f«> f * o a) 

R J J 

o o 




/,,(*) - V, (^) 1/Ml ' 3/a / Xl (* i;/ V (^) 1/Md ' _8/2 / Ad (> 1/ ^) X 


(^ 1/M ’ 1 ) ...(ip) l/ ' a ' p-i 7 |p d^dtnui dy^.dy* dt v ..dt B 

where £ — and 

> o , A = 1, 2 n ; /*'r > 0 , f =• l, 2, ... J> ; »» > 1, n andi> are any 

positive integers or zero, provided that m + n > p. An empty product is 
taken as L 

The v’s can be permuted amongst themselves but it is not permissible to 
permute either \’s or |’s amongst themselves or amongst oneanother. 

The Meliin Barne’s type representation of this tunction has been obtained 
[ 2 ], in the form. 

Iii ■"» Ip 


(2.3) 


v \ x (*) 


L C 

2ti J 


c+tc 0 


>(»-*)* 


c — tco 


m r( 

* 4 v 

r= 1 


15L + -i- + 

2 2 4 /1 


|r/ v r + _£ j 

t 1 v ~ 2 ~ 2 4 / 


r 


1 




t = 1 I r 

l 1 \I 2 


At i i 1 


)j 


.. I f T(£ i |^-P■ , is|2+^'iH + ^^2) \ v ~s 
jl 1 \ T (li/2 4- 4 V/2 - f* 'j/4 + 1/2) j 


* 


where A=’jn+^ 1 +...+f*n~/ i, i~ '... — f*'p and 
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m+n>p. The formula in (2.3), unlike the formula in ( 2 . 2 ), is valid when 
m>0. Also n and p may be any positive integers or zero. Aa empty product is 
taken as 1. Taking m= 0, n = q and replacing x by X 2“ 

(2.3) reduces to 


(2.4) 


2 A/2 


£p 

4 Xat 


(* >) 


_1_ 

2m 


c+ico 




P 

7 r 

.7 = 1 




is 


Thus with bi-\il2 - -/4 4 - 1/2, cj=/W/2, * = 1,..., g and tf]=fj/2 + f*'j/4 + l/2, 
«j=p.'j/2, j‘= 1,..., J> ; the function in (2‘4) is identical to H(*) of (1.1). I gave 
this formula as early as 1958 [2]. Likewise (2.3) can be deduced from (1.1) 
by a little adjustment of the values of the parameters. I have established inde- 
pendently [2], using the theory of convergence in mean square sense, that the 
function in (2.2) is a symmetrical Fourier Kernel. 

Further let m= 1, v 1= v, n = q and replace # by 2 h x, where ‘h’ has the value 
assigned above, then from (2.2) and (2.3) we have the equivalent forms. 


(2-5) 


2^/2 £i> • • •> £p» 




vV 2“ 

; R 


/“/“Vft.. hjp’i—yn) J v (2 h * tvhhv-y*) X 




Xft)^-1/* J dy.-.JM dt v ..dh 


1 

2tti 


/ 


C + tco 


c — *co 



J ^ (frj+gjQ 
l r ctf) 


} 


r fa - fif) 1 x ~ s ds 

r fa — e- } 4 Sji) J 

= u (2x)ld~2 



where q in (l.l) has been replaced by (?+l) aQ d we take b in v/2 + 1/4, 
e #+i = 1/2 ; and b h a 
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i = l,...g ; d' h 1,..., p have the values assigned to them in (2.4). Also h 

and R have values previously assigned to them. 

We may also write (2.5) in a slightly different form by applying a result [3| 5 
2.1, 17, p. 52), similar to Parseval theorem on Mellin transforms, in the form 

j ,-A + ioo 

~2 wTJ F (0 G (s) x " s = [ g (u) f(x]u ) dulu 

k — tco J 

where F (*) and G (?) are Mellin transforms of /(#) and j (a;) respectively. 

This gives another formula for the equivalent forms in (2.5) viz. 


( 2 . 6 ) 


>*/2 


f X; 


Ip 


!'• 


(2 h x) 


-r 


H (#/j>) s/^ J (y) dy/y> 


where now H (x) has q factors in the first product as in (1.1) and not J+1 as in 
the defination of H (#) in (2.5), The values of hi , <?$, (7 and flj, P 

are the same as given in (2.4). 

It is believed that this new representation of the H-function provides an 
extention, to its defination given by Fox, by representing it as a multiple integral 
involving the product of Bessel functions. Thus this representation connects H- 
function to a well known function of Mathematical Physics. 
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ON THE ASYMPTOTIC BEHAVIOUR OF SOLUTIONS OF PARTIAL 
DIFFERENTIAL EQUATIONS OF HYPERBOLIC TYPE 

By 

M. RAMA MOHANA RAO and B. VISWANATHAM 
Department of Mathematics, Osmania University, Hyderabad-7 
[ Received on 17 th December, 1962 ] 

ABSTRACT 

In this note we shall derive some criteria for the boundedness of solutions of partial differential 
equations of hyperbolic type of the form 

Zxy =/(*. y> z ’P>$ 

We shall also consider when a solution of the differential equation is unbounded. 

r 1 In this note we shall derive some criteria for the boundedness of solutions 
of partial differential equations. We shall also consider when a solution of the 

differential equation is unbounded. 

In § 2 we shall derive a lemma which is useful for subsequent discussion. In 
Owe shall find the estimate for the bound to the difference between any two 
solutions of the partial differential equation. In § 4 we consider a comparison 
theorem which follows naturally form the lemma proved in § 2. 

Consider the initial value problem 

(A) z*y = / (*» y> Z ’P> ?) z ( x > °) = <T W > 2 

where <r(o) = T(o) = Zo in tlie re gi° n R:o<x<a=,o<JV< 03 • 

Suppose f(x 1 y,Z,p,q) is a real valued and continuous function defi 
o s? * <« , o s£ y <co , - o o < Z,p,q <Tca . 

By a solution of (A), we mean a real valued continuous function s(x,y)hav' g 
partial derivatives Z x (x,y), Zy (x, y), Zxy {x,y) in the region R. z(x,y) 
integral equation 

(1.1) 2 (x,y) = o(x) + T (y) - Z 0 + J f /(*• ^ ^ Zx ^ ^ ^ ***** 

O O 

We shall say that the function <p (*, y, Z, p, q) possesses the property ( ) 

tl p (x,y Z, p, q ) is a continuous, non-negative and bounded function defined 
for (x, y) e R and non-negative z, p, q » non-decreasing m each of t e varia es Z, , q 
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before proceeding further We shall derive a iemrfia 

§ 2. Lemma : Let (i) a (x % y ), /3 (x $ y), y (#, y) be non-negative >' o, measurable 
functions defined in R such that a is continuous, (3 is nniformely Lipschiz’s continuous 
with respect to y and y is uniformely Lipschiz’s continuous with respect to x and 
satisfy 

x .y 

(2.1) a (* s j)) ^ J J f (s, t, a (s, t), /3 ( s , t), y (j, t) ) dsdt 

O O 

r y 

(2.2) p (*, y) < J f (*. <» a (*> 0. v (*, t) ) dt 

o 

(2.3) y{x,y) < I ? (s, « (*, y), P (*. J'). V (*, j) ) * 

o 

(ii) <f (x, y, Z,p, q) have the property (*) and * (x,y) be the mazimal solution of 


IB) 

then 



t, 4* (s, t), fx (*. t)> 1>J ( s > *)) dsdt 


a (x, y) x (*> y) 


P l*> y) =s£ Xx (*> y) 
y ( x > y) < Xv (*> y) 


Proof: 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 


Define the sequence of successive approzimations to (B) as follows: Let 
2 ° (x,y) = a (x,y), u° (x,y) = p (x,y), »° (x,y) = 7 (#, y) and for k >1 
„y 


Z* (*, 7) *£ J" j~ f {s, t, 2 k - 1 (s, f), u k 1 (f, t), v* 1 ( s , t ) 

O O 

a k (x,y) ^ <P (x } t, £ k_1 (x, t), a 15 ' 1 (s, t), o' 1 ' 1 (*, t) ) 

O 

/ X 

t ( s,y , ■s’ 1-1 y), m ' 1 " 1 («, y), o' 1 ' 1 ( s > y ) ) 


dsdt 


dt 


ds 


The functions z k , u k } v k are defined on R which can be taken independent of L 
The inequalities (2.1), (2.2), (2.3) gives the cases ft = o of 
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(C\ n\ k f" 1 ^ . ^-[-1 k ll-l-i ’ •» 

(*.oj £ ^ 2 , a ^ u , z; ^ z; 1 

The cases when fc>o follows by induction by virtue of monotony of 9, The 
boundedness of 9 implies the uniform boundedness of the functions 

k k k 

Z a a , v hence -as- &-»co. 

k k k 

(2.9) Z — Ztm £ , u, = Ziff* & > ff = /tm p exists in R, It is obvious from (2.4), 

(2.8) and (2.9) that - , . . ; „ ; ; 

(2.10) o<a^,e:,o< (3 ^ u , o 

Lebesque’s theorem on term by term integration, under bounded convergence 


implies 


(2.11) 

r x py 

2 (*. y) = J J 9 (1, t, 2 (f } .4 a (s, t), v (*, t) ) dsdf 

O O 

(2.12) 

« (*>y) = J 9 (*, t, Z (X, t), U {X, t !), P (*, <)) 

0 

(2.13) 

c x 

V ( x > y) — J 9 (s, y, z (r, j), a (t, j>), 0 (s, y)) ds - 


It is clear that Z x = u > Z r = v almost everywhere. Thus the lemma follows from 
(2.10) and with the assumption (ii) of the hypothesis. 

§ 3. Theorem 1. Let (i) 9 (x,y, * y ) have the property (*). 

(ii) the function /of the differential equation (A) satisfy the condition 
I /(*> y, ^ 1 , u\ p 1 ) - f(x, y 3 Z\ u\ o a ) | f (Z, y, | ] ; | u i_ a 2 | f 

| p x -p* j) 

where the functions i 1 , 2% a 1 , a», pi and p 2 being arbitrary; 

(ii*) X (*» j) be the maximal solution of (B), then the difference between the 
two solutions of (A) together with its partial derivatives w. r. t. x and y are ^ x 
( x ,y), X* (#> j)> Xr (*» y) respectively. 

Proof-. Let 2 1 (2,.y), 2 2 (x,y) be two different solutions of (A). 
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Let (x,jy), v 1 {x,y) <uid a* (x,y), v* (x,y) be the functions associated with 
Z l and £ a . We have 

«x r y 


I (*> jO - (*,jO I < f j I / (*, *> *\ - f(s 9 h z\ u\ v ■) | 


J j y 9 (s, t, I Z 1 - Z 2 I , | u 1 - u* | , | v'-D* I ) 


dsdt 


dsdt 


putting 


a = | Z x -& | 

/3 a | a 1 — u* | , 7 = | 0 1 — fl 2 | 


we obtain 

x y 

| | J J 9 (*, <, «, j9, y) dsdt 

o o 

Similarly we get 

/0> 

I «,*-«,> J <J 9 (#> *. <*>P>y) dt 

o 

| V —2 y a | s£ J 9(s,y,y, p,V) ds 

O 

By applying the lemma, we get 

| | 

| Zx 1 — 2 x a [ ^ X* (*» j) 

I Z^-Zj 2 I ^X/(*>^) 


This completes the proof of the theorem. 

Corollary /. If all the functions involved are defined throughout #j-plane and 
if x (x,y) } X* (#> y) Xy (#> f) are bounded as (#, y) -> co if any one solution of the 

differential equation (A) together with its partial derivatives w. r. t. x and^ are 
known to be bounded, then every other solution together with its partial 
derivatives is also bounded. 

Cor : 2. If in the above corollary we suppose that / (x, y, o, o, o) = o then 
every solution together with its partial derivatives w. r. t. x and y is bounded. 

Theorem 2: Let (i) <p (x f y, f t \ jt x , ip y ) have the property (*) and 9 (x, y, $%, $ j ) 
= o if and only if </> = o r = o, $ y ss. o. 
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(ii) the function/ of the differential equation (A) satisfy the condition 
I y (x,y, Z 1 , u\ v l ) - / (x, y, u* s v*) | > y (x,y, | z 1 - Z 2 I , | uf - u? | | v l - v* j ) 

Then if any one solution of (B) together with its partial derivatives w. r. t. x and^ 
is unbounded as (x, y) — > oo , then at least one solution together with its partial 
derivatives w. r. t. x andj> is unbounded. 

Prooj: Let Z\ Z? be any solutions of (A). Let u l , u 1 and v 1 be the functions 

associated with Z 1 and Z* 

z\ y -&x y - fix, y, z\ « J , p 1 ) - fix, y, z\ u 2 , a 2 ), 

putting 

a = Z 1 - Z 2 

P = a 1 — «*, V = v 1 - n 2 
Therefore 

Z 1 xy — 2*iy = fix, y, a + Z 2 , p + a 4 , y + v 2 ) - f{x,y, Z\ a 2 , » 2 ) 

= F (x, j>, a, p, y) (say) 

Hence by condition (ii) 

| F ix.y, a, P, y) | > f ix, y, | a | , | 0 | , | 7 | ) 

Since F ( x,y , «, jS, y) is continuous and <f (x,y, a, P, y) > o, 
the above inequality shows that o, p, y are different from aero, 

Fix,y, a, p,7) ^ o which implies that either 

F (x, y, a, p, y) > f i x ,y, l“i> \ P \ , I Y I ) 
or F {x,y, <X,P, y) « - f ix,y, |«|,lP|.|v|) 

Let us consider the first case: 

By similer argument of lemma in § 2 it can be easily verified that 

« i x ,y) > k ix,y) 

P ix,y) > K (x,y) 
y ix,y) > k T (x,y) 

where k (x, y) is the given unbounded solution of (B). Hence the results follows 
immediately. Similarly we can prove the other case. 

§ 4. In this section we shall compare the solutions of two differential equations. 

Theorem 8. Suppose 

Zxy — f ix,y, Z,p, q '); Z (x, 0 ) = o(x), z ( o,y ) = r (y) 
zxy"= & i x ,y, z, <f, q)> z ix, ») -- v(x), « (o,y) = f iy) 
where o(o) = T(o) - ^ are two differential equations in R, satisfying the follows 
conditions 
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( i ) f ( x > y > P> q) t s i x > y> z, p» 9 ) I 

' < ? ■(#» js I I > I p~ J, I ,-| 2 - 71) 5 

00 V (#, y> &> P> 2) have the property (*) and x ( x >y) * s maximal 
lolution oi (B) and x (#> JV)> X* i x >y)> Xy (#,jv) are bounded as co 

(iii) Any solution of one of the equations together with its partial 
derivatives w. r. t. x and y is bounded then, any solution of the other equation 
together with its partial derivatives w. r. t. x and y is also bounded. 


Prooj : The solutions of the two equations are given by 

r y 


z (#» y) = CT (*) + r iy) - Zo + j J f (s,t, z, P, q ) dsdt 

o o 

X y - 

* <*»j 0 - <K*) + r O) - z 0 + J J g 0, t , z p, q) 


dsdt 


Therefore 


© o 


r x rJ 


I 41 “ « I < J j I / ( s > i»Z,p, q) —g (s, t, z> p, ~ q ) | 


dsdt 


x y 

J j V I Z- z | > I "P~p I * | q 1 ) 


dsdt 


Putting 


* - \ *-z \ > P = \ P- p\ >7 = \ q-q \ 

x y 


dsdt _ 


1 *- 1 1 < / J ^ (*»*»■«. a.?) 

o o 

Similarly l P~ P I ^ j f .(*> t > a > P> V) t// 

O 

I 2 " 9 I *S. J f ( s > J. «. / 3 . v) rfs 


By applying the lemma we get 

I « - * I < X (*, J>) “ - ■ - 

| - 2x | ^ Xx (*»y) ■ • • - ■ 1 • 

I ~ I .*<Xr(*.jO - 

„ and the desired result immediately follows from (ii) and (iii). 

^ Corollary : If g (x y y> o t o, 0 ) — o and x> X*> Xy are bounded as (x, y)-> co every 
solution of (A) with its partial derivatives w. r. t. x and y is bounded as (#, jy) co . 
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ABSTRACT 

The author here establishes two theorems, in which various integral transform are connected 
through a chain of relations, after first proving two lemmas required in their development. Obtaining 
some interesting properties of the MacRobert E-Function, he utilizes them to evaluate certain infinite 
integrals involving that function. 

1. Introduction : In this paper we deal with certain properties of 
MacRobert E-Kunction and utilize them to develop a chain of relations under 
various Integral Transforms. Some new integrals involving E-Function have also 
been evaluated. 


Meijer (4) gave a generalization of the classical Laplace Transform 
co 

* (J>) = P j e~ pt f(t) dt, R (p) > o 


as 


oo (pt) f (t) dt 

4>(J>) = P J {pl)~ k ~ h * 


R (?) > ° 


while Vartna (8) generalize it in the form 

05 (/> 0/(0 dt 


(LI) 

( 1 . 2 ) 

(1.3) 


0»>— * , R (f)>0 

o 

The generalized Stieltjes transform of/ ( t ) of order \ is given by 

CO 

<t>(p) = P j (P+t)~ X fit) dt i arg /> I < vr (1.-4) 

O 

Throughout this note we shall repiesent (1.1). (1*2). (1-3) and (1.4) symbolically as 
fit) = <p(p) 

M 

f(t) > </> (p) 

kjn 

V 

f(t)—><h( P ) 

km 


respectively. 


and fit) > <j> ( p ) 

A 
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2. Before establishing the main theorems we shall prove the following 
lemmas : 

Lemma 1. If g (t) belongs to L (0 S co) and 
M 

if *(*)_> f(j>) ... (i) 

k,m 

/( 0 = ••• 00 

and < V “ 2 <#> (t) = Xi (/>) ••• (iii) 


then 


X lW - f 1 -’'/ «- 2 E [ 2- *+rj t 2 T :#»] « («) * (2.1) 


provided R(p) > o, R (2 -k±m) > o and the integrals involved aie convergent.' 
Proof : With conditions (i) and (ii) above Jaiswal (3) has shown that 

* (?) = ^ty> P J 2 A [ : - V ] *» 4 ” 

© 

provided R (p) > o, R (2 - &±m) > o and the integral on the right is convergent. 
Finding <j> (0 and substituting its value in 

CO 

\l (p)—P f e ~^ i7 ' 2 0 (0 dt 

o 

and interchanging the order of integration which is permissible under condition* 
stated (Fubini’s theorem) we obtain 


X, M - > J ** . <* J ■-« ‘ y -' A [ ■ - f] 

Evaluating the latter integral by (2) p. 212, we obtain the result. 

Lemma 2. : If f a g (0 belongs to L (o, co ) 

M 

and if g(t) — > f(p) 
k>m 

f co = </> (/>) 

^ <f> (t) ~ X x (?) 

M 


co 


dt du 


then 


and A-i (~) > X 2 (P) 

M/ K ,» l 

OS 


x, (?) - p 1 - 7 -’- / ] * «* 
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( 2 . 2 ) 


R (Y) > o, R (2 -k±m) > o, R (Y+Yj-i, ±w x ) ;> o 
Proof : Using the result of Lemma 1, we substitute the value of £i \ x (fit) in 

O 

and interchange the order of integration which is permissible as before* to obtain 

X* (P) = P J «- 2 g (») j i 7+Yl ‘ 1 W^<f?> 

o o 

X E [ 2 — 2A Zm * y : " u lt ] dt. du. 
We can evaluate the latter integral with the help of (5, p. 171) 

du 


C e -ium+ y-% (u) E r«i ; 1 1, 

J k ’ m L A-./G* “ J 

/3 j . . . . y+m — J 


(A) 

“ /->J /->s* / -f'7/* — /h-J— 5 «■» 

and obtain (2.2) 

3, By repeated application of result (A) we are in a position to generalize 
the result obtained in (2.2) and to state it in the form given below. 


Theorem 1 : If t~ * g (t) belongs to L (o, co) 

M 

g (t) — >f(p) 
k t m 

/(0 = * (P) 


and if 


/Y-4 ^ (0 = X, (p) 


t 7 9 \ 


*• 0 ) 

1 


M 

k l m l 

M 


X ip) 


(~) — > x 3 (p) 
1 k%, m s 


VT 


* n 2 A.c-2 ( '7" ) ^ 

^n-2> 
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and t 


Vn-i 


^ ( 7 ) 


u 


then 


kn- 


Aa ( P ) 


1> ^n-l 


X„ ( P ) = 

x f «-2 e f y > y J r a r - k r + m f , y -f a f 

V L 2 - '2k, y + a r -2 k r 

r 

where a f = S lf r . r = 1, 2, ...(«- 1) and R (7) > o, 


/ Wf : />“ ] g («) 

(3.1) 


R (2-A±m) > o, R ( 7 + a r - k f ± m f ) > o, r = 1, 2, ... (n - 1) 
Corollary : Taking k — ± m, k f = the theorem reduces to : 

L (o, co) 

and g(t)=f (p) 

( P ) 

t y ' 2 <j>(t) = X 1 (p) 


1 7l *x (7) = Xa(« 


<y “ 1 Xo_i (7-) = \ n (/>) 

then 

1 03 

An (/>) = P 1 J* u ~* E £ 2, y, y -| - a y : : J g (u) du 

o 

provided 


R (?) > °» R ( V + « r ) > o, a = 2 y 

1 

Example : Taking g (/) = - I £ a > b 

its Meijer’s transform is given by Verma (9, p. 172) as 


; r = 1,2, ... (« - 1 ) 


(3.0) 


f{p) = Ap- l ~ k ~i 
R ( P ) > R (c-a-b) >0, 



1 -8 -a, 1 -8-b, l-k+\ 

*"5j r a b ~ ft, l -\-m — / — m T ^ 


) 


and A = JXf) 

T(a) r(c-a) f(e-b) f(b) 
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R ( Z±m 4-0 +* S+i ) > o, R ( /±n*+^5 + \ ) > o 


The Laplace transform is then (2, p. 222) 


:7J 


t \ A Z+fc _1_ 1-8 — 0, 1 — S — ^5 l-k+.k> Z+4+i \ * *"** 

* (p)=Ap - 8 ,«-.-i- 8 ,i + »+i.l'-«+i / 


R ( l+k+S - 4) > °> R (l+k+a+b 8— « - 4 ) > o, | arg^> | < 5w/2 

l-f-8, 1 -c+a +6+8, | - / 4- m 

S+-0, 5 — J — £►, 4“"Z”+&, 4 l - k 


= A // + & + i G 


3, 4 

4, 4 


( 


) 


since (1, p. 209) 

G W ’ n 
P> 7 




(B) 


y'2 


Finding L. T. of r <j> (t) we get 
Xi(P) = a/ /8-/ -*- Y 


G 3 ’ 5 I L 3 /a -l-k-y, 1 + 8, 1 -«+«+'* +s> i~l+m 

5 > 4 \P 8 + a, 8+6, % -i+fc, i-l -k ■ 


) 


R(y+2*)<o, R Q-i-k-y -«-8)>o, R Q — l—k-y—b-ti) >o 

Now the Meijer’s transform of i Yl X x (-) can be found by using Verroa’s result 


(9, p. 173), so that we have . * 

X 2 (p) = A/ /a - N *~ Y “ Yl 

,3,7 / 1 I y, 1+8, 1 - c + a +6+8,|-Z + m 


xG_ 


' 7 > 5 ^|s+a, 8 + 6, i-i + fc, 4 -z-fe, 3 / 2 -Z-fc-y-Y,+?£i 

R (r+y 1 + 2A:-* 1 ±mi)>o, R (S + a+Z +*-+Y+-yi-fc,±»»j - 4)>o 
Proceeding similarly we get ultimately, 

An 0) = A //“ l-k-y-a n _. 


) 


p3, 2/2+3 
^2 /j + 3, /z + 3 


( 


y 2 -l~k-y-a r 4- fc f ±m r , 3/ a -Z-fe-y, 1+5, 
1 - c+a + A+S, | - l ~ m, \ - 1+ m 
8 +a, 8 + 6, 4 — 3 /j — l—k — y - a r ■+ 2 fcj, 

[ 289 ] 


) 



R(7+2 fc+« r -fc f ±m r ) > o 


* v 

Where « r - 25 Y r , r *» 1, 2, 3 ... (fl — 1) 

R (5 + ^ + /-^+r + a f -fc f ±m r -1) > o, 
r =* 1» 2, ... (« — 1) 

Substituting the value of \ n (p) ia theorem (3.1) we get after a little 
simplification and putting 

y ”j" *■* kf +- TWj* — 07* 
y + « r — k r - m r = b T 


co 

J 


j& + J + 8~ D /a 


y + tf^ -r 2 &t* — 

E r 2 -A±m, y, a 1 ,...a w _ u *1. •••*„_ i 

2 ~ 2k > e x> c n _ i 



IM _ x 
“ r («>■ r (« — r <«— *) r 

?*, 2n+3 £ | I 3 /a ~ l — k~ a r , ^j 2 - l ~ k — by 1 ~f*8j 

2n f 3, «+3 It 1 — c + a + 6+8, £ — J±m 

^ I 8 + fl, 8 + &, J — liik, — 

(3.3) 

r == 1, 2, 3...(«-l) 

R (p) >o, R (c~ a — b) >o, R (/±m +0 + S + i) >°> R (i±m + 6 + 3+£) >o 
R (y + 2 &) <Co, R (8 4 " 0 'M+&+^ — c — ^) o, R (|-Y«/ - fc -8 — 0 ) ^ o 
R (4 — / — fc— Y -8 — ft) >0, R ( 8 + 0 -W+& + 0 f ~ 2 )>°> R (8 + 0 + i + fc + ^ f - i)>° 
R (a, +2A) >o, R ( 6 f +2fc) > o 

4. Theorem 2 : If i~ y j> (f) belongs to L (o, co ) 


and if 


4> (0 

/(0 


V 




\ 


P 2 - y JW 

HP) 
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* y > 


t 


>'2 


Hr) 

*.(4) 




(?) 


^*2 ( P ) 


V 


f 7 ”- 1 X ( T ) * X “‘ l W 

n -*i \ t / U 77 ? — , 


^n-is 

V 


md 


i 7n X (I) » x «0>) 


then 

1 _ x _ „ n 1 f - 7 F r X, 7, 7+2 m, X+«r> X-H*r+2»r . ^ ( u ) d 

Ab(/>)=P . ”r7\ x I u E L y+m-fc+l> ^4- a »- +wr-^r+a J 

^ ; o (4.1) 

where a. = 2 7r- r«= 1, 2,...n. and R (x) > 0 ,. R 00 >° 

1 

R (y+2m) >o, R (X+«r) >°, R 0+«r+W >0 » 

p ro nf : With the first three conditions above Saxena (7, p. 349) has shown that 

X (p s , f a -7 e[ X ’ 7 ’ y+2m :^« ] <^» («) riu 

r(X) J 1 y+m-k+i J 


r (x) >0, R (7) >0, R (y+2m) >0 

Finding t 7l <P ( -y ) and substituting its value in 


x.M = p ]'- lpt *<4 )"' 

o 

we get on changing the order of integration and using the result (A) above 

p 1 ~ y *~ x r - r E r x, r, y+ 2w > x+7 a , x+y 1 + 2 >ni ipu j ^ ( u ) ^ 
Xi.( 6 ) T (X) J *■ y-t-m-4-Hr, X+7i+%~ ^i+i 

R (X) >o, R (7) > 0 , R (7+2m) >o, R (X+7i) >0 R (X+7i+2wi) >*> 
Proceeding similarly we get finally, 
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An (fi) 


where 


. I - X - a a o® g- 

f Tw~~ J “~ r E [ 


‘A, 7} y 2nt ) X -{-* at ,\-ba: *4*2w2 


7 + tn-k+%,\+a r + mr 


x -f-zmr -j 

- k+i 1 p u J 4 * (“) du 


r 

a r 5=8 ^ 7 r ; 7==1 > 2 ‘ > ‘ n * and . R (A) >o, R (y; >o, R >o 

^ R c>, ■ R (x^f***!* -J-»2 m r ) ]>o 

Hence the theorem. 

Corollary : Taking k f =s — m f + \ ; r = 1, 2...« the theorem reduces to ; 


If t ^ (t) belongs to L (o, co) 
V 2^v 

and if <j> (t) > p l J f ( ; p ) 

fc, 7/1 


/ (0 — > * (/>) 

\ 


•Sir? 


then 


‘ ' ‘ * 7 ’ K T ) # Al (/ ° 

^ ( -J ) =r X a (/>) 

i 7n Xn ~ 1 ( y ) = Xo (p) 

Xn (?) = ~ ~ f u“ 7 E [ X ’ 7 ’ 7 + 2m > X +“r :pu 1 <J» (u) du. 

i W ^ L y+m _ fc+% J 


(4.2) 


where « f =U f .f= 1, 2,...n ; R (7+2m) > o and R (x + « r ) > o 

R (X) > o, L< (7) > 0 
Example : Take 

$(t)= a v iC -2m-l F r 4-viA 1 . , 


? r 

5 L «, r — ‘ 


2m 


: - 1 ] 

0 J 


r («) r (c — 2 m> 

so that (7, p. 349) 

1-f-c+v — 2m — 7, 1+v 

X+e+ v “ y ~ + i- V- 

R (X) > o, R (7 +2m — c) >o, R (7 + 2m-c-f|-v±/*) >o 


J2m+y — c— p - X+l 3 o 

Jl p 5 > ^ 

W f (X) T(i-v±^ U 2,3 


( 




) 
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Therefore 


t Ti *(7 ) 


1 ^ 2m — y+Vi — 1 

T ' = r ( a) r 

x G 3 ’ 2 ( «| 1 +<-'2m+v — r, 1+v > 

2, 3 1 ^ 

\ I X+ v +e- 2 m-r, i+/t, 

f \+v+<-2m-y+y 1 - l 

r (\) r ( 4 — vi^t) 


2,3/ l-x-v-c+2m+y, |-,t, 

3 . 2 I a< 

\ y + 2 m-c-v, —v 


on using (B). 


Finding its Varma-transform (6 p. 275), we get 

y - \-v - c +2m+y - y! 

Xi (J>) = foot ft-vi#*) 

2 5/J 1 - X- v-c+2m+y-y lJ 1- \-v~c+2m+Y- yj-2%, 
X G 5* 3\J>\ 1- A- v -*+2m+Y, 1-P. 

\ I Y+2m — c — v, —v, 1 - \ - v — c 4- 2m +* y — Yj+&i - — J 

I ar ? ^ I < ^ 2 ' 7r i (A. + ?i) > o, R (\ -f Yi+c — 2m - y+mj +mj) > o 


Proceeding similarly we obtain 


An (p) “ — 


l_X~v-c4-2m + y~« Q 

~r“5Tf(i-v±'#i) 


^2,2 «+3 / | 1 -\+Y ~c + 2m — v — 1 - X+Y — c+2m- v-a^-2m r * 

2 tz+ 3 , n f 2 1 + 1 -\+Y - c+2m — v, J — /*, I +/ 4 \ 

\r I y + 2m - c« v t —v, 1 — \+y + 2m- v~o^ + ^ ~m^ — J t 


where, « » % y ; r=I, 2 ,.. t «. 

r 1 

Substituting in the theorem we get after a little simplification and putting 
i " v "f J 4 — a \* i — v — /*= a 2 , i+c — k — m = j 3 
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c=b 1 ,e-2m = b^\+ot r = « r ; XH *. + 2 ro f = d f and \+« f - - k f + 1= e f 


/ 


A 


— 




°1> a i> a Z 

_ u 


^‘2* C f * 

: 

E 

f \pu 

bi b. 

a 


y+fl 8 - f> a , ^ 




L- 


... r_W r(M a K^i+^-i) *(«,+,*- i)+y-i, 

r (aj) r foo r («,) F 


|(<Ji+0 s 4-l)+7 -6 2 -e f »%(a,+^a+^) + ' y — v 

i(^i+ flj+ i(l — ^(1 d" a i — °s) \ (4.3) 

|(a, + <j 4 - l) + 7 — t 2 , + 0> K a i+ fl 2-H) / 

+ y-b i -e r ' 

where a =J v . r = l,2,3...n. 
r 1 r ’ 

R (X) >o, R (\ - i 4 ) > o, R (7) > o, R (7-6 2 ) > o 

R (a 1 -i a + y) > o, R (a 2 -6 a +7) > o, R (e r ) > o, R ( d r ) > o 

R (i a +C f -7) > o, R (VM f - 7) >o and 1 arg> | < 2 tt 

The author is thankful to Dr. V. K. Varma for suggesting the problem and to 
Dr. P. L. Srivastava for his valuable guidance. 


x g: 


2,2n+3 
2n + 3, n + 2 
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